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ABSTRACT 


Stability  experiments  are  conducted  in  the  Arizona  State  University  Unsteady 
Wind  Tunnel  on  a 45°  swept  airfoil.  The  pressure  gradient  is  designed  to  pro- 
vide purely  crossflow-dominated  transition;  that  is,  the  boundary  layer  is  subcriti- 
cal  to  Tollmien-Schlichting  disturbances.  The  airfoil  surface  is  hand  polished  to  a 
0.25  pm  rms  finish.  Under  these  conditions,  stationary  crossflow  disturbances  grow 
to  nonuniform  amplitude  due  to  submicron  surface  irregularities  near  the  leading 
edge.  Uniform  stationary  crossflow  waves  are  produced  by  controlling  the  initial  con- 
ditions with  spanwise  arrays  of  micron-sized  roughness  elements  near  the  attachment 
line.  Hot-wire  measurements  provide  detailed  maps  of  the  crossflow  wave  structure, 
and  accurate  spectral  decompositions  isolate  individual-mode  growth  rates  for  the 
fundamental  and  harmonic  disturbances.  Roughness  spacing,  roughness  height,  and 
Reynolds  number  are  varied  to  investigate  the  growth  of  all  amplified  wavelengths. 
The  measurements  show  early  nonlinear  mode  interaction  causing  amplitude  satura- 
tion well  before  transition.  Comparisons  with  nonlinear  parabolized  stability  equa- 
tions calculations  show  excellent  agreement  in  both  the  disturbance  amplitude  and 
the  mode-shape  profiles. 
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0.9  mm.  Rec  — 2.4  x 10°,  [6|  1 2]  roughness 239 

5.90  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.45,  Y — 

0.9  mm.  Rec  = 2.4  x 106,  [6 1 1 2]  roughness 240 

5.91  Individual-mode  disturbance  amplitude  and  amplification  factor  N 

for  Rec  = 2.4  x 106  and  [6|  12]  roughness 241 

5.92  Total  and  individual-mode  disturbance  amplitude  and  amplification 

factor  N for  Rec  = 2.4  x 10(l  and  [6 1 1 2]  roughness 242 

6.1  Streamwise  velocity  contours  at  x/c  — 0.05.  Rec  = 1.6  x 106,  [6|  12] 

roughness 243 

6.2  Streamwise  velocity  contours  at  x/c  — 0.10.  Rec  = 1.6  x 10fi,  [6 1 1 2] 

roughness 244 

6.3  Streamwise  velocity  contours  at  x/c  = 0.15.  Re,  = 1.6  x 10b,  [6 1 1 2] 

roughness 245 

6.4  Streamwise  velocity  contours  at  x/c  = 0.20.  Rec  = 1.6  x 106,  [6 1 12] 

roughness 246 

6.5  Streamwise  velocity  contours  at  x/c  = 0.25.  Rec  = 1.6  x 106,  [6 1 1 2] 

roughness 247 

6.6  Streamwise  velocity  contours  at  x/c  = 0.30.  Rer  = 1.6  x 106,  [6|  12] 

roughness 248 

6.7  Streamwise  velocity  contours  at  x/c  — 0.35.  Rec  — 1.6  x 106,  [6|  12] 

roughness 249 

6.8  Streamwise  velocity  contours  at  x/c  = 0.40.  Rec  = 1.6  x 106,  [6|12] 

roughness 250 

6.9  Streamwise  velocity  contours  at  x/c  = 0.45.  Rec  = 1.6  x 106,  [6|  12] 

roughness 251 

6.10  Streamwise  velocity  contours  at  x/c  = 0.50.  Rec  = 1.6  x 106,  [6(12] 

roughness 252 
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6.11  Streamwise  velocity  contours  at  x/c  = 0.55.  Rec  — 1.6  x 10°,  [6|  12] 

roughness 253 

6.12  Streamwise  velocity  contours  at  x/c  = 0.60.  Rec  = 1.6  x 106,  [6 1 12] 

roughness 254 

6.13  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  --  0.05.  Rec  — 1.6  x 106,  [6|  12]  roughness 255 

6.14  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.10.  Rec  — 1.6  x 106,  [6 1 1 2]  roughness 256 

6.15  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.15.  Rec  — 1.6  x 106,  [6|  12]  roughness 257 

6.16  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.20.  Rec  = 1.6  x 106,  [6 1 1 2]  roughness 258 

6.17  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  — 0.25.  Rec  = 1.6  x 106.  [6|12]  roughness 259 

6.18  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.30.  Rec  = 1.6  x 106,  [6|  12]  roughness 260 

6.19  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at,  x/c  = 0.35.  Rec  = 1.6  x 10°,  [6 1 12]  roughness 261 

6.20  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.40.  Rec  = 1.6  x 10(\  [6 1 12]  roughness 262 

6.21  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.45.  Rer  = 1.6  x 10G,  [6 1 1 2]  roughness 263 

6.22  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  - 0.50.  Rec  — 1.6  x 10fi,  [6 1 12]  roughness 264 

6.23  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  — 0.55.  Rec  — 1.6  x 106,  [6|  12]  roughness 265 

6.24  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.60.  Rec  = 1.6  x 1 0'\  [6|  1 2]  roughness 266 

6.25  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  — 0.05.  Rcr  = 1.6  x 10°,  [6|  12]  roughness 267 
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6.26  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  = 0.10.  Rec  = 1.6  x 106,  [6|  12]  roughness 268 

6.27  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  — 0.15.  Rec  = 1.6  x 106,  [6 1 1 2]  roughness 269 

6.28  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  = 0.20.  Rec  = 1.6  x 106,  [6 1 1 2]  roughness 270 

6.29  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  = 0.25.  Rec  = 1.6  x 106,  [6 1 1 2]  roughness 271 

6.30  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  = 0.30.  Rec  = 1.6  x 106,  [6 j 12]  roughness 272 

6.31  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  = 0.35.  Rec  = 1.6  x 106,  [6|  1 2]  roughness 273 

6.32  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  = 0.40.  Rec  = 1.6  x 106,  [6|  1 2]  roughness 274 

6.33  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  = 0.45.  Rec  = 1.6  x 106,  [6|  1 2]  roughness 275 

6.34  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  = 0.50.  Rec  = 1.6  x 106,  [6 1 1 2]  roughness 276 

6.35  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  = 0.55.  Rec  = 1.6  x 106,  [6|  12]  roughness 277 

6.36  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  xjc  = 0.60.  Rec  = 1.6  x 106,  [6|  12]  roughness 278 

6.37  Disturbance  velocity  contours  at  xjc  = 0.15.  Rec  — 1.6  x 10°,  [6 1 1 2] 

roughness 279 

6.38  Disturbance  velocity  contours  at  xjc,  = 0.20.  Rec  = 1.6  x 10(>,  [6 1 12] 

roughness 280 

6.39  Stationary  crossflow  mode  shapes  for  Rec  = 1.6  x 10fi  and  [6|  12]  roughness.  281 

6.40  Total  disturbance  amplitude  and  amplification  factor  N for  Re,.  = 

1.6  x 106  and  [6 j 12]  roughness 282 

6.41  Streamwise  velocity  contours  at  xjc  = 0.05.  Rer  = 3.2  x 10(),  [6 1 1 2] 

roughness 283 
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6.42  Streamwise  velocity  contours  at  x/c  = 0.10.  Rec  = 3.2  x 10fi,  [6|  12] 

roughness 284 

6.43  Streamwise  velocity  contours  at  x/c  = 0.15.  Rec  — 3.2  x 106,  [6 1 1 2] 

roughness 285 

6.44  Streamwise  velocity  contours  at  x/c  = 0.20.  Rec  = 3.2  x 106,  [6)12] 

roughness 286 

6.45  Streamwise  velocity  contours  at  x/c  = 0.25.  Rec  = 3.2  x 106,  [6|  12] 

roughness 287 

6.46  Streamwise  velocity  contours  at  x/c  — 0.29.  Rer  — 3.2  x 106,  [6|  12] 

roughness 288 

6.47  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.05.  Rec  = 3.2  x 106,  [6 1 1 2]  roughness 289 

6.48  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.10.  Rec  = 3.2  x 106,  [6|  12]  roughness 290 

6.49  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.15.  Rec  = 3.2  x 106,  [6)12]  roughness 291 

6.50  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.20.  Rec  = 3.2  x 106,  [6 1 1 2]  roughness 292 

6.51  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.25.  Rec  — 3.2  x 106,  [6|  12]  roughness 293 

6.52  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.29.  R.ec  = 3.2  x 106,  [6|  12]  roughness 294 

6.53  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  = 0.05.  Rec  = 3.2  x 106,  [6 1 12]  roughness 295 

6.54  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  = 0.10.  Rec  = 3.2  x 106,  [6(12]  roughness 296 

6.55  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  = 0.15.  Rec  = 3.2  x 106,  [6|  12]  roughness 297 

6.56  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  = 0.20.  Rec  = 3.2  x 10°,  [6 1 12]  roughness 298 
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6.57  Spamvise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 


at  x/c  = 0.25.  Rec  = 3.2  x 10° , [6|  12]  roughness 299 

6.58  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  — 0.29.  Rec  = 3.2  x 106,  [6 1 1 2]  roughness 300 

6.59  Stationary  crossflow  mode  shapes  for  Rec  — 3.2  x 106  and  [6|  12]  roughness.  301 

6.60  Total  disturbance  amplitude  and  amplification  factor  N for  Rec  = 

3.2  x 106  and  [6|  12]  roughness 302 

6.61  Spanwise  hot-wire  scan  at  x/c  = 0.05,  Y = 0.3  mm.  Re,.  = 3.2  x 106, 

[6 1 12]  roughness 303 

6.62  Spanwise  hot-wire  scan  at  x/c  = 0.10,  Y = 0.45  mm.  Rec  — 3.2  x 10f>, 

[6|12]  roughness 304 

6.63  Spanwise  hot-wire  scan  at  x/c  = 0.15,  Y - 0.65  mm.  Rec  = 3.2  x 106, 

[6|  12]  roughness 305 

6.64  Spanwise  hot-wire  scan  at  x/c  = 0.20,  Y = 0.7  mm.  Re,.  = 3.2  x 106, 

[6 1 12]  roughness 306 

6.65  Spanwise  hot-wire  scan  at  x/c  = 0.25,  Y = 0.7  mm.  Re,  = 3.2  x 1 06 , 

[6 1 12]  roughness 307 

6.66  Spanwise  hot-wire  scan  at  x/c  = 0.29,  Y = 0.65  mm.  Re,.  - 3.2  x 106, 

[6|  12]  roughness 308 

6.67  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.05,  Y = 

0.3  mm.  Re,.  — 3.2  x 106,  [6|12]  roughness 309 

6.68  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.10,  Y = 

0.45  mm.  Rec  = 3.2  x 10f>,  [6 1 1 2]  roughness 310 

6.69  Power  spectral  density  of  spanwise  hot-wire  scan  at.  x/c  — 0.15,  Y = 

0.65  mm.  Re,.  = 3.2  x 106,  [6|12]  roughness 311 

6.70  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.20,  Y = 

0.7  mm.  Rec  = 3.2  x 106,  [6|  1 2]  roughness 312 

6.71  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.25,  Y = 

0.7  mm.  Rec  = 3.2  x 10G,  [6|  12]  roughness 313 

6.72  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  — 0.29,  Y = 

0.65  mm.  Rev  = 3.2  x 10°,  [6 1 12]  roughness 314 
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roughness 318 

Streamwise  velocity  contours  at  x/c  = 0.10.  Rer  = 2.4  x 10(’,  [6|36] 
roughness 319 
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roughness 320 
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roughness 321 
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roughness 322 
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roughness 323 
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roughness 324 
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roughness 325 

Streamwise  velocity  contours  at  x/c  = 0.45.  Re,  = 2.4  x 10,;.  [6|36] 
roughness 326 

Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 
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Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 
span  of  99  mm  at  x/c  = 0.15.  Rec  = 2.4  x 106,  [6|36]  roughness 329 
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6.88  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.20.  Rer  = 2.4  x 106,  [6|36]  roughness 330 

6.89  Spanwise  array  of  100  mean-flow  bound  ary- layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.25.  Rec  — 2.4  x 106,  [6|36]  roughness 331 

6.90  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  — 0.30.  Rec  = 2.4  x 106,  [6|36]  roughness 332 

6.91  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  — 0.35.  Rec  = 2.4  x 106,  [6|36]  roughness 333 

6.92  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.40.  Rec  = 2.4  x 106,  [6| 36]  roughness 334 

6.93  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a 

span  of  99  mm  at  x/c  = 0.45.  Rec  = 2.4  x 106,  [6|36]  roughness 335 

6.94  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  = 0.05.  Rec  = 2.4  x 106,  [6 j 36]  roughness 336 

6.95  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  = 0.10.  Rec  — 2.4  x 106,  [6 1 36]  roughness 337 

6.96  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  - 0.15.  Rec  - 2.4  x 1 06 , [6 ] 36]  roughness 338 

6.97  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at,  x/c  = 0.20.  Rec  = 2.4  x 10(\  [6|36]  roughness 339 

6.98  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  = 0.25.  Rec  = 2.4  x 106,  [6 1 36]  roughness 340 

6.99  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  = 0.30.  Rec  = 2.4  x 106,  [6] 36]  roughness 341 

6.100  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  = 0.35.  Rec  — 2.4  x 106,  [6 1 36]  roughness 342 

6.101  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  — 0.40.  Rec  — 2.4  x 106,  [6|36]  roughness 343 

6.102  Spanwise  array  of  100  disturbance  profiles  covering  a span  of  99  mm 

at  x/c  = 0.45.  Rec  = 2.4  x 106,  [6|36]  roughness 344 


xxiv 


Page 


Figure 

6.103  Disturbance  velocity  contours  at  x/c  = 0.45.  Rec  = 2.4  x 10(i,  [6|36] 

roughness 345 

6.104  Stationary  crossflow  mode  shapes  for  Rec  — 2.4  x 106  and  [6|36]  roughness.  346 

6.105  Total  disturbance  amplitude  and  amplification  factor  N for  Rec  = 

2.4  x 10b  and  [6|36]  roughness 347 

6.106  Spanwise  hot-wire  scan  at  x/c  = 0.10,  Y = 0.4  mm.  Rec  = 2.4  x 106, 

[6 1 36]  roughness 348 

6.107  Spanwise  hot-wire  scan  at  x/c  = 0.15,  Y — 0.8  mm.  Rec  = 2.4  x 106, 

[6|36]  roughness 349 

6.108  Spanwise  hot-wire  scan  at  x/c  = 0.20,  Y = 0.9  mm.  Rec  = 2.4  x 106, 

[6 1 36]  roughness 350 

6.109  Spanwise  hot-wire  scan  at  x/c  = 0.25,  Y = 1.0  mm.  Rec  = 2.4  x 10b, 

[6|36]  roughness 351 

6.110  Spanwise  hot-wire  scan  at  x/c  = 0.30,  Y = 1.0  mm.  Re,,  = 2.4  x 10b, 

[6 1 36]  roughness 352 

6.111  Spanwise  hot-wire  scan  at  x/c  = 0.35,  Y - 1.0  mm.  Rer  = 2.4  x 10b, 

[6 1 36]  roughness 353 

6.112  Spanwise  hot-wire  scan  at  x/c  = 0.40,  Y = 0.95  mm.  Rec  = 2.4  x 10b, 

[6|36]  roughness 354 

6.113  Spanwise  hot-wire  scan  at  x/c  = 0.45,  Y = 0.9  mm.  Re,  = 2.4  x 106, 

[6|36]  roughness 355 

6.114  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.10,  Y = 

0.4  mm.  Rec  = 2.4  x 106,  [6 136]  roughness 356 

6.115  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.15,  Y — 

0.8  mm.  Rec  = 2.4  x 106,  [6|36]  roughness 357 

6.116  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  — 0.20,  Y = 

0.9  mm.  Rec  = 2.4  x 106,  [6|36]  roughness 358 

6.117  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.25,  Y =- 

1.0  mm.  Rec  — 2.4  x 10b,  [6 1 36]  roughness 359 

6.118  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.30,  Y = 

1.0  mm.  Re,  — 2.4  x 10°,  [6|36]  roughness 360 
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6.119  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  — 0.35,  Y = 

1.0  mm.  Rec  = 2.4  x 106,  [6|36]  roughness 361 

6.120  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  — 0.40,  Y = 

0.95  mm.  Rec  = 2.4  x 106,  [6|36]  roughness 362 

6.121  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.45,  Y = 

0.9  mm.  Rec  = 2.4  x 106,  [6|36]  roughness 363 

6.122  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.45,  Y = 

0.9  mm.  Rec  = 2.4  x 106,  [6| 36]  roughness.  Enhanced  view  showing 
amplitude  of  small-wavelength  disturbances 364 

6.123  Individual-mode  disturbance  amplification  factor  N for  Rec  — 2.4  x 

106  and  [6|36]  roughness 365 

6.124  Individual-mode  disturbance  amplitude  for  Rec  — 2.4  x 10fi  and  [6|36] 

roughness 366 

6.125  Streamwise  velocity  contours  at  x/c  — 0.05.  Rec  = 2.4  x 106,  [18 j 12] 
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NOMENCLATURE 


A disturbance  amplitude 

A0  reference  disturbance  amplitude 

C =1.83  m,  airfoil  streamwise  chord  (along  X axis) 

Cp  pressure  coefficient 

CP2  unswept  pressure  coefficient 

Cp 3 swept  pressure  coefficient 

Cf  hot-wire  temperature-compensation  coefficient 

c = 1.29  m,  airfoil  normal  chord  (along  x axis) 

E hot-wire  anemometer  output  voltage 

Ec  equivalent  anemometer  output  voltage  at  the  hot-wire  calibration  temper- 

ature 

/ frequency  [Hz] 

fc  Nyquist  critical  frequency  [Hz] 

H shape  factor 

h[Q  larger  of  two  heights  where  crossflow  velocity  wt  is  10%  of  maximum  [m] 

k roughness  height  [in] 

M Mach  number;  number  of  FFT  averages 

m mode  number 

N = \n(A/A0),  amplification  factor;  number  of  samples 

n curve-fit  order 

p static  pressure  [Pa] 

Poo  freestream  static  pressure  [Pa] 
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Re,. 


= U^Cjv.  streamwise  chord  Reynolds  number 


Rf'cf  = wtmxxhio/u,  crossflow  Reynolds  number  based  on  maximum  crossflow 

velocity  and  second  10%  point  in  crossflow  profile 

Rek  — U(k)k/u,  roughness  Reynolds  number  based  on  total  velocity  at  top  of 

roughness  and  roughness  height 

Reo  — ut6/u,  attachment-line  Reynolds  number  based  on  edge  velocity  and 

momentum  thickness  of  attachment-line  boundary  layer 

S length  (in  time  [s]  or  space  [m])  of  a data  set;  spanwise  extent  of  a hot-wire 

scan  [m] 

T temperature  [°C] 

Tc  hot-wire  calibration  temperature  [°C] 

U = [a2  + w2)1'2  = U(y),  magnitude  of  total  velocity  [m/s] 

Ut.  boundary-layer  edge  velocity  along  A'  axis  [m/s] 

l x;  freestream  velocity  along  A'  axis  [m/s] 

//,  v,  w velocity  components  in  global  (A".  Y.  Z)  coordinates  [m/s] 

u\  v\  w'  disturbance  velocity  components  in  ( X,Y,Z ) coordinates  [m/s] 

■un,  vn,  wn  velocity  components  in  model-oriented  ( x,y,z ) coordinates  [m/s] 

ut,  vt,  ivt  velocity  components  in  boundary-layer  (xt.  ytjzt)  coordinates  [m/s] 

A . Y Z global  test-section  coordinates:  A'  is  along  the  upstream  flow  axis,  Y is 
normal  to  test-section  side  wall,  Z is  unswept  spanwise  coordinate  (positive 
down) 

x,  y, model-oriented  coordinates:  x is  chordwise,  y is  normal  to  chord  line,  z is 
spanwise  (parallel  to  leading  edge) 

xt,yt,zt  boundary-layer  coordinates:  x,  is  along  outer  streamline,  yt  is  normal  to 
model  surface,  zt  is  crossflow  direction 

r>  airfoil  angle  of  attack;  chordwise  (x)  wavenumber 

3 spanwise  ( z ) wavenumber 
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change  or  difference,  e.g.,  AT  is  a temperature  difference 

sampling  interval 

spectral  frequency  resolution 

spectral  wavelength  resolution 

boundary-layer  thickness  [m] 

displacement  thickness  [m] 

flap  deflection  angle 

sweep  angle 

wavelength  [m] 

crossflow  disturbance  wavelength  measured  normal  to  vortex  axis  [m] 

= A, 

wavelength  measured  in  swept  span  (z)  direction  [m] 

kinematic  viscosity  [m2/s] 

momentum  thickness  [m] 

density  [kg/m3] 

frequency  [rad/s] 

body-intrinsic  coordinates:  £ is  chordwise  and  tangent  to  airfoil  surface,  £ 
is  normal  to  airfoil  surface,  z is  spanwise  (parallel  to  leading  edge) 
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CHAPTER  1 


Introduction 


1.1  Background  and  Motivation 

The  understanding,  prediction,  and  eventual  control  of  the  processes  that  cause  a 
boundary  layer  to  transition  from  laminar  to  turbulent  flow  are  in  the  group  of  the 
most  important  unsolved  problems  in  fluid  mechanics.  This  is  due  in  large  part  to  the 
vast  array  of  practical  engineering  applications  that  depend  strongly  on  the  state  of 
the  boundary  layer.  These  include,  but  certainly  are  not  limited  to,  nose  cone  and  heat 
shield  requirements  on  reentry  vehicles,  efficiency  and  performance  of  turbine  cascades 
and  turbomachinery  systems,  convective  heat  transfer  and  temperature  control,  and 
skin-friction  drag  reduction.  The  last  of  these  areas  has  profound  implications  for 
civil  and  commercial  aviation.  Several  estimates  indicate  that  a 25%  reduction  in  fuel 
consumption  would  be  achieved  by  maintaining  laminar  flow  on  the  wings  of  modern 
transport  aircraft  (Pfenninger  1977;  Thomas  1985;  Saric  1994b).  In  light  of  this,  it 
is  no  surprise  that  boundary-layer  stability  and  transition  have  received  considerable 
attention  throughout  this  century.  Yet  in  spite  of  all  the  theoretical  and  experimental 
efforts,  no  mathematical  model  exists  that  can  predict  the  transition  Reynolds  number 
for  a zero-pressure-gradient  flow  over  a smooth  flat  plate  (Saric  1994c).  This  is  not  to 


say,  however,  that  significant  progress  in  Laminar  Flow  Control  (LFC)  has  not  been 
made.  Wall  suction,  heating  (in  water)  or  cooling  (in  gases),  and  careful  shaping 
of  the  pressure  distribution  can  delay  transition  by  limiting  the  growth  of  unstable 
disturbances.  Thus,  while  important  advances  in  transition  research  have  been  made, 
the  ultimate  success  of  transition  prediction  and  LFC  requires  a more  fundamental 
understanding  of  the  processes  that  lead  to  transition. 

1.1.1  Boundary-Layer  Transition 

Although  the  transition  from  laminar  to  turbulent  flow  is  complicated  by  many  fac- 
tors, the  process  for  boundary  layers  in  external  flows  is  usually  divided  into  three 
phases.  The  first  involves  the  mechanisms  by  which  freestream  disturbances  enter 
the  boundary  layer  and  is  called  receptivity  (Morkovin  1969).  This  process  is  still 
not  well  understood,  but  is  arguably  the  most  important  in  boundary-layer  transition 
as  it  provides  the  critical  initial  amplitude,  frequency,  and  phase  for  unstable  waves 
(Saric  1994c).  The  initial  conditions  are  known  to  come  from  external  disturbances 
in  the  form  of  freestream  fluctuations  (both  acoustic  and  vortical),  surface  roughness, 
and  surface  vibration.  The  appearance  of  a particular  instability  mode,  however,  is 
also  influenced  by  several  other  factors  including  Reynolds  number,  wall  curvature, 
wing  sweep,  and  pressure  gradient.  The  combination  of  these  factors  can  cause  a 
variety  of  instabilities  to  occur,  and  the  receptivity  of  the  boundary  layer  to  these 
various  modes  can  dramatically  affect  the  details  of  transition. 

The  second  phase  of  the  transition  process  concerns  the  initial  growth  of  small 
disturbances  in  the  boundary  layer  and  is  described  by  linear  stability  theory.  The 
growth  of  unstable  waves  is  modeled  with  a sot  of  linear,  unsteady,  disturbance  equa- 
tions obtained  from  the  governing  nonlinear  Navier-Stokes  equations.  The  well-known 
Orr-Sommerfeld  equation  for  incompressible,  parallel  basic  states  is  the  best  example, 
however  similar  equations  can  be  derived  for  more  general  flows.  The  linear  regime  is 
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the  most  studied  stage  of  boundary-layer  transition  and  is,  at  least  conceptually,  well 
understood  for  two-dimensional  flows  (Saric  1992b).  Three-dimensional  flows,  how- 
ever, exhibit  fundamentally  different  stability  characteristics  (Reed  and  Saric  1989) 
and  have  proven  a greater  challenge  for  linear  theory. 

The  third  and  final  phase  of  boundary-layer  transition  is  characterized  by  nonlin- 
ear interactions  among  multiple  instability  modes.  This  occurs  when  the  disturbances 
becomes  large  enough  to  interact  with  each  other  through  the  nonlinear  terms  in 
the  Navier-Stokes  equations.  These  interactions  can  distort  the  basic-state  boundary 
layer  leading  to  the  rapid  growth  of  secondary  instabilities  and  the  onset  of  turbulence. 
Because  the  nonlinear  interactions  are  initially  characterized  by  double  exponential 
growth  (Saric  1992b),  it  is  generally  believed  that  transition  occurs  very  quickly  after 
the  development  of  nonlinear  effects.  As  we  will  see  in  chapters  5 and  6,  this  is  not 
necessarily  true. 

1.1.2  Transition  Prediction 

True  transition  prediction  must  account  for  all  three  stages  of  the  transition  process 
described  above.  However,  due  to  the  relative  simplicity  and  computational  efficiency 
of  the  linear  disturbance  equations,  transition  prediction  schemes  are  typically  based 
on  linear  theory.  By  far  the  most  popular  of  these  techniques  is  the  celebrated  eA 
method  of  Smith  and  Gamberoni  (1956)  and  van  Ingen  (1956).  Complete  reviews 
can  be  found  in  Arnal  (1984,  1992,  1994),  Mack  (1984),  and  Saric  (1992b,  1994c). 
The  basic  assumptions  of  the  method  are  (1)  there  is  some  uniform  norm  of  initial 
amplitude,  (2)  there  exists  a critical  disturbance  amplitude  at  transition,  and  (3) 
this  amplitude  is  achieved  through  the  exponential  growth  described  by  linear  theory. 
Within  this  framework,  the  disturbance  growth  rates  computed  from  linear  theory  are 
integrated  from  the  initial  point  of  instability  to  the  transition  location  (as  provided 
by  experimental  data)  to  give  the  amplitude  ratio  A/A0.  The  natural  log  of  this 
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ratio  is  the  amplification  factor  N (commonly  called  the  A-factor).  For  a known  flow 
situation,  the  A-factor  will  have  some  value  (say,  9)  at  transition.  In  unknown  flow 
situations  or  those  for  which  experimental  data  are  not  available,  the  growth  rates 
are  integrated  up  to  the  location  where  A — 9 (or  whatever  value  of  A was  deemed 
appropriate),  at  which  point  transition  is  assumed  to  occur.  In  this  regard,  the  eN 
method  would  be  better  described  as  providing  transition  correlations  rather  than 
transition  predictions. 

The  eN  method  finds  its  greatest  utility  when  used  as  a comparative  measure  of 
the  role  of  stabilizing  or  destabilizing  effects  within  the  boundary  layer;  i.e.,  heat- 
ing/cooling, suction/blowing,  pressure  gradients,  curvature,  etc.  As  a transition- 
prediction  scheme,  however,  the  limitations  of  the  eN  method  are  obvious.  Since 
linear  theory  can  only  calculate  the  amplitude  ratio  between  two  locations,  initial 
conditions  cannot  be  taken  into  account  and  the  receptivity  process  is  entirely  ig- 
nored. Moreover,  nonlinear  interactions  are  not  considered.  Thus,  the  eA  method 
can  be  expected  to  fail  for  flows  in  which  these  effects  are  important.  For  example, 
Radeztskv  et  al.  (1993a)  show  that  small  changes  in  the  surface  roughness  on  a swept 
airfoil  can  dramatically  change  the  transition  location  as  well  as  the  A-factor  at  tran- 
sition. In  general,  e'v  correlations  work  within  some  error  limits  only  for  flows  with 
identical  disturbance  environments,  and  the  use  of  this  method  without  the  support 
of  experimental  data  is  particularly  dangerous  (Saric  1994c). 

For  flows  in  which  the  eN  method  is  known  to  fail  (such  as  crossflow-dominated 
boundary  layers),  the  recently  developed  parabolized  stability  equations  (PSE)  prom- 
ise' to  be  an  effective  tool  for  transition  modeling.  Herbert  (1994)  gives  a comprehen- 
sive' analysis  of  the  PSE.  The  formulation  results  in  an  initial-boundary-value  problem 
that  can  be'  solved  by  numerical  marching.  The  initial  conditions  must  be  prescribed, 
hence  the  PSE  do  not  address  the  receptivity  problem.  However,  nonparallel  effects 
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are  taken  into  account  and  the  nonlinear  terms  can  be  retained.  Thus,  the  PSE  can 
analyze  the  nonlinear  response  of  forced  modes.  This  represents  a major  step  toward 
the  goal  of  modeling  all  three  phases  of  the  transition  process.  At  the  same  time,  the 
lack  of  initial  conditions  underscores  the  need  for  accurate  and  detailed  experiments. 

1.2  Three-Dimensional  Boundary  Layers 

1.2.1  Swept-Wing  Flows 

The  study  of  three-dimensional  boundary  layers  is  motivated  by  the  need  to  under- 
stand the  fundamental  instability  mechanisms  that  cause  transition  in  swept-wing 
flows.  Research  has  identified  four  types  of  instabilities  for  these  flows:  attachment 
line,  streamwise,  centrifugal,  and  crossflow.  The  attachment-line  problem  is  caused 
by  a basic  instability  of  the  attachment-line  boundary  layer  or  by  its  contamination 
with  turbulent  disturbances  and  develops,  in  general,  on  swept  wings  with  a large 
leading-edge  radius  (Poll  1979,  1984,  1985;  Hall  et  al.  1984;  Hall  and  Malik  1986). 
The  streamwise  instability  is  not  unlike  the  familiar  Tollmien-Schlichting  wave  in 
two-dimensional  flows.  This  mechanism  is  associated  with  the  chordwise  velocity 
component  and  is  generally  stabilized  by  a favorable  pressure  gradient.  Centrifugal 
instabilities  can  appear  over  concave  regions  on  the  surface  and  result  in  the  devel- 
opment of  Gortler  vortices  (Floryan  1991;  Benmalek  and  Saric  1994;  Saric  1994a). 
Crossflow  waves,  on  the  other  hand,  are  an  inviscid  instability  mechanism  caused  by 
the  combined  effect  of  wing  sweep  and  pressure  gradient.  All  of  these  instabilities  can 
appear  individually  or  together  depending  on  the  combination  of  Reynolds  number, 
wall  curvature,  wing  sweep,  pressure  gradient,  and  external  disturbances  (including 
surface  roughness).  Thus,  the  swept  wing  provides  a rich  environment  in  which  to 
study  the  stability  behavior  of  three-dimensional  boundary  layers. 
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1.2.2  Crossflow  Instability 


The  present  experiment  focuses  on  the  crossflow  instability  that  occurs  on  swept 
wings  in  regions  of  strong,  favorable  pressure  gradient.  The  physical  mechanism 
for  the  instability  is  described  as  follows.  The  potential-flow  streamlines  are  highly 
curved  near  the  leading  edge  due  to  the  combined  effect  of  wing  sweep  and  pressure 
gradient.  These  streamlines  are  deflected  as  they  pass  over  the  airfoil,  first  inboard 
near  the  leading  edge,  then  outboard  in  the  pressure  recovery  region  downstream 
of  the  pressure  minimum.  Because  of  the  loss  of  streamwise  momentum  near  the 
surface,  the  deflection  is  greater  within  the  boundary  layer  and  the  total  boundary- 
layer  flow  is  not  in  the  direction  of  the  inviscid  streamline.  The  component  of  flow 
perpendicular  to  the  inviscid  streamline  is  called  the  crossflow  velocity.  The  crossflow 
velocity  satisfies  the  no-slip  condition  and  asymptotically  vanishes  at  the  boundaiy- 
layer  edge,  thus  the  profile  contains  an  inflection  point  and  is  subject  to  an  inviscid 
instability.  This  crossflow  instability  manifests  itself  as  co-rotating  vortices  whose 
axes  are  aligned  roughly  with  the  potential-flow  direction. 

Unlike  Tollrnien-Schlichting  instabilities,  the  crossflow  problem  exhibits  stationary 
(/  = ())  as  well  as  traveling  disturbances  that  are  amplified.  Linear  theory  predicts 
much  larger  growth  rates  for  the  traveling  waves,  however  in  many  experiments  tran- 
sition is  dominated  by  the  stationary  disturbances.  Whether  transition  is  controlled 
by  the  stationary  or  traveling  waves  is  intimately  tied  to  the  receptivity  problem. 
Muller  and  Bippes  (1989),  Bippes  (1990,  1991),  and  Bippes  et  al.  (1991)  have  shown 
that  traveling  waves  are  observed  in  environments  rich  in  unsteady  freestream  dis- 
turbances, whereas  stationary  waves  dominate  transition  in  low-disturbance  environ- 
ments. Since  the  low-disturbance  environment  is  more  characteristic  of  flight,  the 
stationary  waves  are  expected  to  be  more  important.  Under  these  conditions,  the 
disturbance  amplitude  can  be  modulated  by  changing  the  surface  roughness  charac- 


teristics  of  the  model  (Kachanov  and  Tararykin  1990;  Radeztsky  et  al.  1993a,  1994). 
In  light  of  this,  one  should  be  very  suspicious  of  simple  transition  prediction  schemes 
(such  as  the  eA  method)  that  do  not  account  for  initial  conditions. 

The  stationary  waves  (that  is,  the  v'  and  w'  disturbances)  are  typically  very  weak, 
hence  many  theoreticians  insist  that  they  can  be  accurately  modeled  with  linear  the- 
ory. However,  experiments  often  show  evidence  of  strong  nonlinear  effects  (Dagenhart 
et  al.  1989,  1990;  Bippes  and  Nitschke-Kowsky  1990;  Bippes  et  al.  1991;  Devhle  et  al. 
1993;  Radeztsky  et  al.  1994).  The  resolution  of  this  apparent  paradox  lies  in  the 
understanding  of  the  physical  mechanism  by  which  the  stationary  waves  disturb  the 
boundary  layer.  The  key  to  the  stationary  disturbance  is  that  the  wave  fronts  are 
fixed  with  respect  to  the  model  and  nearly  aligned  with  the  potential-flow  direction 
(i.e.,  the  wavenumber  vector  is  nearly  perpendicular  to  the  inviscid  streamline).  Con- 
sequently, although  the  (?/,  w')  motion  of  the  wrave  is  weak,  its  stationary  nature 
produces  an  integrated  effect  that  causes  a strong  u'  distortion  in  the  streamwise 
boundary-layer  profile.  In  simple  terms,  the  weak  stationary  wave  “works”  on  the 
same  fluid  to  produce  a large  u'  disturbance  by  convecting  low-speed  fluid  away  from 
the  surface  (where  v'  > 0)  and  high-speed  fluid  toward  the  surface  (where  v'  < 0). 
This  integrated  effect  and  the  resulting  distortion  of  the  mean  boundary  layer  leads 
to  the  modification  of  the  basic  state  and  the  early  development  of  nonlinear  effects. 

An  interesting  side  effect  of  the  stationary  crossflow  waves  is  the  destabilization  of 
secondary  instabilities.  The  v!  distortions  created  by  the  stationary  wave  are  spatial. 
resulting  in  a spanwise  modulation  of  the  mean  streamwise  velocity  profile.  As  the 
distortions  grow,  the  boundary  layer  develops  an  alternating  pattern  of  accelerated, 
decelerated,  and  doubly  inflected  profiles.  The  inflected  profiles  are  inviscidly  un- 
stable and,  as  such,  are  subject  to  a high-frequency  secondary  instability  (Kohama 
et  al.  1991).  This  secondary  instability  is  highly  amplified  and  leads  to  rapid  local 
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breakdown.  Because  transition  develops  locally,  the  transition  front  is  nonuniform  in 
span  and  characterized  by  a saw-tooth  pattern  of  turbulent  wedges. 

1.3  Review  of  Recent  Results 

The  need  to  better  understand  the  transition  processes  for  swept-wing  flows  has 
sparked  strong  interest  in  three-dimensional  boundary  layers  over  the  last  50  years. 
Significant  theoretical  advancements  have  been  made,  however  there  have  been  rel- 
atively few  experiments  detailed  enough  to  validate  the  theory.  The  combination  of 
complex  geometries,  multiple  instability  mechanisms,  and  observed  nonlinear  effects 
has  proven  a formidable  challenge  and  hence  many  issues  are  still  unresolved. 

1.3.1  Literature  Surveys 

There  is  no  shortage  of  publications  in  the  field  of  boundary-layer  stability  and  tran- 
sition; certainly  more  than  can  be  discussed  in  detail  here.  Comprehensive  reviews 
for  both  two-  and  three-dimensional  flows  are  given  by  Arnal  (1984,  1986),  Mack 
(1984),  Poll  (1984),  Saric  (1992b),  and  Reshot ko  (1994).  Reed  et  al.  (1996)  give  an 
up-to-date  discussion  of  effectiveness  and  limitations  of  linear  theory  in  describing 
boundary-layer  instabilities.  The  reader  is  referred  to  these  reports  for  overviews  of 
much  of  the  early  work  in  stability  and  transition.  In  particular,  the  treatise  of  Mack 
(1984)  provides  the  primary  source  of  fundamental  information  on  stability  theory. 

Several  key  papers  provide  in-depth  reviews  of  stability  and  transition  research  in 
three-dimensional  boundary  layers  and,  in  particular,  swept-wing  flows.  Much  of  the 
early  theoretical  and  experimental  work  is  discussed  by  Reed  and  Saric  (1989).  Swept- 
wings.  rotating  disks,  axisymmetric  bodies  (rotating  cones  and  spheres),  corner  flows, 
and  attachment-line  instabilities  are  reviewed,  as  well  as  the  stability  of  flows  for  other 
three-dimensional  geometries.  This  paper  gives  an  excellent  overview  of  the  unique 
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stability  problems  in  three-dimensional  flows.  For  swept  wings,  a historical  account 
of  the  early  investigations  concerning  the  crossflow  instability  is  given,  along  with  a 
detailed  literature  survey.  Poll  (1984)  and  Arnal  (1986)  also  give  extensive  reviews 
of  transition  in  three-dimensional  flows.  Arnal  (1992)  and  Dagenhart  (1992)  focus 
on  swept-wing  flows  and  give  updated  references  for  the  period  between  1989  and 
1992.  Radeztsky  (1994)  gives  a detailed  review  of  the  latest  developments  related 
to  the  crossflow  problem.  Recent  theoretical  and  experimental  efforts  concerning 
nonparallel  and  curvature  effects,  nonlinear  techniques,  secondary  instabilities,  and 
receptivity  issues  are  discussed.  Radeztsky  concentrates  heavily  on  the  work  since 
1990,  and  as  such  his  review  is  still  current. 

These  references  (Poll  1984;  Arnal  1986,  1992;  Reed  and  Saric  1989;  Dagenhart 
1992;  Radeztsky  1994)  combine  to  provide  a thorough  review  of  the  literature  on 
stability  theory  and  transition  with  emphasis  on  three-dimensional  flows  and  the 
crossflow  instability.  Instead  of  repeating  this  material  here,  the  following  discussion 
will  highlight  the  important  developments  directly  related  to  the  specifics  of  this 
experiment. 

1.3.2  Experimental  Investigations 

Although  crossflow  disturbances  have  been  observed  experimentally  since  the  early 
1950s  (Gray  1952),  much  of  the  important  experimental  investigations  have  occurred 
in  the  last  ten  years.  Saric  and  Yeates  (1985)  studied  crossflow  waves  on  a 25° 
swept  flat  plate  using  a “wall  bump”  to  generate  the  necessary  pressure  gradient. 
In  other  experiments,  Poll  (1985)  used  a variable-sweep  cylinder  and  Michel  et  al. 
(1985)  used  a swept  wing  with  a high  aspect  ratio  to  approximate  an  infinite  span 
in  the  measurement  region.  The  Saric  and  Yeates  work  is  significant  in  that  the  first 
superharmonic  of  the  fundamental  stationary  mode  was  observed,  providing  early 
experimental  evidence  of  nonlinear  effects. 
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The  DLR  experiments  of  Bippes  and  co-workers  provide  important  results  con- 
cerning nonlinear  effects  and  initial  conditions.  The  primary  findings  are  reported  by 
Nitschke-Kowsky  and  Bippes  (1988),  Muller  and  Bippes  (1989),  Bippes  (1990,  1991), 
Bippes  and  Muller  (1990),  Bippes  and  Nitschke-Kowsky  (1990),  Muller  (1990),  Muller 
et  al.  (1990),  and  Bippes  et  al.  (1991).  Recent  results  are  summarized  by  Deyhle  et  al. 
(1993),  Lerche  and  Bippes  (1995),  and  Deyhle  and  Bippes  (1996).  These  experiments 
measure  both  stationary  and  traveling  crossflow  waves,  however  their  relative  impor- 
tance in  influencing  the  details  of  transition  is  found  to  depend  on  the  freestream 
turbulence  level.  Muller  and  Bippes  (1989)  describe  a series  of  comparative  exper- 
iments using  the  same  swept  flat  plate  in  both  low-  and  high-turbulence  tunnels. 
The  stationary  waves  are  found  to  dominate  transition  in  the  low-disturbance  en- 
vironment, however  in  the  high-turbulence  tunnel  both  the  growth  rate  and  final 
amplitude  of  the  stationary  disturbance  are  reduced.  At  the  same  time,  the  traveling 
waves  show  larger  growth  rates  and  dominate  transition.  It  is  interesting  that  transi- 
tion is  reported  to  occur  at  slightly  higher  Reynolds  numbers  in  the  higher-turbulence 
environment. 

In  these  experiments,  the  growth  of  the  stationary  and  traveling  waves  shows 
initial  qualitative  agreement  with  linear  theory,  however  the  disturbance  amplitude 
saturates  due  to  nonlinear  effects.  Also,  the  amplitude  of  the  traveling  waves  shows 
a spanwise  modulation  indicating  nonlinear  interactions  with  the  stationary  modes. 
The  observed  wavelength  for  the  stationary  wave  is  in  general  agreement  with  linear 
theory  and  is  independent  of  the  freestream  turbulence  level.  However,  it  is  reported 
that  superposing  a spanwise  periodicity  on  the  flow  fixes  the  wavelength  of  the  station- 
ary disturbance.  Another  important  result  is  that  the  stationary  waves  are  observed 
to  remain  fixed  relative  to  the  model.  This  indicates  that  the  stationary  instability  is 
sensitive  to  initial  conditions  provided  by  surface  roughness,  and  provides  a precursor 
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to  the  ASU  experiments  (discussed  below). 

The  swept-wing  experiment  of  Arnal  et.  al.  (1984)  provides  important  information 
regarding  the  chordwise  evolution  of  the  stationary  wavelength.  In  this  investigation, 
the  crossflow  wavelength  is  observed  to  increase  with  increasing  distance  from  the 
leading  edge,  forcing  individual  vortices  to  “drop-out”  or  vanish  in  order  to  accom- 
modate the  growth  in  wavelength.  Other  crossflow  experiments  at  ONER  A/CERT 
are  reviewed  in  Arnal  and  Juillen  (1987)  and  Arnal  et  al.  (1990).  In  the  latter  work, 
the  difficulties  of  applying  the  e ^ method  to  three-dimensional  flows  are  discussed. 

Using  a swept  flat  plate  and  wall  bump  essentially  identical  to  the  experimental 
configuration  of  Saric  and  Yeates  (1985),  Kachanov  and  Tararykin  (1990)  investigated 
the  effect  of  various  surface  disturbances  on  the  growth  of  stationary  crossflow  waves. 
Spanwise  periodic  (but  constant  in  time)  blowing/suction  as  well  as  isolated  and 
periodic  roughness  were  observed  to  enhance  the  local  distortion  of  the  streamwise 
boundary-layer  velocity.  Increasing  the  height  of  isolated  roughness  increased  the 
local  disturbance  amplitude.  Some  agreement  with  linear  theory  was  achieved  by 
superposing  the  computed  solutions  for  the  fundamental  and  its  first  three  harmonics. 

Continued  investigations  concerning  the  effects  of  initial  conditions  are  reported 
by  Ivanov  and  Kachanov  (1994),  Kachanov  and  Michalke  (1994),  Gaponenko  et  al. 
(1995a,  1995b),  and  Kachanov  (1995).  Gaponenko  et  al.  (1995b)  concentrate  on  the 
receptivity  of  crossflow  disturbances  to  surface  vibrations  produced  with  a metallic 
membrane  oscillated  by  a variable  magnetic  field.  A complex  receptivity  function  is 
defined,  and  it  is  found  that  the  receptivity  “amplitudes”  are  about  twice  as  large  for 
the  most  unstable  crossflow  modes  as  for  the  quasi-two-dimensional  modes.  Kachanov 
(1995)  reports  that  the  disturbance  frequencies,  spanwise  wavenumbers,  and  propa- 
gation angles  are  independent  of  the  properties  of  the  disturbance  generators  (i.e., 
suction/blowing  or  roughness).  Nonlinear  interactions  between  the  stationary  and 
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traveling  crossflow  waves  are  also  examined,  and  a “pumping”  of  energy  either  to 
both  modes  or  to  the  stationary  disturbance  is  observed  prior  to  transition. 

A comprehensive  study  of  stationary  crossflow  waves  is  contained  in  the  experi- 
ments of  Saric  and  co-workers  at  Arizona  State  University.  Saric  et  al.  (1990)  review 
the  design  of  the  experiments,  which  use  a low-aspect-ratio,  45°  swept  wing.  Wall  con- 
tours in  the  test  section  are  used  to  simulate  an  infinite  span,  and  the  low-turbulence 
wind  tunnel  (see  chapter  2)  ensures  that  the  stationary  crossflow  waves  dominate 
transition.  Dagenhart  et  al.  (1989,  1990)  and  Dagenhart  (1992)  report  the  findings 
for  the  original  investigations.  In  these  experiments,  measurements  are  obtained  for 
both  stationary  and  traveling  waves  in  a crossflow-dominated  boundary  layer  under 
the  conditions  of  natural  surface  roughness.  Both  the  growth  rates  and  wavelengths 
for  the  stationary  disturbances  are  found  to  be  smaller  than  predicted  by  linear  the- 
ory. In  contrast  to  Arnal  et  al.  (1984),  no  “drop-outs”  or  other  adjustments  to  the 
stationary  vortex  spacing  are  observed. 

Later  work  by  Radeztsky  et  al.  (1993a)  investigates  the  sensitivity  of  stationary 
crossflow  waves  to  roughness-induced  initial  conditions  by  introducing  micron-sized 
artificial  roughness  elements  near  the  leading  edge.  These  experiments  show  that 
a single  three-dimensional  roughness  element  can  cause  early  local  transition  and 
dramatically  decrease  the  transition  Reynolds  number.  Radeztsky  et  al.  (1994)  con- 
tinued this  work  under  conditions  where  the  natural  roughness  did  not  induce  measur- 
able stationary  crossflow  waves.  In  these  experiments,  spanwise  arrays  of  distributed 
roughness  are  used  to  control  the  disturbance  wavenumber  spectrum.  Even  for  these 
weak  waves,  no  agreement  is  found  with  linear  theory  predictions.  The  experimental 
setup,  however,  was  such  that  very  large  roughness  (k  = 70  150  /mi)  were  required  to 
excite  the  stationary  instability.  Since  the  large  roughness  (Rek  « 100)  probably  pro- 
duced local  nonlinearities,  the  disparity  between  the  experimental  and  linear  growth 
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rates  is  perhaps  no  surprise. 

In  other  experiments  at  ASU,  Kohama  et  al.  (1991)  showed  that  when  the  bound- 
ary layer  is  dominated  by  the  stationary  crossflow  instability,  transition  is  caused  by  a 
high-frequency  secondary  instability.  This  instability  results  from  the  local  distortion 
of  the  mean  streamwise  boundary-layer  profile  by  the  stationary  disturbance.  Thus, 
in  contrast  to  the  conjecture  of  Muller  and  Bippes  (1989)  who  argued  for  the  im- 
portance of  the  traveling  wave,  the  stationary  wave  provides  the  important  physical 
mechanism  that  ultimately  leads  to  transition  in  low-disturbance  environments. 

1.3.3  Theoretical  Developments 

One  common  theme  runs  through  all  of  the  experiments  discussed  above:  the  im- 
portance of  nonlinear  effects  and  receptivity/initial  conditions.  These  characteristics 
of  swept-wing  instabilities  have  motivated  extensive  theoretical  and  computational 
efforts  in  an  attempt  to  provide  improved  transition  prediction  for  three-dimensional 
flows.  Several  recent  advances  have  been  made  that  have  a direct  impact  on  crossflow- 
dominated  boundary  layers.  These  developments  address  nonparallel  effects,  stream- 
line and  body  curvature,  nonlinear  growth,  secondary  instabilities,  and  receptivity. 
The  reader  is  referred  to  Radeztsky  (1994)  for  a general  review  of  recent  theoreti- 
cal and  computational  efforts  in  these  areas.  The  present  discussion  concentrates  on 
nonlinear  techniques  and  receptivity  issues  as  they  apply  to  the  present  experiment. 

An  important  theoretical  result  concerning  stationary  crossflow  waves  is  reported 
by  Reed  (1988),  who  included  the  primary  distortion  of  the  basic  state  and  pre- 
dicted the  spanwise  wavenumber  doubling  observed  in  the  Saric  and  Yeates  (1985) 
experiments.  More  recently,  direct  numerical  simulations  (DNS)  and  the  parabolized 
stability  equations  (PSE)  represent  important  advances  in  stability  and  transition 
modeling.  These  methods  account  for  nonparallel  and  nonlinear  effects,  which  al- 
low for  the  correct  spatial  evolution  of  the  stationary  crossflow  wave  as  well  as  its 
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distortion  of  the  basic  state. 

DNS  have  historically  been  constrained  by  computer  resources  and  algorithmic 
limitations,  however  some  successes  have  been  achieved  in  relation  to  the  station- 
ary crossflow  problem.  Reed  and  Lin  (1987)  and  Lin  (1992)  performed  numerical 
simulations  for  stationary  waves  on  an  infinite-span  swept  wing  similar  to  the  ASU 
experiments  discussed  above.  Meyer  and  Kleiser  (1990)  investigated  the  disturbance 
interactions  between  stationary  and  traveling  crossflow  modes  using  Falkner-Scan- 
Cooke  similarity  profiles  for  the  basic  state.  The  results  were  compared  to  the  experi- 
ments of  Miiller  and  Bippes  (1989).  With  an  appropriate  initial  disturbance  field,  the 
nonlinear  development  of  stationary  and  traveling  crossflow  modes  was  simulated  rea- 
sonably well  up  to  transition.  Wintergerste  and  Kleiser  (1995)  continue  this  work  by 
using  DNS  to  investigate  the  breakdown  of  crossflow  vortices  in  the  highly  nonlinear 
final  stages  of  transition. 

With  the  continued  development  of  new  and  powerful  computers  and  numerical 
methods,  DNS  are  playing  an  increasingly  important  role  in  transition  modeling. 
Kleiser  (1991)  reviews  the  literature  and  Reed  (1994)  discusses  the  details  of  spatial 
DNS.  Reed  covers  nonparallel,  nonlinear,  and  three-dimensional  effects,  as  well  as 
considerations  for  compressibility,  pressure  gradient,  and  surface  geometry.  Important 
discussions  concerning  receptivity  are  also  included. 

The  recently  developed  PSE  appear  poised  to  replace  traditional  linear  theory  as 
the  state-of-the-art  tool  for  boundary-layer  stability  analyses.  Herbert  (1994)  gives 
a detailed  description  of  the  PSE.  The  formulation  results  in  a system  of  parabolic 
differential  equations  describing  the  disturbance  motion.  This  allows  for  the  proper 
spatial  evolution  of  disturbance  modes  as  opposed  to  the  traditional  patching  of  lo- 
cal solutions.  Moreover,  the  nonlinear  terms  can  be  retained  to  provide  full  nonlinear 
stability  analyses.  The  parabolic  nature  of  the  PSE  allows  the  use  of  computationally 
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efficient  marching  algorithms,  hence  PSE  solutions  are  generated  in  a fraction  of  the 
time  required  by  DNS.  On  the  down  side,  the  initial  conditions  must  be  specified,  thus 
the  PSE  do  not  address  the  receptivity  problem.  However,  the  nonlinear  response  of 
forced  modes  can  by  studied  by  varying  the  initial  conditions.  With  proper  guidance 
from  careful  experiments,  the  PSE  have  the  ability  to  accurately  model  nonlinear 
effects  in  three-dimensional  boundary  layers.  For  swept-wing  flows,  nonlinear  PSE 
calculations  exhibit  the  disturbance  amplitude  saturation  characteristic  of  the  DLR 
and  ASU  experiments.  Wang  et  al.  (1994)  investigate  both  stationary  and  traveling 
crossflow  waves  for  the  swept  airfoil  used  in  the  ASU  experiments  and  predict  non- 
linear amplitude  saturation  for  both  types  of  disturbances.  It  is  suggested  that  the 
interaction  between  the  stationary  and  traveling  waves  is  an  important  aspect  of  the 
transition  process.  Other  examples  of  PSE  applied  to  swept  wings  can  be  found  in 
Stuckert  et  al.  (1993),  Schrauf  et  al.  (1995),  and  Haynes  and  Reed  (1996). 

1.4  Experimental  Goals 

As  the  discussions  in  this  chapter  have  illustrated,  stability  and  transition  in  three- 
dimensional  boundary  layers  is  a complicated  process  with  many  unanswered  ques- 
tions. In  crossflow-dominated  flows,  the  importance  of  nonlinear  effects  and  their 
dependence  on  initial  conditions  is  not  well  understood.  These  issues  must  be  ad- 
dressed for  the  future  development  of  transition  prediction  techniques  and,  ultimately. 
LFC.  Thus,  the  goals  of  the  present  investigation  are  (1)  to  study  the  fundamental 
physics  of  the  stationary  crossflow  instability  that  leads  to  transition  on  swept  wings 
in  low-disturbance  environments,  and  (2)  to  provide  a detailed  and  accurate  experi- 
mental database  for  the  development  of  current  analytical  tools  such  as  the  nonlinear 
PSE. 

The  focus  of  the  present  investigation  is  to  study  the  effects  of  distributed  sur- 
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face  roughness  on  the  (nonlinear)  development  of  stationary  crossflow  waves.  Not 
only  will  this  provide  important  input  for  receptivity  studies,  but  the  controlled  ini- 
tial conditions  and  detailed  disturbance  measurements  will  supply  theoreticians  with 
critical  experimental  data  for  code  validation.  The  airfoil  model  and  test  condi- 
tions (discussed  in  chapter  3)  are  chosen  so  that  the  boundary  layer  is  subcritical  to 
attachment-line,  Tollmien-Schlichting,  and  Gortler  instabilities,  while  crossflow  waves 
are  strongly  amplified.  The  extremely  low  turbulence  levels  of  the  ASU  Unsteady 
Wind  Tunnel  ensure  that  the  stationary  waves  dominate  the  transition  process.  In 
order  to  investigate  the  effects  of  controlled  roughness,  the  aluminum  surface  of  the 
model  is  hand  polished  to  a 0.25  //m  rms  finish.  This  allows  the  use  of  micron-sized 
artificial  roughness  elements  to  control  the  wavenumber  spectrum  of  the  stationary 
disturbance  without  saturating  the  initial  disturbance  amplitude. 

The  experimental  methods  are  designed  with  two  objectives  in  mind.  The  first  is 
to  document  the  detailed  structure  of  the  stationary  waves.  This  provides  important 
information  on  the  global  flowfield,  including  total  disturbance  mode  shapes  and 
amplitude  distributions.  The  second  objective  of  the  measurements  is  to  isolate  and 
track  the  growth  of  individual  crossflow  modes.  These  data  are  used  to  investigate 
nonlinear  interactions  among  various  modes,  and  allow  accurate  comparisons  with 
single-wavelength  stability  calculations. 


1.5  Outline 

The  details  of  the  current  experiment  are  presented  in  the  remaining  chapters.  Chap- 
ter 2 describes  the  ASU  Unsteady  Wind  Tunnel  facility,  including  the  instrumentation 
and  measurement  devices  used  for  data  acquisition.  Details  of  the  experimental  de- 
sign and  configuration  are  presented  in  chapter  3.  This  chapter  also  discusses  the 
test  conditions,  for  which  basic-state  and  linear  stability  calculations  are  presented. 
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The  measurement  techniques  used  to  acquire  and  process  the  data  are  described  in 
chapter  4.  Hot-wire  calibration  techniques  and  wind-tunnel  control  methods  are  also 
discussed.  Chapters  5 and  6 present  the  experimental  data  in  detail  and  discuss  the 
results.  Chapter  5 focuses  on  the  baseline  experimental  configuration  and  includes 
comparisons  with  both  linear  and  nonlinear  theoretical  predictions.  Transition  data 
and  basic-state  measurements  are  also  presented.  The  effect  of  Reynolds  number, 
roughness  spacing,  and  roughness  height  on  the  growth  of  the  stationary  waves  are 
investigated  in  chapter  6.  Chapter  7 gives  the  conclusions. 


CHAPTER  2 


Wind-Tunnel  Facility 


2.1  Unsteady  Wind  Tunnel 

The  present  experiment  is  conducted  in  the  Unsteady  Wind  Tunnel  at  Arizona  State 
University.  The  tunnel  was  originally  built  in  1970  by  Dr.  Philip  Klebanoff  and 
calibrated  by  Dr.  James  McMichael  at  the  National  Bureau  of  Standards  in  Gaithers- 
burg, Maryland.  In  1984,  the  wind  tunnel  was  relocated  to  Arizona  State  University 
under  the  direction  of  Dr.  William  Saric.  It  became  operational  in  1987,  after  exten- 
sive modifications  designed  to  improve  the  flow  quality.  In  its  present  configuration 
(figure  2.1),  the  facility  operates  as  a low-speed,  low-turbulence,  closed-circuit,  atmo- 
spheric wind  tunnel  in  which  the  stability  and  transition  of  laminar  boundary  layers 
are  investigated.  Saric  (1992a)  gives  a detailed  description  of  the  facility. 

A 150  hp,  variable-speed,  DC  motor  powers  the  wind-tunnel  fan.  The  1.83  m- 
diameter,  single-stage,  axial  fan  (Buffalo  Forge  model  G300F)  has  nine  adjustable- 
pitch  blades  and  eleven  stators.  The  maximum  rated  fan  speed  is  1350  rpm,  giving  a 
maximum  test-section  velocity  of  35  m/s  for  this  experiment.  The  motor  is  controlled 
by  a Mentor  II  digital  DC  drive.  This  controller  is  equipped  with  a R S/232  serial 
interface,  which  is  converted  to  GPIB  for  digital  communication  with  the  wind-tunnel 
computers  (discussed  below).  The  desired  rotational  speed  is  set  by  writing  into  a 
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16-bit.  register,  giving  a speed  resolution  of  0.02  rpm  for  the  G300F  fan.  Once  this 
register  is  loaded,  the  speed  reference  is  exact  and  the  accuracy  is  determined  by 
the  feedback  system.  With  a 1000-line  optical  encoder  and  a PID  control  loop,  the 
Mentor  II  maintains  the  motor  speed  to  within  0.01%  of  the  set  point. 

Several  key  design  features  result  in  very  low  freestream  turbulence  levels.  The 
tunnel  is  lengthened  by  5 meters  over  the  original  design,  allowing  the  primary  diffuser 
to  be  extended.  The  return  loop  is  contoured  to  provide  a smooth  transition  into  the 
fan  inlet.  All  four  corners  are  fitted  with  turning  vanes  (item  4a’  in  figure  2.1). 
The  turning  vanes  are  circular-arc  airfoils  with  a 50  mm  chord  and  40  mm  spacing. 
Upstream  of  the  contraction  cone  the  flow  passes  through  a 76  rnm-thick  aluminum 
honeycomb  Avail”  (item  ‘b1  in  figure  2.1).  The  honeycomb  consists  of  6.35  mm 
hexagonal  cells.  Immediately  downstream  of  the  honeycomb  are  seven  stainless  steel 
screens  (item  h*'  in  figure  2.1).  The  first  five  screens  have  an  open-area  ratio  of  0.7; 
the  last  two  are  seamless  with  a 0.65  open-area  ratio.  The  screens  are  constructed 
with  0.165  mm  wire  on  a 30/inch  mesh  and  are  separated  by  230  mm.  (Additional 
screens  are  placed  in  the  diffuser  and  test-section  recovery  region  to  prevent  stall.) 
Following  the  screens  is  a 1.64  m settling  chamber  where  viscosity  dissipates  the  small- 
scale  turbulent  fluctuations.  The  steel-reinforced  contraction  cone  (5.33:1  contraction 
ratio)  follows  a 5th-order  polynomial  to  eliminate  curvature  discontinuities  at  the 
contraction  entrance  and  exit. 

Special  precautions  are  taken  to  minimize  the  motor-  and  fan-generated  turbu- 
lence. The  aft  end  of  the  motor  is  fitted  with  a nacelle  to  reduce  wake  turbulence. 
Behind  the  nacelle  are  splitter  plates,  which  reduce  the  large-scale  vortical  motion 
created  by  the  fan.  Screens  are  placed  immediately  after  the  splitter  plates  and  at  the 
downstream  end  of  the  diffuser.  These  measures  are  necessary  to  prevent  a diffuser 
stall.  To  minimize  mechanical  vibrations,  the  fan  housing  is  connected  to  the  wind 
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tunnel  with  flexible  rubber  couplings  and  is  supported  on  a concrete  pad  that  is  iso- 
lated from  the  building  foundation.  Similar  measures  isolate  the  test  section.  Finally, 
the  fan/motor  and  test  section  are  on  opposite  sides  of  a sound-insulated  wall  that 
divides  the  building. 

These  design  features  reduce  the  test-section  turbulence  to  exceptionally  low  lev- 
els. Hot-wire  measurements  in  the  freestream  show  that  u' /U^  < 0.02%  (20  m/s, 
2 Hz  high-pass).  More  information  on  the  calibration  of  the  Unsteady  Wind  Tunnel 
is  given  by  Saric  et  al.  (1988)  and  Mousseux  (1988). 

The  Unsteady  Wind  Tunnel  is  equipped  with  two  complete  and  interchangeable 
test  sections,  each  measuring  1.4  m x 1.4  m x 4.9  m.  As  mentioned  above,  flexible 
couplings  provide  the  only  physical  connection  between  the  test  section  and  the  wind 
tunnel.  By  simply  removing  these  couplings,  the  test  section  can  be  easily  rolled  out 
of  the  tunnel  and  the  second  test  section  can  be  rolled  into  place  with  a different 
experiment.  Thus,  one  experiment  can  be  configured  in  the  work  area  while  another 
is  in  the  tunnel. 

Unsteady  flow's  are  generated  using  a unique  double-duct  design.  Opening  a trap 
door  in  the  plenum  diverts  air  from  the  primary  duct  (i.e.,  contraction  cone  and  test 
section)  into  a secondary  duct  located  above  the  test  section.  These  passages  recom- 
bine in  the  recovery  region  downstream  of  the  test  section.  Immediately  before  this 
point,  each  duct  contains  a set  of  rotating  shutters.  The  shutters  in  the  secondary 
duct  lag  those  in  the  primary  by  90°,  allowing  up  to  100%  velocity  fluctuations  at 
25  Hz  while  maintaining  (relatively)  steady  loading  on  the  fan.  Any  number  of  shut- 
ters in  either  duct  can  be  disengaged  to  reduce  the  unsteady  amplitude.  To  function 
as  a steady  wind  tunnel,  the  primary  shutters  are  locked  open  and  the  secondary  duct 
is  closed.  The  present  experiment  operates  in  this  steady  configuration. 


2.2  Computer  Systems 


Every  effort  is  made  to  automate  the  experimental  procedures  at  the  Unsteady  Wind 
Tunnel.  This  improves  operational  efficiency  and  increases  consistency  and  reliability 
by  minimizing  the  ever-present  subjectivity  of  the  researcher.  Central  to  this  task  are 
several  computer  systems,  which  oversee  all  wind-tunnel  operations. 

Table  2.1  lists  the  features  and  capabilities  of  the  Unsteady  Wind  Tunnel  computer 
systems.  All  wind-tunnel  operations  are  controlled  by  the  Sun  SPARCstation  20. 
This  multiprocessor  workstation  is  equipped  with  a GPIB  interface  board,  National 
Instruments  model  GPIB-SPRC-B.  This  versatile  device-communication  bus  allows 
the  Sun  to  interface  with  all  computer-controlled  instrumentation  while  retaining  the 
flexibility  and  expansibility  to  grow  with  changing  needs.  With  this  arrangement,  all 
instrumentation  devices  are  “external”  in  that  they  do  not  reside  within  the  computer 
but  communicate  with  it  digitally  via  the  GPIB.  Consequently,  signal  degradation 
is  minimized  since  the  data-acquisition  equipment  is  placed  close  to  the  experiment, 
thereby  eliminating  the  need  for  long  analog  cables  between  the  wind-tunnel  room 
and  the  control  room.  Moreover,  the  modularity  of  the  entire  system  allows  virtually 
any  component,  including  the  computer,  to  be  replaced  or  upgraded  independently. 

The  SPARCstation  runs  version  2.4  of  the  Solaris  operating  system:  Sun’s  UNIX 
based  on  AT&T’s  System  V,  Release  4.  This  multiuser,  multitasking  environment 
allows  simultaneous  data  acquisition,  analysis,  wind-tunnel  control,  and  program  de- 
velopment. OpenWindows  provides  a modern,  X Windows-based,  graphical  user 
interface  displayed  on  a 20  inch,  25G-color  console  terminal. 

All  data-acquisition,  analysis,  and  wind-tunnel-control  programs  are  written  “in 
house”  using  C,  C++,  or  Lab  VIEW.  A set  of  custom  object-code  libraries  simpli- 
fies the  programming  task  by  providing  a standardized  interface  to  the  acquisition 
and  control  instruments  (Reibert  1996).  Tecplot  is  available  for  plotting,  and  UT^X 
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Table  2.1:  Unsteady  Wind  Tunnel  computer  systems. 


Feature 

Sun 

PC 

Digital 

Macintosh 

System  and  Operating  Environment 

System 

SPARC  20 

ACT486 

DEC  5000 

Quadra 

Model 

612MP 

50 

200 

650 

OS 

Solaris 

SCO  UNIX 

Ultrix 

Mac  OS 

OS  Version 

2.4 

2.0 

4.41 

7.5.3 

GUI 

Open  Windows 

OSF/Motif 

OSF/Motif 

N/A 

Processor 

CPU 

SuperSPARC+ 

80486DX 

R3000 

MC68040 

No.  of  CPU 

2 

1 

1 

1 

Speed  [MHz] 

60 

50 

25 

33 

MIPS 

167a 

12b 

24 

8b 

MFLOPS 

36. 6a 

2b 

3C 

N/A 

Graphics 

System 

TurboGX 

Trident  8900 

8-plane 

Built-in 

Monitor  [inch] 

20 

17 

19 

17 

Resolution 

1024  x 768 

1024  x 768 

1024  x 864 

832  x 624 

Colors 

256 

256 

256  (gray) 

65536 

Memory  and  Storage 

RAM  [MB] 

64 

20 

16 

8 

Disk  [MB] 

3150 

425 

2430 

230 

Tape  [MB] 

5000 

2000 

100 

8000 

Floppy 

• 

• 

• 

CD-ROM 

• 

• 

I/O  Buses  and  Connectivity 

SCSI 

• 

• 

• 

• 

GPIB 

• 

• 

Ethernet 

• 

• 

• 

• 

aValue  for  each  SuperSPARC+  processor  in  a single-CPU  system. 
b Estimated. 
c Double  precision. 


(along  with  the  standard  complement  of  supporting  programs)  is  installed  for  tech- 
nical document  production.  The  SPARCstation  also  functions  as  a server  on  the 
World  Wide  Web.  The  ASU  Wind  Tunnel  Complex  home  page  can  be  reached  at 
http:  / / wtsun.eas.asu.edu . 

A PC  compatible  serves  as  a back-up  data-acquisition  system.  This  computer 
runs  the  Santa  Cruz  Operation’s  Open  Desktop  Server  System  version  2.0.  This  mul- 
tiuser UNIX  environment  is  based  on  AT&T’s  System  V,  version  3.2.4.  The  system 
features  OSF /Motif  X Windows  graphics  displayed  on  a 17  inch,  256-color  console  ter- 
minal. This  machine  is  also  equipped  with  a GPIB  controller,  National  Instruments 
model  AT-GPIB/TNT.  This  allows  the  PC  to  take  over  all  data-acquisition  and  wind- 
tunnel-control  functions  simply  by  moving  the  GPIB  cable  from  tin*  SPARCstation 
to  the  PC.  To  facilitate  this,  the  custom  programming  libraries  used  to  access  the 
wind-tunnel  instrumentation  devices  are  maintained  on  the  PC.  Thus,  the  researcher 
needs  only  to  move  the  GPIB  cable  from  the  SPARCstation  to  the  PC  and  recompile 
his  programs  to  bring  the  back-up  data-acquisition  computer  on  line. 

A DECstation  5000  provides  additional  workstation  capabilities  and  data  storage 
space.  This  unit  runs  Ultrix  version  4.41:  Digital’s  UNIX  based  on  BSD.  The 
OSF/Motif  X Windows  graphics  system  is  displayed  on  a 19  inch,  256-shade,  gray- 
scale console  terminal. 

A Macintosh  Quadra  650  is  available  for  general-purpose  computing.  This  com- 
puter runs  System  7.5.3  of  the  Macintosh  OS  displayed  on  a 17  inch,  64k-color  ter- 
minal. A wide  array  of  software  is  installed  including  drafting  and  word-processing 
applications. 

High-quality  printed  output  is  provided  by  a Hewlett-Packard  HP4M  600  dpi 
PostScript  laser  printer.  A Hewlett-Packard  1200C/PS  300  dpi  PostScript  inkjet 
printer  is  available  for  color  output. 
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All  wind-tunnel  computers  and  printers  are  connected  via  the  Ethernet.  A subnet 
router  connects  the  Unsteady  Wind  Tunnel  to  the  campus  backbone,  which  is  in 
turn  connected  to  the  world-wide  Internet.  The  router  functions  as  a gateway  bv 
differentiating  between  local  transmissions  and  those  intended  for  machines  outside 
the  Unsteady  Wind  Tunnel.  This  creates  a fast  local  network  by  isolating  local  traffic 
from  the  global  network. 

2.3  Instrumentation 

Freestream  flow  conditions  are  determined  by  pressure  and  temperature  measure- 
ments near  the  test-section  entrance  plane.  A 10  torr  differential  pressure  transducer 
(MKS  model  398HD)  measures  dynamic  pressure  from  a Pitot-static  tube.  The  static 
side  of  the  probe  is  also  connected  to  a MKS  390HA  1000  torr  absolute  pressure 
transducer.  Both  temperature-compensated  transducers  are  connected  to  MKS  270B 
14-bit.  signal  conditioners.  These  provide  visual  displays  as  well  as  digital  and  analog 
output  signals,  the  latter  of  which  are  interfaced  with  the  data-acquisition  system. 
The  test-section  temperature  is  measured  with  a thin-film  RTD.  The  calibrated  ana- 
log output  is  also  wired  into  the  data-acquisition  system. 

Hot-wire  anemometry  provides  accurate  boundary-layer  velocity  measurements. 
The  system  consists  of  Dantec  55P15  boundary-layer  probes  and  two  Dantec  55M01 
constant-temperature  anemometers  equipped  with  55M10  CTA  standard  bridges. 
The  hot-wire  probes  use  5 /mi  platinum-plated  tungsten  wires.  The  probe  tines  are 
1.25  mm  apart  and  are  offset  3 mm  from  the  probe  axis  to  facilitate  measurements 
close  to  the  model  surface. 

A two-channel  filter/amplifier  (Stewart  model  VBF44)  provides  analog  signal  con- 
ditioning for  AC  measurements.  This  unit  has  two  high-pass  and  two  low-pass  filter 
responses.  The  low-pass  filters  have  an  AC-couple  option.  Cutoff  frequencies  range 
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from  1 Hz  to  255  kHz,  and  pre-  and  post-filter  gains  provide  a maximum  amplification 
of  70  (IB.  All  features  are  remotely  programmable  through  a RS/232  serial  interface. 
A three-channel  Tektronix  AM502  filter/amplifier  provides  additional  analog  filtering, 
and  a GPIB-eontrolled  Stanford  Research  Systems  SR.530  lock-in  amplifier  measures 
amplitude  and  phase  data. 

All  analog  signals  are  digitized  with  two  IOtech  ADC488/8SA  analog-to-digital 
(A/D)  converters.  Each  A/D  converter  can  simultaneously  sample  and  hold  up  to 
eight  differential  signals  with  16-bit  resolution.  Connected  in  a master/slave  arrange- 
ment, the  two  units  use  a common  clock  trigger  to  provide  a total  of  sixteen  channels 
of  simultaneous  A/D  conversion.  The  input  voltage  range  for  each  channel  is  indepen- 
dently programmable  between  ±1  and  ±10  volts,  thus  “small”  signals  can  be  resolved 
to  33  //\\  The  aggregate  sampling  rate  varies  discretely  from  0.02  Hz  to  100  kHz. 
These  external  A/D  converters  communicate  with  the  data-acquisition  computer  via 
the  GPIB. 

An  IOtech  DAC488HR/4  digital-to-analog  converter  is  available  for  source  signal 
generation.  This  four-channel,  16-bit.  unit  provides  synchronous  analog  output  at 
a maximum  update  rate  of  100  kHz  per  channel.  The  unit  can  also  operate  as  a 
waveform  generator  and  a precision  DC  voltage  source.  All  options  are  remotely 
accessible  via  the  GPIB. 

Other  typical  laboratory  electronic  equipment  is  available,  including  an  eight- 
channel  Tektronix  5440  oscilloscope,  Fluke  8050A  digital  multimeters,  and  various 
signal  generators  and  power  supplies. 

2.4  Traverse  and  Sting 

A computer-controlled,  three-dimensional  instrumentation  traverse  allows  accurate 
positioning  of  the  hot-wire  probes  within  the  test  section.  A detailed  description  of 
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the  system  is  given  by  Radeztsky  (1994).  The  important  features  are  highlighted 
below. 

The  traverse  carriage  (figures  2.2  and  2.3)  is  located  outside  the  test  section  and 
moves  in  the  X (streamwise)  direction  on  two  stainless  steel  Thompson  rails.  Only  the 
instrumentation  sting  penetrates  the  flowfield.  This  is  a critical  design  feature  nec- 
essary to  minimize  the  intrusiveness  of  the  measurement  system.  A moving  traverse 
within  the  test  section  will  cause  local  pressure  field  variations  and  global  flow  adjust- 
ments that  can  dramatically  alter  the  results  of  boundary-layer  stability  experiments 
(Saric  1990). 

The  hot-wire  sting  reaches  into  the  measurement  region  through  a slotted  plex- 
iglass window  contained  within  the  test-section  wall.  A zipper  automatically  opens 
and  closes  the  slot  around  the  sting  when  the  carriage  moves  in  the  X direction.  An 
air-tight  plexiglass  outer  wall  (item  ‘f’  in  figure  2.2)  encloses  the  entire  system  in  a 
pressure  box.  This  equalizes  the  pressure  across  the  interior  test-section  wall,  virtu- 
ally eliminating  any  transverse  loading  on  the  slotted  window.  Mass  transfer  through 
the  small  gaps  in  the  slot  surrounding  the  sting  is  also  minimized. 

The  Y (wall-normal)  and  Z (vertical)  motion  subsystems  are  entirely  contained 
on  the  traverse  carriage.  The  sting  mounts  to  a small  aluminum  sub-carriage  (item  ‘d’ 
in  figures  2.2  and  2.3),  which  is  supported  on  two  parallel  rails.  A high-resolution 
lead  screw  moves  the  Y carriage  normal  to  the  test-section  wall.  Vertical  motion  is 
provided  by  moving  the  Y carriage  with  twin  lead  screws  and  rails,  shown  as  items  ‘b’ 
and  ‘c’  in  figure  2.3.  Another  set  of  twin  lead  screws  and  rails  moves  the  slotted 
window  in  conjunction  with  the  vertical  motion  of  the  Y carriage.  These  movements 
are  coordinated  to  ensure  the  sting  is  always  centered  in  the  slotted  window. 

All  lead  screws  are  driven  by  high-resolution  Compumotor  microstepping  motors. 
1000-line  Renco  optical  encoders  provide  digital  position  feedback  for  all  axes  includ- 
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Table  2.2:  Traverse  system  capabilities. 


Direction 

X 

Y 

Z 

Total  Travel 
Minimum  Step 

1.25  m 
12  /mi 

100  mm 
0.7  fin l 

175  mm 
1.3  fin i 

ing  the  slotted  window.  Quadrature  increases  the  effective  encoder  resolution  by  a 
factor  of  four,  giving  the  minimum  step  sizes  shown  in  table  2.2.  A four-axis  digital 
motion  controller  (Compumotor  model  CM4000)  directs  all  traverse  movements.  This 
microprocessor-based  controller  internally  governs  all  aspects  of  the  motion  control 
including  the  encoder  feedback  loop. 

New  software  algorithms  have  been  implemented  that  greatly  improve  the  accu- 
racy and  reliability  of  the  traverse  system.  The  CM4000  firmware  contains  a BASIC- 
like  programming  language.  Although  awkward  to  use  for  complex  motion  control,  it 
is  sufficient  to  communicate  with  the  data-acquisition  computer  via  the  GPIB.  Dur- 
ing an  experiment,  a resident  program  on  the  CM4000  waits  for  commands  indicating 
a traverse  move  is  desired.  When  these  commands  are  received,  the  controller  initi- 
ates the  motion,  monitors  the  encoder  feedback  and  applies  any  necessary  “correction 
moves” , then  indicates  to  the  computer  when  all  axes  are  within  the  dead-band  tol- 
erance. This  process  repeats  until  the  computer,  at  the  end  of  the  experiment,  tells 
the  CM4000  to  shutdown  the  traverse  system.  The  essential  feature  of  this  method 
is  the  resident  program  on  the  CM4000.  Without  it,  the  controller  cannot  monitor 
the  encoder  feedback  and  position  accuracy  is  potentially  compromised.  This  was  the 
case  with  all  previous  experiments. 

The  hot-wire  sting  (figure  2.4)  is  the  same  used  in  earlier  swept-wing  experiments 
at  the  Unsteady  Wind  Tunnel  (Dagenhart  et  al.  1989,  1990;  Kohama  et  al.  1991; 
Radeztsky  et  al.  1993a,  1994).  The  streamlined,  carbon-composite  body  attaches  to 
the  traverse  via  an  aluminum  mounting  strut.  Two  Dantec  probe  tubes  are  mounted 
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oil  the  composite  section.  The  probe  mount  at  the  tip  of  the  sting  allows  rotation 
about  the  Z axis.  This  motion,  coupled  with  probe  rotation  about  the  tube  axis,  is 
necessary  to  accurately  position  the  hot-wire  relative  to  the  three-dimensional  surface 
of  the  swept  wing.  Radeztsky  (1994)  gives  further  details  concerning  the  sting. 


CHAPTER  3 


Test  Model  and  Experimental  Design 


This  chapter  reviews  the  experimental  design  and  setup.  The  swept-wing  test  model 
and  its  configuration  within  the  test  section  are  described.  The  test  conditions  are 
chosen,  and  stability  calculations  are  presented  for  those  conditions. 

3.1  Background 

As  discussed  in  section  1.3,  the  early  ASU  experiments  (Dagenhart  et  al.  1989,  1990; 
Saric  et  al.  1990)  investigated  the  stability  and  transition  of  swept-wing  boundary 
layers  dominated  by  the  crossflow  instability  mechanism.  In  these  experiments,  the 
initial  conditions  for  the  disturbance  amplitude  came  from  the  unknown  natural 
roughness  of  the  surface.  Later  experiments  by  Radeztsky  et  al.  (1993a)  studied 
the  sensitivity  to  isolated  leading-edge  roughness.  This  work,  however,  focused  on 
transition  location  and  in  large  part  ignored  the  details  of  the  disturbance  growth. 
It  was  not  until  the  experiments  of  Radeztsky  et  al.  (1994)  that  a detailed  and  sys- 
tematic investigation  of  the  effects  of  distributed  surface  roughness  on  (very  weak) 
stationary  crossflow  waves  was  documented. 

The  present  investigation  returns  to  the  configuration  of  the  original  experiments 
in  order  to  study  the  effects  of  distributed  surface  roughness  in  a crossflow-dominated 
boundary  layer.  The  following  sections  discuss  the  details  of  the  experimental  design. 


3.2  Model  Configuration 


3.2.1  Airfoil 

Th<’  NLF(2)-0415  airfoil  (Somers  and  Horstmann  1985)  is  the  same  used  in  all  previous 
swept-wing  work  at  the  Unsteady  Wind  Tunnel  (Dagenhart  et  al.  1989,  1990;  Saric 
et  al.  1990;  Kohama  et  al.  1991;  Radeztsky  et  al.  1993a,  1994).  The  NLF(2)-0415 
is  designed  as  an  unswept,  natural-laminar-flow  airfoil  for  use  on  general  aviation 
aircraft.  The  airfoil  cross  section  and  pressure  distribution  for  the  design  angle  of 
attack  of  0°  are  shown  in  figure  3.1.  The  favorable  pressure  gradient  back  to  the 
pressure  minimum  at  x/c  = 0.71  is  designed  to  maintain  laminar  flow  on  the  upper 
surface'  by  controlling  the  Tollmien-Schlichting  (T-S)  instability. 

The  airfoil  is  swept  45°  for  the  ASU  experiments,  creating  a test  model  well- 
suited  for  the  study  of  three-dimensional  boundary  layers.  The  small  leading-edge 
radius  eliminates  attachment-line  instabilities  for  the  Reynolds  number  range  of  the 
Unsteady  Wind  Tunnel  [Rcq  — 44  at  a = —4°  and  Rec  = 2.4  x 106),  and  the 
absence  of  concave  regions  on  the  upper  surface  suppresses  the  Gortler  instability. 
The  result  is  a nearly  ideal  platform  for  the  investigation  of  crossflow  and/or  T-S 
instabilities.  At  small  negative  angles  of  attack,  the  favorable  pressure  gradient  from 
the  attachment  line  to  the  pressure  minimum  at  x/c  = 0.71  produces  strong  crossflow 
and  stabilizes  the  T-S  modes.  At  o = 0°,  the  pressure  gradient  is  weakly  favorable 
back  to  x/c  = 0.71.  Under  these  conditions,  both  crossflow  and  T-S  disturbances  are 
weakly  amplified.  At  small  positive  angles  of  attack,  the  pressure  minimum  moves 
forward  to  x/c  = 0.02.  and  the  adverse  pressure  gradient  leads  to  strong  growth  of  T-S 
waves.  A 20%-chord.  trading-edge  flap  with  a maximum  deflection  of  ±20°  allows 
further  contouring  of  the  pressure  distribution.  Dagenhart  (1992)  gives  a detailed 
review  of  the  operating  range  for  the  NLF(2)-0415. 
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3.2.2  Test-Section  and  Wall  Liners 

The  swept  NLF(2)-0415  airfoil  is  mounted  vertically  in  a dedicated  test  section  mea- 
suring 1.4  m x 1.4  m x 4.9  m.  The  vertical  orientation  simplifies  instrumentation 
access  since  the  test  surface  (i.e.,  airfoil  upper  surface)  faces  the  front  wall  of  the  test 
section.  The  model  attaches  to  the  test  section  via  a shaft  and  thrust  bearing.  The 
axis  of  the  shaft  is  parallel  to  the  leading  edge  and  is  located  at  xjc  = 0.25.  This  com- 
bination allows  angle-of-attack  rotations  about  the  1/4-chord  line  from  —4°  to  4-4°  in 
1°  increments.  Since  it  is  good  experimental  practice  to  avoid  symmetry  planes,  the 
thrust  bearing  is  positioned  76  mm  off  center  placing  the  pivot  point  610  mm  from 
the  test-section  front  wall  and  760  mm  from  the  rear  wall.  In  addition,  moving  the 
model  closer  to  the  front  wall  reduces  the  span  of  the  instrumentation  sting,  which 
helps  minimize  probe  vibration. 

In  the  interest  of  detailed  measurements,  the  swept  airfoil  chord  of  1.83  m is  chosen 
to  allow  significant  boundary-layer  growth  (<5  ^ 4 mm  in  the  mid-chord  region  for 
moderate  chord  Reynolds  numbers).  Of  course,  wall-interference  effects  cannot  be 
ignored  when  a model  this  size  is  placed  in  a 1.4  m-square  test  section.  One  way 
to  handle  these  effects  is  to  include  the  test-section  walls  in  all  theoretical  models. 
This  is  relatively  straightforward  for  the  front  and  rear  walls  (i.e.,  those  opposite  the 
airfoil  upper  and  lower  surfaces).  However,  including  the  effects  of  the  test-section 
floor  and  ceiling  significantly  complicates  the  computational  effort.  The  simplifying 
assumption  of  spamvise  invariance  cannot  be  used,  and  a fully  three-dimensional  code 
is  required  for  both  the  basic  state  and  the  stability  calculations.  On  the  other  hand, 
the  flowfield  is  spamvise  invariant  (i.e.,  the  boundary-layer  and  stability  characteristics 
are  invariant  along  lines  of  constant  chord)  if  the  airfoil  is  infinite  in  span.  Under 
these  conditions,  the  problem  can  be  modeled  in  two  dimensions  with  the  addition  of  a 
constant  spanwise  velocity  [/ooSin(A)  in  the  z direction.  This  results  in  dramatically 
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more  efficient  computational  methods.  Just  as  important,  the  infinite-span  swept 
wing  produces  a benchmark  configuration  for  crossflow  studies  (not  unlike  the  flat 
plate  for  Tollmien-Schlichting  instabilities). 

In  light  of  this,  the  challenge  now  falls  to  the  experimentalist  to  create  an  infinite- 
span  airfoil  in  a finite  test  section.  At  the  Unsteady  Wind  Tunnel,  this  is  accomplished 
by  contouring  the  test-section  floor  and  ceiling  with  end  liners.  These  liners  create  a 
stream  surface  that  follows  the  inviscid  streamlines  for  an  infinite-span  airfoil.  This 
experiment  uses  the  same  end  liners  constructed  for  the  previous  experiments  at 
n = —4°.  The  reader  is  referred  to  Dagenhart  (1992)  for  a discussion  the  liner 
design  and  Radeztsky  (1994)  for  a detailed  description  of  the  construction  technique. 
The  NLF(2)-0415  airfoil  and  end  liners  for  the  present  configuration  are  shown  in 
figure  3.2. 

The  floor  and  ceiling  liners  are  130  mm  thick  at  the  test-section  entrance.  This 
necessitates  the  installation  of  contraction-cone  fairings,  which  reduce  the  contrac- 
tion exit  area  to  match  the  smaller  test-section  entrance.  The  fairings  begin  at  the 
inflection  point  in  the  original  contraction  contour  and  follow  a 5t.h-order  polynomial 
to  avoid  curvature  discontinuities.  The  contraction  ratio  is  6.55:1  with  the  fairings 
installed.  Symmetric  entrance  flow  is  maintained  since  the  floor  and  ceiling  liners  are 
the  same  thickness  at  the  test-section  entrance. 

3.3  Test  Conditions 

Many  factors  influence  the  choice  of  the  experimental  test  conditions,  but  by  far  the 
most  important  are  the  stability  characteristics  of  the  boundary  layer.  Consequently, 
stability  calculations  are  an  integral  part  of  the  design  process  for  this  type  of  exper- 
iment. Traditionally,  an  exhaustive  study  is  performed  in  which  the  boundary-layer 
behavior  is  analyzed  over  the  entire  parameter  space  of  the  proposed  experiment.  The 
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test  conditions  are  then  chosen  to  provide  the  desired  stability  features.  After  the 
experiment  is  conducted,  the  results  are  compared  to  the  predictions  and  conclusions 
are  drawn  as  to  the  applicability  of  the  theoretical  model. 

The  original  experiments  of  Dagenhart  et  al.  (1989,  1990),  as  well  as  those  of 
Radeztsky  et  al.  (1993a,  1994),  proceeded  in  this  fashion.  The  basic  state  was  com- 
puted with  the  Kaups  and  Cebeci  (1977)  boundary-layer  code,  with  edge  conditions 
supplied  by  the  MCARF  code  (Stevens  et  al.  1971).  Linear,  parallel  stability  predic- 
tions were  then  obtained  using  the  SALLY  (Srokowski  and  Orszag  1977)  and  MARIA 
(Dagenhart  1981)  codes,  from  which  the  test  conditions  were  chosen.  However,  the 
disturbance  measurements  from  these  experiments  bear  little  resemblance  to  the  the- 
oretical predictions.  While  this  is  certainly  instructive  in  its  own  right,  the  failure 
of  linear  theory  under  these  conditions  does  little  to  assist  the  experimentalist  in 
designing  future  tests. 

Fortunately,  the  present  work  has  the  benefit  of  hindsight,  and  can  rely  on  the  pre- 
vious experiments  for  assistance  in  choosing  appropriate  test  conditions.  When  used 
in  conjunction  with  theoretical  predictions,  this  dramatically  improves  the  ability  of 
the  experimentalist  to  pick  operating  conditions  that  exhibit  the  desired  boundary- 
layer  stability  characteristics.  For  the  present  experiment,  this  means  that  the  com- 
bination of  angle  of  attack,  Reynolds  number,  and  surface  roughness  distribution  can 
be  chosen  a priori  with  little  or  no  “guesswork” . 

3.3.1  Angle  of  Attack 

Since  the  present  investigation  focuses  on  the  crossflow  instability,  the  angle  of  attack 
is  set  to  —4°  and  the  flap  is  not  deflected  (<5 p = 0°).  Figure  3.3  shows  the  unswept 
airfoil  contour  and  upper-surface  Cp  distribution  for  this  configuration.  With  the 
continuously  negative  pressure  gradient  from  the  attachment  line  to  x/c  = 0.71, 
the  boundary  layer  is  subcritical  to  T-S  waves  and  transition  is  dominated  by  the 
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stationary  crossflow  instability.  This  is  predicted  by  Dagenhart  (1992)  and  confirmed 
experimentally  by  Dagenhart  et  al.  (1989,  1990)  and  Koharna  et  al.  (1991).  These 
early  experiments  show  that  the  45°-swept  NLF(2)-0415  at  a = -4°  is  a near  perfect 
crossflow  generator. 

3.3.2  Reynolds  Number 

Although  several  factors  influence  the  choice  of  Reynolds  number,  the  final  decision 
represents  a tradeoff  between  crossflow  disturbance  growth,  transition  location,  and 
wind-tunnel  heating.  With  the  present  configuration,  the  maximum  chord  Reynolds 
number  is  3.6  x 106,  however  this  cannot  be  sustained  over  long  runs  due  to  wind- 
tunnel  heating  (see  section  4.1.2).  Fortunately,  at  a = —4°  the  crossflow'  disturbance 
growth  is  strong  even  at  moderate  chord  Reynolds  numbers  so  it  is  not  necessary 
to  maximize  the  tunnel  speed.  The  baseline  Reynolds  number  for  this  experiment 
is  Re.,.  = 2.4  x 106.  This  is  large  enough  to  generate  significant  crossflow,  yet  small 
enough  to  minimize  wind-tunnel  heating  effects.  Moreover,  the  disturbance  growth 
can  be  analyzed  in  detail  since  laminar  flow  is  maintained  beyond  50%  chord.  When 
the  effect  of  Reynolds  number  on  the  disturbance  amplitude  is  desired,  measurements 
are  also  taken  at  Rec  = 1.6  x 10fl  and  3.2  x 1 0h . 

3.3.3  Roughness  Elements 

Whereas  traveling  crossflow'  disturbances  are  influenced  by  freestream  turbulence 
(Midler  and  Bippes  1989;  Muller  1990;  Bippes  et  al.  1991;  Deyhle  et  al.  1993;  Deyhle 
and  Bippes  1996),  the  development  of  stationary  waves  depends  strongly  on  surface 
roughness  near  the  attachment  line  (Radeztskv  et  al.  1993a,  1994).  As  outlined  in 
chapter  1.  the  primary  goal  of  this  experiment  is  to  study  the  growth  of  station- 
ary crossflow  disturbances  and  their  dependence  on  initial  conditions.  It  is  therefore 
essential  to  carefully  control  and  document  the  surface  roughness  distribution. 
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The  first  step  in  this  process  is  to  reduce  the  natural  roughness  of  the  model,  thus 
the  aluminum  surface  of  the  NLF(2)-0415  airfoil  is  highly  polished.  Figure  3.4  shows 
a profilometer  measurement  of  the  surface  finish.  The  0.12  /mi  mis  finish  is  typical  of 
the  roughness  level  near  the  mid-chord  region.  Near  the  leading  edge  the  roughness 
level  is  0.25  /mi  rms.  This  very  smooth  surface  provides  an  ideal  environment  within 
which  the  effect  of  roughness  on  stationary  crossflow  waves  is  investigated. 

Following  Radeztsky  et  al.  (1993a),  the  initial  conditions  are  controlled  by  apply- 
ing roughness  elements  to  the  airfoil  surface  near  the  attachment  line.  The  ability 
of  artificial  roughness  to  induce  stationary  crossflow  waves  is  well-documented  by 
Radeztsky  et  al.  (1993a,  1994).  Two  different  roughness  elements  are  used.  The  first 
are  Geotvpe  #GS-104  circular  “rub-down”  dots  common  in  the  graphic  arts  industry. 
These  3.7  nun-diameter,  dry-transfer  dots  produce  a 6 /un-thick  roughness  element 
when  rubbed  onto  the  airfoil  surface.  Although  the  surface  of  the  element  is  somewhat 
uneven,  the  edges  are  clean  and  the  dots  can  be  stacked  with  little  compression  of  the 
lower  layers.  Careful  profilometer  measurements  indicate  a mean  thickness  of  6 //m 
per  layer.  The  practical  limit  on  stacking  the  dots  is  3 to  4 layers  based  on  the  ability 
to  maintain  a well-defined  edge.  For  the  present  experiment,  one-  and  three-layer  dots 
are  used  to  produce  k = 6 pm  and  18  //in  roughness.  Thicker  roughness  elements  are 
die-cut  from  3M  #850  industrial  polyester  tape.  The  diameter  of  these  elements  is 
also  3.7  mm,  and  the  total  height  including  the  adhesive  is  48  pm  (measured  with  a 
profilometer).  When  stamped  carefully  from  the  #850  tape,  these  circular  roughness 
elements  have  exceptionally  clean  edges  and  an  absolutely  uniform  thickness. 

The  dots  are  applied  in  full-span  arrays  along  the  z axis  at  x/c  = 0.023.  This 
location  is  near  the  neutral  point  for  the  stationary  crossflow  instability,  and  has  been 
shown  by  Radeztsky  et.  al.  (1993a)  to  maximize  the  influence  of  the  roughness.  The 
fundamental  disturbance  wavelength  is  fixed  by  the  spanwise  spacing  of  the  elements. 
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Two  different  spacings  are  used  for  this  experiment:  12  mm  and  36  mm.  As  shown 
in  chapter  5,  the  12  mm  spacing  is  chosen  to  amplify  the  dominant  wavelength  that 
appears  in  the  absence  of  artificial  roughness.  The  36  mm  spacing  allows  the  study 
of  the  interaction  between  the  fundamental  and  multiple  harmonic  disturbances. 

The  size  of  the  elements  with  respect  to  the  boundary  layer  is  measured  by  the 
roughness  Reynolds  number: 


Rek 


U(k)k 

v 


(3.1) 


where  k is  the  dimensional  roughness  height  and  U(k)  is  the  total  boundary-layer 
velocity  at  the  top  of  the  element.  Clearly,  Rek  depends  on  the  freestream  veloc- 
ity and  the  roughness  location.  Since  the  boundary  layer  is  too  thin  to  measure 
at,  x/c  = 0.023,  the  theoretical  boundary-layer  solution  is  used  to  calculate  Rek. 
Table  3.1  lists  the  Rek  values  for  the  conditions  of  this  experiment.  Also  listed  is 
the  roughness  height  normalized  by  the  boundary-layer  thickness  <5  and  displace- 
ment thickness  S*  at  x/c  = 0.023.  The  Rek  values  are  well  below  the  Braslow  limit 
for  three-dimensional  roughness  (von  Doenhoff  and  Braslow  1961;  Juillen  and  Arnal 
1990),  hence  the  elements  do  not  trip  the  boundary  layer  or  induce  a local  turbulent 
wedge. 


3.4  Coordinate  Systems 

The  correct  interpretation  of  the  theoretical  and  experimental  results  requires  an 
understanding  of  the  coordinate  systems  in  which  the  data  are  presented.  Figures  3.5 
and  3.6  show  the  common  coordinate  systems  used  to  describe  a swept  wing.  The 
freestream  flow  is  from  left  to  right  as  indicated.  The  global  test- section  coordinates 
(AM',  Z)  are  aligned  with  the  tunnel.  The  streamwise  coordinate  A'  is  in  the  flow 
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Table  3.1:  Roughness  element,  measures  at  x/c.  - 0.023. 


Roughness  Configuration 

Type 

Layers 

k [//m] 

k/S 


k/8* 


Rek 


Rec  = 1.6  x 106 
(<5  = 0.92  mm,  5*  - 0.30  mm) 


GS-104 

1 

6 

0.0065 

0.0203 

0.061 

GS-104 

3 

18 

0.0196 

0.0610 

0.55 

3M  850 

1 

48 

0.0522 

0.1627 

3.8 

Rec  = 2.4  x 106 

(5  = 0.75  mm,  8*  = 0.24 

mm) 

GS-104 

1 

6 

0.0080 

0.0249 

0.11 

GS-104 

3 

18 

0.0240 

0.0749 

1.0 

3M  850 

1 

48 

0.0640 

0.1992 

7.0 

Rec  = 3.2  x 106 

(5  = 0.65  mm,  8*  = 0.21 

mm) 

GS-104 

1 

6 

0.0092 

0.0288 

0.17 

GS-104 

3 

18 

0.0277 

0.0865 

1.5 

3M  850 

1 

48 

0.0738 

0.2308 

10.7 

direction,  Y is  normal  to  the  front  wall  of  the  test  section,  and  the  vertical  coordinate 
Z completes  the  right-handed  system.  The  velocity  components  in  this  system  are 
denoted  by  (u,  v,w).  The  model-oriented  coordinate  system  ( x,y,z ) is  attached  to 
the  wing-chord  plane.  The  x coordinate  is  normal  to  the  leading  edge,  y is  normal  to 
the  chord  line,  and  z is  parallel  to  the  leading  edge  in  the  swept  span  direction.  The 
velocity  components  in  this  system  are  ( un , vn,  iun).  The  third  system  is  the  boundary- 
layer  coordinates  (xt,  yt,  Zt).  In  this  system,  xt  is  tangent  to  the  inviscid  streamline,  tjt 
is  normal  to  the  model  surface,  and  zt  completes  the  right-handed  orthogonal  system. 
The  velocity  components  in  the  boundary-layer  system  are  (ut,  vt,  wt).  The  crossflow 
direction  is  defined  by  zt,  and  wt  is  the  crossflow  velocity  component.  Figure  3.7 
shows  a typical  swept-wing  boundary  layer  in  the  boundary-layer  coordinate  system. 
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With  A positive  as  shown  in  figure  3.5,  the  crossflow  component  is  negative  upstream 
of  the  pressure  minimum.  Body-intrinsic  coordinates  shown  in  figure  3.6, 

represent  yet  another  coordinate  system  often  used  in  computational  analyses.  This 
system  is  similar  to  model-oriented  coordinates,  except  the  origin  is  attached  to  the 
airfoil  surface  rather  than  the  wing-chord  plane.  The  £ coordinate  is  tangent  to  the 
model  surface  and  normal  to  the  leading  edge.  The  coordinate  £ = yt  is  normal  to 
the  surface. 

The  differences  in  these  coordinate  systems  must  he  considered  when  compar- 
ing theoretical  and  experimental  results.  Computations  are  typically  performed  in 
model-oriented  or  body-intrinsic  coordinates,  however,  the  axes  of  the  instrumenta- 
tion traverse  are  aligned  with  the  global  (X,  Y,  Z)  coordinates.  While  it  is  certainly 
possible  to  program  the  three-dimensional  traverse  to  move  in  another  coordinate  sys- 
tem, practical  limits  on  the  step  size  may  prevent  this.  For  instance,  a step  AC  normal 
to  the  airfoil  surface  decomposes  into  the  steps  A A'  = AC  sin(0)  and  AY  = AC  cos(0), 
where  0 is  the  angle  between  C and  ) (see  figure  3.6).  However,  over  most  of  the 
model  0 is  only  a few  degrees.  Thus,  a reasonably  small  AC  (say,  30  pm)  would 
require  a A A'  that  is  smaller  than  the  minimum  step  size  in  the  X direction.  On 
the  other  hand,  moves  in  the  z direction  pose  no  problems  since  the  step  sizes  are 
typically  much  larger  and  decompose  equally  into  the  A and  Z directions. 

Because  of  these  restrictions,  no  adjustments  are  made  to  the  Y motion  of  the 
traverse.  Consequently,  the  hot-wire  scanning  techniques  discussed  in  section  4.4  ac- 
quire data  in  the  (V',  z)  plane.  Within  this  plane,  the  boundary-layer  hot-wire  probe 
support  is  rotated  about  the  Z axis  by  an  angle  0 to  account  for  local  surface  cur- 
vature' effects.  The  hot-wire  probe  itself  must  then  be  rotated  about  its  longitudinal 
axis  by  an  angle  a so  that  the  tines  arc  equidistant  from  the  model  surface.  These 
rotations  are  shown  in  figure  3.8,  and  the  angles  3 and  a are  listed  in  table  3.2  for 
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Table  3.2:  Boundary-layer  hot-wire  probe  rotation  angles. 


all  measurement  locations.  The  hot  wire  measures  the  total  velocity  normal  to  the 
element  after  both  rotations.  Accurate  comparisons  with  theory  depend  strongly  on 
the  ability  of  the  computational  results  to  be  cast  in  this  field.  For  example,  if  the 
computations  are  generated  in  body-intrinsic  coordinates,  the  velocity  profile  in  the 
(£,  C)  plane  must  be  projected  onto  the  plane  normal  to  the  hot-wire  element  and  then 
interpolated  from  Q to  Y . Only  then  can  the  boundary-layer  profiles  and  disturbance 
measurements  be  consistently  compared  with  this  experiment.  The  computational 
results  to  which  the  experimental  data  are  compared  are  transformed  in  this  manner. 


3.5  Theoretical  Predictions 


This  section  describes  the  theoretical  basis  for  the  experimental  test  conditions.  The 
basic-state  boundary-layer  solutions  are  presented,  as  are  stability  calculations  for 
the  stationary  crossflow  disturbances. 

3.5.1  Basic  State 


The  first  step  in  a stability  analysis  is  to  compute  the  basic-state  boundary  layer. 
This  requires  appropriate  edge  conditions,  which  for  this  experiment  are  generated 
by  the  MCARF  code  (Stevens  et  al.  1971).  MCARF  computes  the  Cp  for  an  unswept 


airfoil  in  a duct,  hence  it  correctly  handles  the  fiat  front  and  rear  walls  of  the  test 
section.  The  effect  of  wing  sweep  is  accounted  for  by  transforming  the  MCARF 
pressure  distribution  into  a “three-dimensional”  or  “swept”  Cp  according  to 

CP3  = CP2  cos2 (A).  (3.2) 

This  represents  a renormalization  by  the  dynamic  pressure  in  the  freestream  (A") 
direction.  Figure  3.9  shows  the  swept  airfoil  contour  and  upper-surface  Cp  for  the 
NLF(2)-0415  at  a = —4°  and  8f  = 0°  in  the  Unsteady  Wind  Tunnel.  There  are  three 
differences  between  this  plot  and  figure  3.3. 

1.  The  airfoil  cross  section  is  viewed  in  the  (X,  y)  plane  and  nondimensionalized 
by  the  swept  chord  C,  hence  it  appears  thinner.  (In  the  coordinate  systems 
described  in  section  3.4,  x/c  and  X/C  are  equivalent  but  y/c  and  y/C  are  not.) 

2.  The  swept  CP3  defined  bv  equation  (3.2)  is  plotted. 

3.  The  pressure  distribution  is  computed  with  the  front  and  rear  test-section  walls 
in  place. 

The  primary  effect  of  the  wind-tunnel  walls  is  to  accelerate  the  flow  to  a slightly  lower 
minimum  pressure.  This  will  increase — although  not  dramatically — the  crossflow  dis- 
turbance growth  rates. 

The  MCARF  code  is  run  in  an  inviscid  mode,  i.e.,  no  corrections  for  displace- 
ment thickness  are  applied.  The  cross-sectional  area  of  the  Unsteady  Wind  Tunnel 
test  section  increases  by  3.7%  to  account  for  wall  and  model  boundary-layer  growth. 
Although  only  approximate,  Dagenhart  (1992)  has  shown  that  this  correction  is  suf- 
ficiently accurate  for  this  experiment. 
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Using  the  Cp  distribution  of  figure  3.9,  the  laminar  boundary  layers  are  computed 
with  the  boundary-layer  code  of  Haynes  (Haynes  and  Reed  1996).  This  code  is  de- 
signed for  infinite-span  swept  wings  and  employs  the  conical-flow  assumption.  The 
boundary-layer  profiles  at  several  chord  locations  for  Rec  = 2.4  x 10f)  are  shown  in 
figures  3.10-3.12.  The  velocity  component  tangent  to  the  inviscid  streamline  (fig- 
ure 3.10)  accelerates  continuously  under  the  presence  of  the  favorable  pressure  gradi- 
ent. The  crossflow  component  (figure  3.11)  shows  rapid  initial  growth  due  to  the  large 
pressure  gradient  near  the  leading  edge,  followed  by  a more  moderate  and  sustained 
development  in  the  mid-chord  region.  The  streamwise  velocity  component  measured 
by  the  hot  wires  is  plotted  in  figure  3.12.  These  boundary-layer  profiles  are  typical 
for  swept-wing  flows. 

The  crossflow  Reynolds  number  is  based  on  the  maximum  crossflow  velocity  and 
the  larger  of  the  two  heights  where  the  crossflow  velocity  is  10%  of  the  maximum, 
thus 


Rec{  = (3.3) 

Historically,  the  crossflow  Reynolds  number  has  been  used  for  transition  correlations 
based  purely  on  basic-state  boundary-layer  characteristics  (see,  for  example,  Reed 
and  Haynes  1994).  For  the  present  experiment,  Rec f varies  from  5.7  at  x/c  = 0.005 
to  270  at  x/c  = 0.60  for  Rec  = 2.4  x 106.  Figure  3.13  shows  the  streamwise  variation 
of  Rec f for  several  chord  Reynolds  numbers. 

3.5.2  Linear  Stability  Calculations 

The  linear  stability  behavior  of  the  above  boundary  layers  aids  in  the  choice  of  exper- 
imental test  conditions  and  provides  a benchmark  to  which  the  experimental  results 
are  compared.  Traditional  analyses  of  this  type  involve  solving  the  Orr-Sommerfeld 


equation  for  the  disturbance  growth  rates  and  normal  modes  (discussed  below).  Many 
codes  such  as  SALLY  (Srokowski  and  Orszag  1977)  and  COSAL  (Malik  1982)  are  gen- 
erally available  for  this  task.  More  recently,  the  parabolized  stability  equations  (PSE) 
have  become  very  popular  as  an  alternate  method  for  predicting  stability  behavior 
(see  the  discussion  in  section  1.3.3). 

Although  the  experimental  results  are  compared  to  both  linear  and  nonlinear  PSE 
calculations  in  chapter  5,  the  present  discussion  focuses  on  spatial  stability  analysis 
using  linear,  parallel  theory.  While  these  restrictions  certainly  limit  the  class  of  flows 
to  which  the  analysis  applies,  they  allow  us  to  effectively  illustrate  the  fundamental 
ideas  of  stability  theory.  Complete  reviews  of  this  topic  can  be  found  in  many  places. 
In  particular,  the  reader  is  referred  to  Mack  (1984),  Arnal  (1992,  1994),  and  Saric 
(1992b,  1994c)  for  detailed  discussions.  The  summary  presented  here  follows  that  of 
Saric  (1994c). 


The  analysis  begins  by  assuming  a parallel  basic  state  given  by 

U = U(y),  V = 0,  W = W{y),  (3.4) 

where  y is  normal  to  the  surface  and  U,  V,  and  W are  the  chordwise  (:?;),  wall-normal 
(//),  and  spanwise  (2)  velocity  components,  respectively.  This  assumption  immedi- 
ately forces  a local  analysis  in  that  the  stability  characteristics  at  each  streamwise 
location  determined  using  the  local  velocity  distributions — are  obtained  indepen- 
dently of  all  others. 


The  total  field  quantities  consist  of  the  basic  state  plus  small,  three-dimensional 
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disturbances  and  are  written  as 


'f-i 

+ 

II 

3 

(3.5a) 

v = v'(x,  y,  2,  t) 

(3.5b) 

w — W + w'(x,  y,  z,  t) 

(3.5c) 

p = P + p'(x,y,z,t) 

(3.5d) 

Equations  3.5  are  substituted  into  the  incompressible  Navier-Stokes  equations,  the 
basic-state  solution  is  removed  (the  basic  state  itself  satisfies  the  Navier-Stokes  equa- 
tions), and  products  of  the  small  disturbance  quantities  are  neglected  (e.g.,  u'u'x  < 
u').  The  resulting  linear  disturbance  equations  are 


ux  -\-  v'y  + w'z  — 0 (3-0) 

U/  + T Uj.  + IJyV  + \\  u'z  + px  — ^2  u!/R  — 0 (3.7) 

v[  + Uv'x  + Wv'z  + p'y  - V2v'/R  = 0 (3.8) 

w[  + Uw'x  + Wyv'  + Ww'z  + p'z  - V2w'/R  = 0 (3.9) 


where  subscripts  denote  partial  differentiation  and  the  equations  have  been  appropri- 
ately nondimensionalized  ( R is  the  Reynolds  number). 

Equations  (3.6)-  (3.9)  are  reduced  to  ordinary  differential  equations  with  the  in- 
troduction of  the  normal  mode 


<]'(■>:,  y,  z,  t)  = q(y)  ^*+4*-^.)  + c c (3 .10) 

where  q'  represents  any  one  of  the  disturbance  quantities  and  C.C.  stands  for  com- 
plex conjugate.  The  form  of  the  normal  mode  is  suggested  by  the  linearity  of  the 
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disturbance  equations  and  the  fact  that  the  coefficients  (i.e.,  basic-state  velocities) 
are  functions  only  of  y.  Equation  (3.10)  represents  a single  sinusoidal  wave  (i.e., 
one  Fourier  mode)  with  chordwise  wavenumber  a,  spanwise  wavenumber  ft,  and  fre- 
quency u).  The  amplitude  function  q(y)  is  complex  but  q'  is  real  by  the  inclusion  of 
the  complex  conjugate. 

Substituting  equation  (3.10)  into  equations  (3.6)“(3.9)  produces  a 6th-order  sys- 
tem of  ordinary  differential  equations  describing  the  disturbance  motion.  This  system 
can  be  combined  into  a single  4th-order  equation  known  as  the  Orr-Sommerfeld  equa- 
tion: 

| (/T  - k2f  - iR  [(at/  + ftW  - u>)  (J D 2 - k2)  - a ( D2U ) - ft  ( D*W )]  J 0 = 0, 

(3.11) 

where  k2  = cv"  + ft2,  D = d/dy,  and  0 = v represents  the  normal-mode  amplitude 
function  for  the  v'  disturbance.  (The  variable  change  is  simply  to  remain  consistent 
with  conventional  notation.)  The  disturbances  are  zero  at.  the  wall  and  must  vanish 
in  the  farfield,  hence  the  boundary  conditions  are 

(j)( 0)  = D<j)( 0)  = 0 and  0 0.  (3.12) 

The  condition  Dft  = 0 is  a statement  of  the  no-slip  condition  [Dv  = 0 in  equa- 
tion (3.6)]. 

Equations  (3.11)  and  (3.12)  form  a linear,  homogeneous  system  for  the  normal- 
mode amplitude  function  0 (=  v)  in  terms  of  the  basic-state  velocity  profiles  and  the 
parameters  a.  ft,  u>,  and  R.  The  system  defines  an  eigenvalue  problem  (by  virtue  of 
its  linearity  and  homogeneity),  thus  solutions  to  equation  (3.11)  are  obtained  only 
for  certain  combinations  of  the  parameters.  This  combination  is  expressed  bv  the 
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dispersion  relationship 


F{a,{3,ijj,R)  = 0.  (3.13) 

For  spatial  stability  analysis,  a and  /i  are  complex  and  uj  is  real.  Thus,  equa- 
tion (3.13)  represents  two  equations  in  six  unknowns.  Actually,  the  Reynolds  number 
is  known  and  the  frequency  is  typically  specified,  leaving  the  wavenumbers  ( ar  and  fir) 
and  the  growth  rates  (a,  and  j3i ) as  the  four  undetermined  quantities.  Consequently, 
additional  constraints  on  a and/or  3 are  required  in  order  to  solve  equation  (3.13). 
Many  such  constraints  have  been  proposed  and  are  reviewed  by  Arnal  (1994).  Typi- 
cally they  involve  a somewhat  ad  hoc  assumption  concerning  the  direction  of  growth, 
the  direction  of  propagation,  or  the  disturbance  wavelength.  The  present  analysis 
follows  Mack  (1988)  and  assumes,  quite  simply,  that  (3r  is  fixed  (and  specified)  and 
fti  = 0.  This  amounts  to  nothing  more  than  a statement  of  the  infinite-span  as- 
sumption and  is  verified  experimentally  in  chapter  6.  With  uj  = 0 for  stationary 
disturbances  and  3r  input,  equation  (3.13)  is  solved  for  the  chordwise  wavenumber 
olt  and  growth  rate  —a*.  The  solution  is  strictly  local,  and  must  be  re-evaluated  at 
each  chord  location  using  the  new  local  conditions  as  input. 

The  disturbance  amplitude  ratio  between  two  locations  is  computed  by  integrating 
the  spatial  growth  rate,  giving 


A 

T0 


= e 


N 


(3.14) 


where  the  amplification  factor  or  integrated  growth  rate  N (also  called  the  “A-factor” ) 
is  given  by 


N = 


a dx. 


(3.15) 
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In  general,  the  growth  rate  a is  some  combination  of  a,  and  $ meant  to  represent 
the  total  amplification  of  the  disturbance.  As  mentioned  above,  however,  # = 0 for 
the  infinite-span  airfoil.  In  this  case,  the  integrated  growth  rate  N reduces  to 

N = f —aldx  = In—.  (3.16) 

The  iV-factor  provides  the  basis  for  the  celebrated  eN  transition  correlation  method 
(Smith  and  Gamberoni  1956;  van  Ingen  1956).  See  Saric  (1992b)  and  Arnal  (1992, 
1994)  for  current  reviews. 

Using  these  techniques,  linear  stability  predictions  for  stationary  crossflow  distur- 
bances are  computed  for  the  experimental  test  conditions  outlined  in  section  3.3.  The 
calculations  are  generated  with  the  stability  code  of  Haynes  (Haynes  and  Reed  1996). 
Figure  3.14  shows  N versus  x/c  for  Rec  = 2.4  x 106.  Amplitude  curves  are  shown  for 
several  values  of  the  spanwise  wavelength  Xz  = 2rr / 1%.  The  .A- factors  are  computed 
relative  to  x/c  — 0.05.  The  data  show  typical  linear  behavior  for  stationary  crossflow 
modes  in  a swept-wing  boundary  layer.  Short  wavelength  disturbances  grow  early  in 
the  region  of  rapid  boundary-layer  growth,  while  long-wavelength  modes  are  amplified 
at  larger  x/c  after  a region  of  initial  decay.  The  modes  with  the  largest  -V-factor  are 
those  that  arc  short  enough  to  grow  early  but  long  enough  to  remain  unstable  until 
tin1  pressure  minimum  at  x/c  = 0.71.  In  agreement  with  the  SALLA  calculations 
used  by  Dagenhart  (1992)  and  Radeztskv  (1994),  the  A,  = 12  mm  mode  is  the  most 


unstable. 


CHAPTER  4 


Experimental  Methods 


The  experimental  procedures  at  the  Unsteady  Wind  Tunnel  are  designed  to  provide 
high-quality,  reliable  data  in  each  test.  This  requires  a clear  understanding  of  the 
stability  problem,  as  well  as  certain  wind-tunnel  effects  that  may  influence  the  results. 
This  chapter  discusses  these  issues,  and  describes  the  specific  data-acquisition  and 
post-processing  techniques  used  for  the  present  experiment. 

4.1  Special  Considerations 

4.1.1  Stationary  Disturbances 

The  stationary  crossflow  problem  on  a swept  wing  presents  several  unique  challenges 
to  the  experimentalist.  Unlike  the  rotating  disk  analogue,  the  stationary  waves  are 
fixed  in  space.  Thus,  the  measurement  probe  must  be  moved  through  the  flowfield  to 
capture  the  stationary  structure.  This  requires  special  steps  to  ensure  the  positioning 
accuracy  of  the  instrumentation  traverse  system,  and  corrective  actions  are  usually 
needed  to  compensate  for  misalignment  between  the  traverse  and  the  model. 

Stationary  disturbances  also  require  mean-flow,  or  DC,  measurements.  At  a given 
measurement  point  in  the  boundary  layer,  the  time-dependent  u'  fluctuations  are 
due  to  traveling  disturbances,  whereas  the  mean  of  the  measurement  represents  the 
stationary  component.  Consequently,  one  cannot  AC-couple  the  signal  and  amplify 
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the  fluctuations.  This  process  would  remove  the  stationary  component  of  interest  and 
leave  only  the  traveling  waves.  Instead,  the  stationary  disturbance  component  must 
be  determined  by  looking  for  spatial  variations  in  the  total  boundary-layer  velocity 
a.  This  works  well  downstream  where  the  distortions  of  the  mean  boundary  layer  are 
large,  however  near  the  leading  edge  it  becomes  increasingly  difficult  to  extract  small 
disturbance  quantities  from  0(1)  measurements.  Unfortunately,  the  critical  initial 
disturbance  amplitude  must  be  determined  under  these  conditions,  which  underscores 
the  need  for  accurate  measurements. 

4.1.2  Wind-Tunnel  Heating 

Additional  complications  are  caused  bv  the  fan-generated  and  frictional  heating  of 
a closed-circuit  wind  tunnel.  This  problem  can  be  minimized  or  even  removed  by 
the  introduction  of  a heat  exchanger.  These  systems,  however,  add  mechanical  com- 
plexities to  the  tunnel,  create  a large  drag  penalty,  and  may  alter  the  flow  quality 
(Rae  and  Pope  1984).  In  addition,  the  large  thermal  time  constant  characteristic 
of  large  wind  tunnels  may  affect  the  accuracy  to  which  the  test-section  temperature 
can  be  controlled.  Thus,  active  cooling  systems  are  not  necessarily  the  panacea  for 
wind-tunnel  heating. 

The  Unsteady  Wind  Tunnel  is  not  equipped  with  a cooling  system,  hence  the 
test-section  temperature  is  at  the  mercy  of  the  natural  heating  effects.  The  warm 
desert  environment  at  Arizona  State  University  compounds  the  problem.  The  ambi- 
ent temperature  frequently  reaches  45  °C  during  the  summer  months.  As  a result, 
the  test-section  temperature  can  increase  by  20  °C  during  a high-Reynolds-number 
experiment.  On  the  other  hand,  low-speed  runs  that  continue  into  the  evening  are 
often  characterized  by  an  initial  increase,  then  decrease  in  the  test-section  temper- 
ature. Velocity  measurements  using  hot-wire  anemometry  must  account  for  these 


temperature  fluctuations. 


4.2  Hot-Wire  Techniques 
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Experiments  at  the  Unsteady  Wind  Tunnel  use  hot-wire  anemometry  for  all  quanti- 
tative velocity  measurements  within  the  boundary  layer.  This  topic  is  well  studied, 
with  no  dearth  of  publications  describing  various  calibration  and  signal-analysis  tech- 
niques. The  approach  at  the  Unsteady  Wind  Tunnel  focuses  on  empirical  accuracy. 
Although  the  methods  account  for  the  dominant  behavior  of  the  hot  wire,  no  par- 
ticular physical  model  for  the  velocity  or  temperature  response  is  used.  Computer 
analysis  and  automated  procedures  are  emphasized.  This  eliminates  the  need  for 
analog  equipment  (such  as  temperature  compensators  and  linearizers),  but  requires 
additional  software  development. 

4.2.1  Velocity  Calibration 

For  low-speed  flows,  the  voltage  output  from  a constant-temperature  anemometer 
(CTA)  is  dominated  by  the  fluid  velocity  and  temperature  difference  between  the 
hot-wire  element  and  the  fluid  (Perry  1982).  Thus, 

U — T{E,  AT),  (4.1) 

where  the  velocity  is  isolated  as  the  dependent  variable  since  it  is  the  desired  quantity. 
Assume  for  now  the  temperature  difference  between  the  wire  and  the  fluid  (AT)  is 
fixed.  Then  the  classic  model  of  the  anemometer’s  velocity  response  is  given  by  King’s 
Law  (King  1914,  1915): 


U = (P  + QE2)2 , 


(4.2) 


where  P and  Q are  parameters  involving  the  physical  properties  of  the  hot-wire 
element  and  the  fluid.  The  Unsteady  Wind  Tunnel  model  generalizes  King’s  Law  by 
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using  a simple  polynomial  fit: 

U = j^AkEk.  (4.3) 

fc= o 

The  coefficients  Ak  are  determined  such  that  equation  (4.3)  best  fits  (in  the  least- 
squares  sense)  a moderate  number  of  voltage/ velocity  data  points.  A 4t.h-order  fit 
(n  = 4)  is  chosen  to  match  King’s  Law  to  leading  order.  The  polynomial  fit,  however, 
is  more  general  in  that  the  E and  E3  terms  not  present  in  King’s  Law  allow  a more 
robust  approximation  of  the  data  points. 

The  hot  wires  are  calibrated  in  situ  by  monitoring  the  anemometer  output  and 
Pitot-static  velocity  as  the  tunnel  is  increased  though  a set  of  predetermined  calibra- 
tion speeds.  Figure  4.1  shows  a typical  calibration  curve  fit.  The  calibration  points 
must  entirely  encompass  the  operating  speeds  of  the  hot  wire  since  extrapolating  a 
polynomial  fit  is  potentially  unstable.  The  procedure  is  fully  automated  and  takes  less 
than  5 minutes,  allowing  the  hot  wires  to  be  calibrated  at  virtually  any  time  during 
the  experiment.  As  a standard  practice,  the  calibration  is  repeated  each  morning  to 
protect  against  possible  “calibration  drift”  experienced  by  tungsten  hot  wires  (Perry 
1982). 

4.2.2  Temperature  Compensation 

The  above  discussion  has  neglected  the  effect  of  fluid  temperature  changes  on  the 
CTA  output.  This  cannot  be  ignored  since  the  test-section  temperature  varies  during 
an  experiment  (see  section  4.1.2). 

Numerous  methods  for  CTA  temperature  compensation  have  been  proposed.  The 
current  Unsteady  Wind  Tunnel  approach  follows  Bearman  (1971),  and  improves  the 
technique  of  Radeztsky  et  al.  (1993b)  to  include  the  velocity  dependence  of  the 
temperature-compensation  coefficient.  As  with  the  velocity  calibration,  the  emphasis 
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is  on  empirical  data  fits  and  computer  automation.  The  model  assumes  the  square  of 
the  CTA  voltage  varies  linearly  with  temperature  for  a fixed  velocity,  that  is 

E2c=E2  + Ct(Tc-T),  (4.4) 

where  Ec  is  the  CTA  “equivalent”  voltage  at  the  calibration  temperature  Tc.  The 
compensation  coefficient  Ct  is  a function  of  velocity,  as  indicated  in  figure  4.2.  These 
data  are  obtained  by  monitoring  the  test-section  temperature,  anemometer  output, 
and  Pitot-static  velocity  while  the  tunnel  is  pre-heated.  (The  pre-heat  also  mini- 
mizes temperature  changes  during  the  experiment.)  The  speed  dependence  of  the 
compensation  coefficient  is  modeled  well  with  a 2nd-order  polynomial  curve  fit, 

Cr  = E B,U\  (4.5) 

3=0 

whose  coefficients  Bj  are  determined  by  the  least-squares  technique.  It  is  important 
to  note  that  the  compensation  coefficients  are  obtained  before  the  velocity  calibration 
discussed  in  the  preceding  section.  This  allows  the  velocity  calibration  to  be  corrected 
for  any  temperature  changes,  effectively  producing  a constant-temperature  hot-wire 
calibration. 

For  data-acquisition  purposes,  equations  (4.5)  and  (4.3)  are  substituted  into  equa- 
tion (4.4)  to  give  an  implicit  relation  for  the  temperature-compensated  CTA  voltage: 

El  = E2  + f CT.-T).  (4.6) 

j= 0 L U=0  / 

Given  a temperature  T and  an  anemometer  voltage  E,  equation  (4.6)  is  implicitly 
solved  for  Ec.  This  corrected  voltage  is  then  converted  to  an  accurate  hot-wire  velocity 
using  equation  (4.3).  Figure  4.3  illustrates  the  accuracy  obtainable  with  this  method. 
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The  test-section  temperature  is  more  than  10  °C  above  the  calibration  temperature  for 
all  speeds,  and  the  uncorrected  hot-wire  velocities  are  plagued  with  errors  exceeding 
15%.  On  the  other  hand,  the  error  in  the  temperature-compensated  measurements  is 
less  than  1%,  over  the  entire  speed  range. 

The  calibration  and  data-acquisition  techniques  outlined  above  are  implemented 
through  a custom  set  of  programming  libraries  (Reibert  1996).  In  addition  to  stan- 
dardizing the  methods,  this  allows  acquisition  program  to  be  developed  and  updated 
efficiently. 


4.3  Wind-Tunnel  Speed  Control 

Wind-tunnel  heating  can  also  affect  the  freestream  velocity  during  an  experiment. 
More  precisely,  changes  in  kinematic  viscosity  caused  by  wind-tunnel  heating  (or 
cooling)  may  necessitate  a change  in  freestream  velocity  in  order  to  keep  a relevant 
nondiinensional  parameter  constant.  For  example,  consider  an  experiment  in  which 
the  test-section  temperature  rises  from  30  °C  to  50  °C  (not  uncommon  at  the  Unsteady 
Wind  Tunnel).  If  the  velocity  remains  fixed,  the  Reynolds  number  will  decrease  by 
11%  (air  at  1 atmosphere).  This  may  have  a significant  impact  on  the  stability 
characteristics  of  the  experiment.  Consequently,  the  freestream  velocity  must  be 
increased  to  compensate  for  the  increase  in  the  kinematic  viscosity. 

To  coordinate  these  velocity  adjustments,  the  Unsteady  Wind  Tunnel  is  controlled 
from  the  Sun  SPARCstation  data-acquisition  computer  using  a “cruise  control”  pro- 
gram (figure  4.4).  This  custom  LabVIEW  virtual  instrument  continually  monitors 
the  tunnel  conditions  (test-section  static  pressure,  dynamic  pressure,  temperature, 
and  fan/motor  rpm)  and  adjusts  the  motor  speed  to  maintain  the  desired  control 
value.  Three  control  parameters  are  available:  freestream  velocity,  Reynolds  num- 
ber, and  nondiinensional  frequency.  When  the  control  parameter  or  control  value 
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is  changed,  the  program  commands  the  data-acquisition  system  to  provide  continu- 
ous updates  as  tunnel  approaches  the  new  conditions.  When  the  difference  between 
the  measured  and  target  freestream  velocity  is  within  tolerance  (1%  of  the  target  or 
0.1  m/s,  whichever  is  smaller),  the  program  enters  a “slow  maintenance”  mode  in 
which  the  tunnel  conditions  are  updated  less  frequently  (typically  every  15  seconds). 
This  has  two  benefits.  First,  the  fan  motor  is  not  unnecessarily  chasing  continuous 
but  insignificant  rpm  adjustments.  Instead,  the  program  simply  applies  a “drift  cor- 
rection” to  compensate  for  the  slowly  varying  temperature.  Second,  the  A/D  system 
is  available  for  other  tasks.  This  is  essential  since  the  A/D  converters  are  a shared 
system  resource  and  must  be  available  for  experimental  data  acquisition. 

The  present  work  uses  the  Reynolds  number  control  parameter  to  maintain  a 
constant  chord  Reynolds  number  during  each  experiment. 

4.4  Disturbance  Measurement  Techniques 

The  majority  of  the  data  presented  in  chapters  5 and  6 are  obtained  using  the  two 
measurement  techniques  discussed  below.  Although  each  methods  is  optimized  for 
a specific  task,  the  common  goal  is  to  quantify  the  stationary  crossflow  disturbance 
amplitude  by  accurately  measuring  spatial  variations  in  the  mean  boundary-layer 
velocity. 

4.4.1  Boundary-Layer  Profile  Scans 

Multiple  wall-normal  boundary-layer  scans  provide  a detailed,  two-dimensional  map 
of  the  stationary  structure.  These  maps  are  constructed  by  taking  a spanwise  series 
of  mean-flow  boundary-layer  profiles  at  constant  x/c.  A typical  set  consists  of  100 
profiles,  each  separated  by  1 mm  in  the  swept  span  direction.  The  technique  is  similar 
to  that  of  Radeztsky  (1994). 
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Before  the  scan  begins,  the  “BL”  hot  wire  (positioned  at  the  end  of  the  sting,  see 
figure  2.4)  is  adjusted  for  local  surface  curvature  effects.  This  is  necessary  since  the 
hot  wires  are  aligned  to  measure  the  strearnwise  velocity  u and  not  the  component 
■un  normal  to  the  leading  edge.  The  adjustment  involves  moving  the  hot  wire  very 
close  to  the  model  and,  under  magnification,  rotating  the  probe  about  its  axis  so  that 
the  tines  are  equidistant  from  the  surface.  This  ensures  that  the  hot  wire  element  is 
not  averaging  across  the  boundary  layer.  The  probe  is  then  moved  to  the  edge  of  the 
boundary  layer  and  the  scan  is  started. 

Each  individual  profile  is  fully  self  contained.  The  BL  hot  wire  measures  the 
boundary-layer  velocity  and  the  “FS”  hot  wire  (located  mid-span  on  the  sting,  see 
figure  2.4)  tracks  the  external-flow  velocity  at  the  same  x/cJ  The  scan  starts  with 
the  BL  hot  wire  positioned  at  the  edge  of  the  boundary  layer.  This  reference  point  is 
acquired,  and  the  boundary-layer-edge- to-external-flow  velocity  ratio  is  recorded.  The 
probes  are  then  stepped  toward  the  model,  acquiring  the  mean  velocity  from  both  hot 
wires  at  each  point.  The  boundary-layer  velocity  is  normalized  by  the  instantaneous 
edge  velocity  l’e.  This  is  computed  by  scaling  the  FS  hot-wire  velocity  bv  the  edge-to- 
external-flow  velocity  ratio  from  the  initial  reference  point.  This  instantaneous  local 
normalization  is  essential  since  the  tunnel  speed  is,  in  general,  changing  due  to  heating 
('fleets  (see  section  4.3).  As  the  scan  progresses,  the  step  size  is  scaled  with  u/Ue  to 
provide  finer  resolution  near  the  airfoil  surface.  When  u/Ue  reaches  a predetermined 
threshold,  the  profile  is  terminated  and  the  hot  wires  are  moved  to  the  boundary-layer 
edge  at  the  next  span  location.  A new  boundary-layer-edge-to-external-flow  velocity 
ratio  is  obtained,  and  the  entire  procedure  is  repeated.  Figure  4.5  shows  a typical 
profile. 


'The  FS  hot  wire  measures  the  velocity  in  the  external  flowfield  approximately  150  mm  “above” 
the  boundary-layer  probe. 


Once  the  profiles  are  acquired,  the  data  produce  a two-dimensional  map  of  the 
mean-flow  boundary  layer,  complete  with  any  amplified  stationary  disturbance  struc- 
ture. From  this  velocity  field  it  is  relatively  straightforward  to  generate  velocity 
contours,  disturbance  profiles,  and  stationary  wave  mode  shapes  from  which  the  dis- 
turbance amplitude  is  obtained. 

Like  all  experimental  measurement  techniques,  this  method  has  both  advantages 
and  disadvantages.  On  the  plus  side,  the  individual  profiles  are  self-consistent  and  self- 
aligning. The  location  of  the  airfoil  surface  is  determined  by  extrapolating  each  profile 
to  zero  velocity.  The  slight  pressure-gradient-induced  curvature  of  the  boundary 
layer  is  ignored,  and  a straight  line  is  fit  through  the  lower  portion  of  the  profile. 
Dagenhart  (1992)  and  Radeztsky  (1994)  successfully  used  this  technique  in  earlier 
swept- wing  experiments  at  ASU.  With  a known  reference  point  for  each  profile,  the 
scans  are  easily  assembled  to  produce  the  two-dimensional  velocity  field.  Moreover, 
traverse/model  alignment  concerns  virtually  disappear  since  each  profile  is  aligned 
with  the  airfoil  surface  independently. 

The  foremost  disadvantage  of  this  method  is  the  time  required  for  a complete 
set  of  profiles.  The  0(1)  spanwise  gradients  make  it  impossible  to  use  the  highly 
optimized  technique  of  Radeztsky  (1994),  hence  each  individual  boundary-layer  profile 
is  obtained  without  assistance  from  the  previous  scan.  Moreover,  the  strong  distortion 
of  the  boundary-layer  flow  demands  high  resolution  in  the  wall-normal  direction.  A 
typical  profile  in  the  mid-chord  region  (<5  ~ 4 mm)  contains  on  average  60  data  points. 
Each  data  point,  in  turn,  requires  approximately  3 seconds  to  acquire  (2  seconds  of 
sampling  at  1 kHz  and  1 second  for  the  traverse  move  and  overhead).  Thus,  5 hours 
are  needed  to  obtain  100  profiles.  Consequently,  this  technique  can  quickly  become 
prohibitively  time  consuming. 
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4.4.2  Spanwise  Scans  at  Constant  Y 


One  of  the  primary  goals  of  this  experiment  is  to  accurately  determine  the  wavenum- 
ber content  of  the  amplified  stationary  crossflow  disturbances.  This  will  not  only 
provide  single-wavelength  growth  rates  for  comparison  with  linear  theory,  but  will 
also  guide  nonlinear  PSE  computations  by  providing  initial  conditions  for  individual 
modes.  To  quantify  this  analysis,  some  sort  of  spectral  decomposition  must  be  ap- 
plied to  the  data.  The  spectral  method,  in  turn,  will  impose  certain  restrictions  on 
the  data  that  will  ultimately  require  a second  measurement  technique. 

Spectral  Considerations 

As  with  any  experiment  that  involves  analysis  in  the  spectral  domain,  frequency 
resolution  issues  are  of  paramount  importance.  For  discrete  experimental  data,  the 
sampling  parameters  must  be  chosen  appropriately  to  yield  the  desired  spectral  char- 
acteristics. The  present  experiment  adds  a slight  twist:  the  acquired  signal  is  a 
function  of  space  rather  than  time.  In  this  context,  inverse  wavelength  takes  on  the 
role  of  "frequency”,  and  the  sampling  parameters  determine  the  wavelength  resolu- 
tion in  the  spatial  power  spectrum.  The  translation  from  the  time  domain  to  the 
spatial  domain  illuminates  the  requirements  that  must  be  satisfied  by  the  measure- 
ment technique.  This  is  briefly  reviewed  below.  A complete  discussion  of  discrete 
spectral  methods  can  be  found  in  many  sources.  In  particular,  Press  et  al.  (1992) 
give  a succinct  explanation  of  several  techniques. 

In  general,  a discrete  signal  contains  N evenly  spaced  measurements  over  a length 
S (in  time  or  space).  The  sampling  interval  As  is 


A,  = 


(4.7) 


The  sampling  frequency  / = 1/AS.  For  a time-domain  signal,  / is  in  Hertz.  If,  on 
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the  other  hand,  S (and  therefore  A.,)  are  spatial  quantities,  / is  in  cycles  per  length. 
This  corresponds  to  inverse  wavelength  (not  wavenumber). 

The  one-sided  power  spectral  density  (PSD)  is  defined  at  the  N/2  + 1 discrete 
frequencies 


fk  = 


k — 0 


(4.8) 


giving  a frequency  resolution  of  A/  = 1 /(NAS).  The  Nyquist  critical  frequency 
is  fc  = 1/(2AS).  For  a spatial  spectrum,  the  wavelength  resolution  AA  is  more 
important.  Since  A = 1//, 


AA 


d A 
df 


A2  A2 
NAS  ~ ~S' 


(4.9) 


Several  key  features  of  the  spatial  spectrum  are  exposed  by  equation  (4.9). 


1.  The  wavelength  resolution  is  a function  of  A,  and  increases  with  the  square  of  the 
wavelength.  Thus,  small  wavelengths  are  resolved  better  than  large  wavelengths. 


2.  Increasing  the  length  of  the  data  set  S decreases  the  wavelength  resolution. 
Thus,  all  wavelengths  are  best  resolved  by  maximizing  the  spatial  extent  of 
the  sample.  Increasing  the  number  of  samples  N without  changing  S does  not 
affect  AA;  the  “extra”  information  goes  into  increasing  the  Nyquist  frequency 
by  reducing  the  sampling  interval  A,,. 


3.  The  smallest  wavelength  for  which  a power  estimate  is  obtained  is  A = 2A.,. 
This  is  used  to  determine  the  sampling  interval  for  the  measurements. 

In  short,  one  must  sample  longer  to  improve  the  spectral  frequency  resolution.  Sam- 
pling faster  only  serves  to  increase  the  Nyquist  frequency. 
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Another  important  issue  is  the  choice  of  kernel  used  to  generate  the  PSD.  The 
Fourier  transform  is  certainly  the  most  common,  due  in  no  small  part  to  the  advent 
of  the  fast  Fourier  transform  (FFT).  However,  many  alternate  methods  are  avail- 
able, each  with  characteristic  advantages  and  disadvantages  (Kay  and  Marple  1981). 
Previous  ASU  swept-wing  experiments  (Radeztsky  et  al.  1994)  used  the  maximum 
entropy  method  (MEM),  which  generates  the  PSD  in  terms  of  a Laurent  series.  The 
result  is  a continuous  spectrum  that  is  w'ell-suited  for  “peaky”  spectra  and  small  data 
sets.  The  major  drawback,  however,  is  that  the  user  must  (somewhat  arbitrarily) 
choose  the  number  of  poles  in  the  Laurent  series  expansion.  If  too  few  are  chosen, 
certain  features  in  the  spectrum  may  be  suppressed.  On  the  other  hand,  too  many 
poles  can  create  spurious  peaks  or  even  split  sharp  peaks.  In  addition,  spectral  peaks 
can  shift  frequency  depending  on  the  phase.  As  pointed  out  by  Radeztsky  (1994), 
practical  application  of  the  MEM  requires  subjective  interaction  from  the  user  and  is 
therefore  difficult  to  automate. 

To  avoid  these  problems,  the  present  analysis  relies  strictly  on  FFT-based  spectral 
methods.  The  subjectivity  of  the  MEM  is  removed,  allowing  the  PSD  generation  to  be 
automated.  However,  a new  set  of  concerns  arises  as  a result  of  the  discrete  spectrum: 

1.  How  much  confidence  does  one  have  in  the  FFT-based  power  estimates,  i.e.,  is 
the  PSD  estimate  for  each  frequency  “bin”  correct  on  average? 

2.  Can  the  frequency  (or  wavelength)  resolution  be  improved  when  the  length  of 
the  data  series  is  limited? 

In  general,  these  are  nontrivial  issues  that  are  discussed  in-depth  in  the  literature 
(see,  for  example,  Press  (4.  al.  1992  and  Kay  and  Marple  1981).  Here  the  topics  will 
be  briefly  addressed  in  the  context  of  the  present  experiment. 

The  first  problem  is  cast  in  terms  of  spectral  leakage  and  variance.  Leakage  (of 
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energy  into  adjacent  frequency  “bins”)  is  caused  by  the  abrupt  turning  on  and  off  of 
the  data  series.  The  solution  to  this  problem  is  to  window  the  data  with  a transfer 
function  that  changes  more  gradually  from  zero  to  a maximum  and  back  to  zero  over 
the  length  of  the  data.  Many  standard  windows  are  commonly  available  (Press  et  al. 
1992).  The  Welch  window  is  used  for  this  experiment.  This  quadratic  window  is 
computationally  efficient  and  provides  a reasonable  trade-off  between  the  sidelobe 
fall-off  and  the  highest  sidelobe  level.  The  variance  issue  is  handled  by  averaging. 
This  process  involves  splitting  the  data  series  into  smaller  segments,  computing  the 
PSD  for  each  segment,  then  averaging  the  individual  spectra  to  form  the  final  power 
spectrum.  This  presents  little  difficulty  in  the  time  domain,  where  it  is  easy  to  acquire 
multiple  back-to-back  data  segments  by  simply  sampling  longer.  For  spatial  spectra, 
however,  the  length  of  the  sample  is  typically  fixed  by  the  physical  constraints  of 
the  system.  In  this  case,  one  can  repeat  the  experiment  multiple  times,  but  this  is 
often  too  time  consuming  due  to  overhead  associated  with  setting  up  each  run.  A 
better  approach  is  to  increase  the  sampling  frequency  by  a factor  M representing  the 
desired  number  of  averages.  The  data  stream  can  then  be  demultiplexed  into  M data 
segments — each  covering  the  total  length  S — which  provide  the  necessary  averages. 
This  technique  preserves  the  desired  wavelength  resolution  since  the  spatial  extent  of 
each  data  segment  is  not  compromised. 

The  second  point  concerns  improving  the  wavelength  (or  frequency)  resolution 
without  increasing  the  sample  length  S.  This  is  particularly  important  for  spatial 
spectra  due  to  the  physical  limits  of  the  measurement  region.  At  issue  here  is  whether 
or  not  the  data  or  the  FFT  can  be  modified  so  that  the  PSD  estimates  better  ap- 
proximate the  “true”  power  spectrum  of  the  (windowed)  data.  The  answer  is,  quite 
simply,  yes  and  no.  The  pessimist  will  point  out  that  once  the  data  are  acquired,  the 
“information  content”  of  the  signal  is  fixed  and  cannot  be  magically  increased  by  the 
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FFT  procedure.  In  particular,  the  fundamental  frequency  resolution  is  entirely  de- 
termined by  the  number  of  samples  and  the  sampling  interval,  and  cannot  be  altered 
by  the  FFT.  However,  all  hope  is  not  lost.  It  is  entirely  possible  to  increase  the  num- 
ber of  discrete  spectral  modes  used  in  the  Fourier  expansion.  This  does  not  improve 
the  fundamental  frequency  resolution,  but  instead  allows  the  FFT  to  interpolate  ad- 
ditional PSD  estimates  to  better  approximate  the  available  information  contained 
within  the  data  (Kay  and  Marple  1981).  This  is  accomplished  by  zero-padding  each 
data  segment  before  it  is  converted  by  the  FFT.  Care  must  be  taken  to  add  the  zeros 
after  the  window  has  been  applied,  else  the  effect  of  the  window  will  be  seriously 
compromised.  When  applied  correctly,  zero-padding  has  the  quite  desirable  effect  of 
“smoothing''  an  otherwise  coarse  power  spectrum. 

All  spectra  used  in  this  work  have  mean  square  power  spectral  density  plotted 
on  the  ordinate.  Thus,  by  Parseval’s  theorem  the  square  root  of  the  integral  of  the 
spectrum  equals  the  rms  of  the  original  signal. 

Scanning  Technique 

Spatial  spectra  for  stationary  crossflow  waves  are  best  obtained  from  a single  span- 
wise  trace  of  u/Ue  at  constant  Y.  The  sampling  requirements  for  these  measurements 
are  obtained  directly  from  the  previous  discussion.  The  velocity  profile  must  cover 
the  entire  spanwise  extent  of  the  measurement  region  in  order  to  provide  the  best 
possible  wavelength  resolution  in  the  spectral  domain.  For  the  present  experiment, 
that  distance  is  240  mm.  With  256  data  points  along  this  span,  the  wavelength  res- 
olution will  be  A A — 0.6  mm  at  A 12  mm,  and  wavelengths  between  1.9  mm  and 
240  mm  will  be  resolved  by  the  spectrum.  Four-times  over-sampling  (acquiring  1024 
points  along  the  total  span)  allows  four  averages,  each  with  256  points  covering  the 
total  span  of  240  mm. 

One  way  to  generate  the  u/Ue  versus  span  profile  is  to  take  a spanwise  slice 
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(at  constant  Y)  across  a set  of  wall-normal  boundary-layer  profiles.  This  has  the 
distinct  advantage  that  all  heights  within  the  boundary  layer  are  available  for  analysis. 
Unfortunately,  at  nearly  3 minutes  per  boundary  layer  it  would  take  over  two  days  of 
continuous  data  acquisition  to  obtain  the  1024  profiles  needed  to  meet  the  spectral 
requirements!  Consequently,  this  approach  is  abandoned. 

To  reduce  the  data-acquisition  time,  the  measurements  are  limited  to  a single 
spanwise  scan  at  a constant  height  above  the  airfoil  surface.  As  with  the  wall-normal 
profiles,  the  spanwise  scan  begins  with  the  BL  hot  wire  positioned  at  the  edge  of 
the  boundary  layer.  This  reference  point  is  acquired,  and  the  boundary-layer-edge- 
to-external-flow  velocity  ratio  is  recorded.  The  probes  are  then  stepped  toward  the 
model  until  the  BL  hot  wire  is  a user-specified  distance  from  the  surface.  The  scan 
proceeds  with  the  probes  being  moved  in  the  swept  span  (2)  direction,  acquiring  the 
mean  velocity  at  1024  points  along  the  240  mm  span.  The  boundary-layer  velocity 
is  nondimensionalized  with  the  instantaneous  edge  velocity  in  the  same  manner  as 
the  wall-normal  scans.  At  each  measurement  point.,  the  V'  position  of  the  hot  wires 
is  adjusted  to  compensate  for  any  misalignment  between  the  model  and  the  traverse 
(more  on  this  later).  Figure  4.6  shows  a typical  scan.  This  velocity  profile  is  well- 
suited  for  the  spatial  spectral  techniques  discussed  above,  and  the  peaks  in  the  PSD 
are  integrated  to  obtain  the  disturbance  energy  in  individual  crossflow  modes. 

Traverse  Alignment 

The  stationary  crossflow  disturbance  amplitude  is  a strong  function  of  Y.  It  is 
therefore  critical  to  remain  a fixed  distance  above  the  airfoil  surface  during  a spanwise* 
scan.  This  is  accomplished  by  prefacing  the  spanwise  scan  with  a traverse-alignment 
scan.  The  purpose  of  this  preliminary  scan  is  to  locate  the  airfoil  surface  in  the 
coordinate  system  in  which  the  traverse  moves. 

In  previous  ASU  experiments,  Radeztsky  (1994)  aligned  the  traverse  by  moving 
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along  the  span  and  recording  the  Y position  (in  traverse  coordinates)  corresponding 
to  u/Ue  = 0.75.  This  method  is  relatively  fast  and  works  well  provided  T = Y{ujUe ) 
is  one-to-one  at  all  span  locations.  The  primary  disadvantage  is  that  the  physical 
location  of  the  model  surface  is  never  actually  determined.  Only  a relative  correc- 
tion is  obtained  allowing  the  spanwise  scan  to  proceed  at  a nominal  u/Ue.  For  the 
present  experiment,  disturbance  amplitudes  computed  from  the  spanwise  scans  will 
be  directly  compared  to  those  obtained  from  wall-normal  boundary-layer  profiles. 
The  quality  of  these  comparisons  depends  strongly  on  being  able  to  specify  the  exact 
height  above  the  airfoil  surface  for  the  spanwise  scan.  Moreover,  Y = Y(u/Ue ) is 
not  necessarily  one-to-one  due  to  the  large  distortions  of  the  mean  boundary  layer. 
Consequently,  the  technique  of  Radeztsky  cannot  be  used. 

The  location  of  the  model  surface  can  be  accurately  determined,  however,  by 
extrapolating  a boundary-layer  profile  to  zero  velocity.  This  idea — borrowed  from 
the  wall-normal  scanning  technique — provides  the  basis  for  the  traverse-alignment 
procedure.  A small  number  of  wall-normal  boundary-layer  profiles  (usually  25)  are 
taken  over  the  240  mm  extent  of  the  spanwise  scan.  These  profiles  are  optimized  to 
concentrate  data  points  near  the  airfoil  surface  at.  the  expense  of  detail  high  in  the 
boundary  layer.  This  minimizes  the  time  required  for  each  profile  while  maintaining 
an  accurate  extrapolation  to  zero  velocity.  The  surface  coordinate  is  recorded  at  each 
span  location,  and  the  data  are  fit  with  a low-order  polynomial.  Figure  4.7  shows 
a typical  traverse-alignment  profile.  The  curve  fit.  effectively  represents  the  physical 
location  of  the  model  surface  as  a function  of  span.  This  relation  is  used  by  the 
spanwise  scan  to  maintain  the  desired  height  above  the  airfoil  surface,  resulting  in 
the  velocity  profile  shown  in  figure  4.6. 
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4.5  Flow  Visualization 

Stationary  crossflow  waves,  through  a spanwise  modulation  of  the  mean  boundary- 
layer  flow,  produce  a surface  shear  stress  pattern  that  can  be  identified  with  a shear- 
stress-sensitive  flow-visualization  technique.  Naphthalene  is  an  effective  medium  for 
this  type  of  visualization;  at  room  temperature  it  sublimes  at  a rate  proportional  to 
shear  stress. 

To  apply  the  naphthalene  to  the  wing,  it  is  first  dissolved  in  1 , 1 , 1-trichloroethane. 
The  saturated  solution  is  sprayed  onto  the  airfoil  surface  using  a compressed-air 
sprayer.  The  solvent  quickly  evaporates,  leaving  a thin  coating  of  naphthalene  on 
the  model.  Care  must  be  taken  to  avoid  spraying  the  leading  edge  since  the  naphtha- 
lene adds  roughness  to  the  surface.  As  a general  rule,  the  naphthalene  is  not  applied 
forward  of  .T fc  = 0.20.  According  to  Radeztsky  et  al.  (1993a),  this  is  well  beyond  t he 
point  where  surface  roughness  influences  the  stationary  crossflow  instability. 

With  the  naphthalene  coating  on  the  airfoil,  the  wrind  tunnel  is  started  and  brought 
to  the  desired  test  condition.  The  naphthalene  quickly  sublimes  in  regions  of  high 
shear,  but  remains  where  the  shear  stress  is  low.  The  stationary  crossflow  waves 
are  visible  as  streaks  of  naphthalene.  This  pattern  is  caused  by  the  transposition  of 
high-  and  low-momentum  fluid  within  the  boundary  layer  resulting  in  the  spanwise 
modulation  of  the  surface  shear  stress.  The  naphthalene  also  quickly  sublimes  in  the 
turbulent  region,  leaving  a clear  picture  of  the  transition  front. 


CHAPTER  5 


Results — Part  1:  Baseline  Configuration 


5.1  Overview 

The  experimental  results  and  discussion  are  divided  into  two  chapters,  of  which  this 
is  the  first.  Stationary  crossflow  disturbance  measurements  for  the  baseline  test  con- 
dition (described  below)  are  presented  in  detail.  For  this  data  set,  comparisons  with 
theoretical  predictions  provide  valuable  assistance  not  only  in  analyzing  the  data, 
but  in  understanding  the  underlying  physical  mechanisms  of  the  crossflow  instability. 
Transition  data  and  basic-state  measurements  are  also  presented  in  this  chapter. 

The  following  chapter  continues  the  investigation  by  describing  the  experimental 
results  under  varying  freestream  arid  surface-roughness  conditions.  From  the  reference 
condition  discussed  in  this  chapter,  the  chord  Reynolds  number,  roughness  spacing, 
and  roughness  height  are  varied  independently  to  cover  the  parameter  space.  Chap- 
ter 6 concludes  by  addressing  several  summary  topics  and  specific  questions  raised 
during  the  experiment. 

In  both  this  and  the  next  chapter,  all  of  the  crossflow'  disturbance  amplitude 
measurements  are  acquired  using  the  two  hot-wire  scanning  techniques  described  in 
section  4.4.  The  results  are  grouped  into  the  six  data  sets  summarized  in  table  5.1. 
The  chord  Reynolds  number  is  shown  in  the  column  marked  Re,.,  while  the  roughness 
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Table  5.1:  Experimental  data  set  configuration. 


Data 

Set 

ReJ  10° 

Roughness 

Scan  Type 

k [/im] 

A2  [mm 

BL 

Span 

A 

2.4 

6 

12 

• 

• 

B 

1.6 

6 

12 

• 

C 

3.2 

6 

12 

• 

• 

V 

2.4 

6 

36 

• 

• 

£ 

2.4 

18 

12 

• 

• 

T 

2.4 

48 

12 

• 

• 

height  and  spanwise  spacing  are  given  in  the  k and  A,  columns,  respectively.  Data  are 
obtained  using  the  wall-normal  boundary-layer  profile  and  spanwise  scan  at  constant 
Y techniques  as  indicated  by  a bullet  (•)  in  the  BL  or  Span  column.  This  chapter 
concentrates  on  the  baseline  or  control  configuration  defined  by  data  set  A.  Data 
sets  B through  T are  presented  in  chapter  6. 

The  roughness  elements  are  applied  in  full-span  arrays  at  x/c.  = 0.023  as  out- 
lined in  section  3.3.3.  Different  roughness  configurations  are  achieved  by  varying  the 
height  or  spanwise  spacing  of  the  elements,  but  the  chord  location  remains  fixed  and 
the  roughness  always  covers  the  entire  span.  Therefore,  a statement  of  the  element 
height  and  spanwise  spacing  unambiguously  defines  the  roughness  configuration.  The 
shorthand  notation  [A;|A2]  will  be  used  for  this  purpose.  Thus,  [6|36]  denotes  k = 6 /in l 
roughness  with  a 36  mm  spanwise  spacing,  while  [48|12]  means  the  roughness  elements 
are  48  /mi  thick  and  12  mm  apart.  The  absence  of  artificial  roughness  is  denoted  by 
[010]- 

Following  the  convention  commonly  used  in  CFD,  crossflow  modes  are  expressed 
as  (/.  m).  In  this  notation,  / is  the  disturbance  frequency  (zero  for  stationary  waves) 
and  m is  the  mode  number  defined  as  the  disturbance  wavenumber  normalized  by 
the  wavenumber  of  the  fundamental  mode.  Thus,  m = 1 denotes  the  fundamental 
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disturbance,  while  m — 2 represents  the  first  superharmonic  with  a wavenumber  equal 
to  twice  that  of  the  fundamental  (i.e.,  half  the  wavelength).  The  spanwise-invariant 
disturbance — usually  called  the  “mean-flow  distortion”  mode — is  denoted  by  (0,0). 
This  should  not  be  confused  with  the  local  distortions  of  the  mean  boundary  layer 
caused  by  the  stationary  crossflow'  waves.  Where  ambiguity  in  the  terminology  may 
result,  the  meaning  should  be  clear  from  the  context. 

5.2  Flow  Visualization 

Surface  shear-stress  patterns  and  transition  locations  are  obtained  using  naphthalene 
flow'  visualization  (see  section  4.5).  This  technique  has  been  calibrated  against  hot- 
wire and  hot-film  measurements  (Dagenhart  1992)  and  has  been  successfully  used  by 
Dagenhart  et  al.  (1989,  1990)  and  Radeztskv  et  al.  (1993a)  as  a transition  indicator. 
One  of  the  more  significant  results  to  come  out  of  Radeztskv  et.  al.  (1993a)  was  the  fact 
that  micron-sized  roughness  is  not  effective  beyond  10%  chord.  Thus,  the  inherent 
roughness  associated  with  the  deposition  of  naphthalene  is  benign  if  applied  away 
from  the  attachment  line. 

Figures  5.1  and  5.2  show  the  naphthalene  pattern  at  Rec  = 2.2  x 106  and  3.2  x 10(i. 
The  roughness  configuration  is  [6)12]  for  both  cases.  The  flow'  is  left  to  right,  and  the 
photographs  show  the  airfoil  upper  surface  as  viewed  through  a glass  window  in  the 
test-section  front  wall.  Lines  of  constant  chord  drawn  on  the  model  with  a felt-tipped 
pen  appear  as  black  lines  parallel  to  the  leading  edge.  The  white  numerals  indicate 
x/c  in  percent.  The  stationary  crossflow'  vortex  pattern  is  clearly  indicated  by  the 
alternating  light  and  dark  streaks.  Strictly  speaking,  these  streaks  show'  the  spanwise 
modulation  of  the  surface  shear  stress  caused  by  the  transposition  of  high-  and  low- 
momentum  fluid  within  the  boundary  layer  (more  on  this  later).  In  this  sense,  the 
naphthalene  pattern  provides  a visual  image  of  the  “footprint”  of  the  stationary  wave 
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structure. 

The  region  within  which  the  boundary  layer  is  turbulent  is  marked  by  the  absence 
of  naphthalene.1  The  jagged  transition  front  that  develops  as  a series  of  turbulent 
wedges  is  typical  for  swept-wing  boundary  layers  dominated  by  the  stationary  cross- 
flow instability.  It  is  interesting  to  note  that  the  turbulent  wedge  pattern  is  somewhat 
nonuniform  is  span  even  when  the  initial  conditions  are  carefully  controlled  with  ar- 
tificial roughness.  This  feature  becomes  more  pronounced  as  the  growth  rates  are 
decreased  (i.e.,  the  transition  pattern  is  more  regular  for  Rer  — 3.2  x 106  than  for 
Re,,.  = 2.2  x 10°).  Without  a detailed  receptivity  study,  it  is  unclear  whether  this 
nonuniformity  is  caused  by  submicron  imperfections  in  the  height  or  edges  of  tin' 
roughness  elements  or  is  due  to  small  perturbations  in  the  spanwise  spacing  (or  both). 

Some  light  is  shed  on  this  phenomenon  by  repeating  the  flow  visualizations  with  a 
different  roughness  spacing.  Figures  5.3  and  5.4  show  the  same  experiments  with  the 
roughness  spacing  increased  to  36  mm  (i.e.,  [6|36]  roughness).  As  with  the  12  mm 
spacing,  the  turbulent  wedge  pattern  is  more  uniform  for  the  larger  Reynolds  number. 
However,  for  both  Re,  = 2.2  xlO6  and  3.2  xlO6,  the  jagged  transition  front  is  strikingly 
more  regular  for  the  36  mm  roughness  spacing  than  for  the  12  mm  spacing.  Clearly, 
the  spectral  content  of  the  surface  roughness  distribution  can  significantly  influence 
the  local  transition  characteristics  of  the  boundary  layer.  It  is  also  evident  that  the 
36  mm  spacing  excites  multiple  disturbance  modes.  In  both  figures  5.3  and  5.4, 
the  turbulent  wedges  are  36  mm  apart  corresponding  to  the  roughness  spacing,  but 
the  naphthalene  clearly  shows  vortex  streaks  at,  larger  spanwise  wavenumbers  (i.e., 
smaller  wavelengths) . 

Although  the  turbulent  wedges  that  comprise  the  transition  front  are  well-defined, 

1 This  does  not  apply  to  the  leading-edge  region,  where  the  naphthalene  sublimes  quickly  due  to 
the  high  shear  of  the  thin,  laminar  boundary  layer.  Moreover,  naphthalene  is  not  applied  upstream 
of  x/c  0.20  to  avoid  adding  surface  roughness  in  this  region. 
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Table  5.2:  Approximate  transition  location  determined  by  naphthalene  flow  visual- 
ization. 


Roughness 

Transition 
Location  [x/c] 

A:  [//m] 

X2  [mm] 

6 

12 

0.52 

0.32 

6 

36 

0.49 

0.30 

18 

12 

0.51 

0.30 

48 

12 

0.49 

0.28 

Rec 

2.4  x 10fi 

3.2  x 106 

the  jagged  pattern  and  potential  spanwise  nonuniformity  create  some  ambiguity  in 
determining  the  transition  location.  Dagenhart  (1992)  defined  the  transition  location 
as  the  average  of  the  beginning  and  ending  locations  of  the  turbulent  wedges.  This  was 
necessitated  by  the  random  distribution  of  “natural”  surface  roughness  that  influenced 
transition  in  his  experiments.  The  present  work  ignores  the  downstream  edge  of  the 
saw-tooth  pattern,  and  defines  the  transition  location  as  the  average  starting  location 
of  the  turbulent  wedges.  This  definition  is  based  on  the  simple  premise  that  the 
onset  of  transition  signifies  the  loss  of  laminar  flow.  In  addition,  considering  only  the 
vertices  of  the  wedges  helps  remove  some  of  the  ambiguity  when  the  transition  front 
is  nonuniform  in  span.  With  this  in  mind,  the  average  transition  location  for  several 
test  conditions  is  shown  in  table  5.2.  As  indicated  in  figures  5.1  5.4,  increasing  the 
roughness  spacing  moves  transition  forward . On  the  other  hand,  fixing  the  spanwise 
spacing  at  12  mm  and  increasing  the  roughness  height  from  6 //m  to  48  /an  (an  Re * 
increase  of  nearly  two  orders  of  magnitude)  has  surprising  lit  tle  effect  on  the  transition 
location.  This  is  in  contrast  to  the  results  of  Radeztsky  et  al.  (1993a),  which  show 
that  the  vertex  of  the  transition  wedge  behind  an  isolated  roughness  element  moves 
forward  by  nearly  10%  chord  when  the  roughness  height  is  increased  from  6 //in  to 
18  //m. 
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5.3  Basic  State 

The  basic  state  is  documented  with  pressure  measurements  and  boundary-layer  pro- 
files. The  former  provide  information  on  the  global  flowfield,  while  the  latter  allow 
specific  details  of  the  basic-state  boundary  layer  to  be  analyzed. 

5.3.1  Cp  Measurements 

The  pressure  distribution  is  measured  with  two  rows  of  pressure  ports  mounted  in 
the  test  surface  of  the  NLF(2)-0415.  Each  row  consists  of  23  pressure  taps  aligned 
with  the  freestream  (Ar)  direction.  The  ports  extend  from  x/c  = 0 to  x/c  = 0.78 
and  are  clustered  near  the  leading  edge.  One  row  of  ports  is  located  near  the  ceiling 
or  “upper”  end  of  the  model,  while  the  other  is  located  near  the  floor  or  “lowrer” 
end.  The  ports  are  647  mm  (measured  parallel  to  the  leading  edge)  from  the  center 
span  of  the  model.  The  pressure  coefficient  is  determined  by  measuring  p — with 
a differential  pressure  transducer  and  dividing  by  the  freestream  dynamic  pressure 
pU^/2.  This  normalization  generates  the  swept  Cp ( defined  by  equation  (3.2). 

Figures  5.5  5.7  show  the  pressure  distribution  obtained  from  both  sets  of  pressure 
ports  for  Rt'.(.  = 1.6  x 10°,  2.4  x 10G,  and  3.2  x 106.  Also  plotted  is  the  theoretical  Cv 
computed  using  the  MCARF  code  (Stevens  et  al.  1971).  As  discussed  in  section  3.5.1, 
the  computations  account  for  the  front  and  rear  test-section  walls  but  do  not  apply 
any  displacement  thickness  corrections.  The  MCARF  solution  under-predicts  the 
measured  Cv  for  all  Reynolds  numbers.  For  the  upper  ports,  the  difference  is  largest 
in  the  range  0.05  < x/c  < 0.35.  Dagenhart.  (1992)  indicates  that  the  test-section 
ceiling  liner  is  slightly  too  thin  in  this  region,  to  which  he  attributes  the  difference 
in  pressure.  For  the  lower  ports,  the  difference  between  the  measured  and  theoretical 
Cp  increases  with  increasing  x/c.  This  may  be  caused  by  the  inclined  floor  of  the  test 
section,  which  drops  by  50  mm  over  its  length  (4.9  m)  to  approximately  account  for 
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wall  boundary-layer  growth  (see  section  3.5.1).  The  proximity  of  the  pressure  ports 
to  the  test-section  end  liners  also  compromises  the  comparison  with  MCARF.  All 
of  the  ports  are  within  50-  150  mm  of  the  liner  surface,  which  magnifies  the  effect  of 
small  errors  in  the  end  liners  with  respect  to  their  influence  in  the  mid-span  region. 
Boundary-layer  profiles  presented  in  the  following  section  confirm  that  the  basic  state 
is  spanwise  independent  in  the  test  region  of  the  model. 

The  pressure  distribution  has  a slight  Reynolds-number  dependence  as  shown  in 
figures  5.8  and  5.9.  This  is  not  a compressibility  effect  (M  <0.1  for  all  operating 
speeds),  but  instead  is  due  to  changes  in  the  displacement  thickness  of  the  boundary 
layers  on  the  test-section  walls  and  the  airfoil  model.  This  is  not  unexpected  since 
the  inclination  of  the  test-section  floor  is  technically  optimum  for  only  one  operating 
condition.  The  Reynolds-number  effect,  however,  is  weak  and  can  be  ignored  without 
loss  of  accuracy. 

In  short,  although  there  are  slight  differences  between  the  measured  and  theoret- 
ical pressure  distribution,  the  generally  good  agreement  shows  that  the  experimental 
fiowfield  reasonably  matches  the  CFD  design. 

5.3.2  Boundary-Layer  Profiles 

As  discussed  in  chapter  3,  the  45°-swept  NLF(2)-0415  airfoil  at  a = —4°  produces 
large-amplitude  stationary  crossflow  waves  for  moderate  chord  Reynolds  numbers. 
While  ideal  for  the  study  of  this  instability  mechanism,  the  strong  distortions  of 
the  mean  boundary  layer  make  it  impossible  to  measure  the  basic-state  profiles  for 
most  of  the  test  conditions  of  this  experiment.  Even  in  the  absence  of  artificial 
surface  roughness,  the  nonlinear  mean-flow  distortion  [(0,0)  mode]  is  observed  for 
Rec  < 1.8  x 10’.  In  light  of  this,  basic-state  boundary  layers  are  measured  oniv 
at  Rec  = 1.6  x 106.  At  this  Reynolds  number,  the  stationary  wave  growth  rates 
are  sufficiently  small  and  the  mean-flow  distortion  is  negligible.  Although  this  does 
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not  allow  a direct  comparison  with  theory  for  the  majority  of  the  test  cases,  it  does 
provide  a reference  point  at  which  the  agreement  between  theory  and  experiment  can 
be  assessed. 

Basic-state  profiles  are  obtained  with  the  boundary-layer  scanning  technique  de- 
scribed in  section  4.4.1.  Figure  5.10  shows  the  profiles  at  xjc  = 0.20.  The  plot 
actually  contains  100  profiles  each  separated  by  1 mm  in  the  swept  span  (z)  direc- 
tion. Each  profile  contains  approximately  45  measurement  points  in  the  wall-normal 
direction.  The  profiles  are  virtually  identical,  and  show  no  evidence  of  stationary 
crossflow  waves.  More  importantly,  the  absence  of  any  variation  confirms  that  the 
flowfield  is  spanwise  invariant  in  the  test  region  of  the  model.  Thus,  the  infinite-span 
assumption  is  verified.  This  plot  also  answers  two  questions  about  the  measurement 
technique.  First,  the  boundary  layer  is  sufficiently  linear  near  the  wall  to  allow  a 
straight-line  extrapolation  to  locate  the  airfoil  surface.  Second,  because  the  airfoil 
surface  is  located  by  each  profile  independently,  the  method  is  “self-aligning”  in  that 
the  scans  are  well-matched  at  the  wall  and  the  boundary-layer  edge. 

The  boundary-layer  profiles  acquired  under  the  same  conditions  at  xjc  = 0.60 
are  plotted  in  figure  5.11.  These  100  mean-flow  profiles  show  a very  slight  spanwise 
modulation  due  to  the  presence  of  weak  stationary  crossflow  waves.  Because  the 
spanwise  variation  is  small  (less  than  1%  rrris),  the  average  of  the  100  profiles  is 
taken  as  the  basic-state  boundary  layer.  In  fact,  at  all  chord  locations  the  basic-state 
profile  is  defined  as  the  spanwise  average  of  100  mean-flow  boundary-layer  scans  taken 
over  99  mm  of  span.  These  profiles  are  shown  in  figure  5.12  for  0.05  < xjc  < 0.60. 
The  number  of  measurement  points  in  the  wall-normal  direction  varies  with  xjc,  but 
ranges  from  about  30  data  points  at  xjc  — 0.05  to  more  than  60  at  xjc  — 0.60. 

The  measured  profiles  are  compared  to  the  theoretical  basic-state  solutions  in 
figures  5.13  5.24.  The  computational  results  are  provided  by  Haynes  using  the 


boundary-layer  code  of  Haynes  and  Reed  (1996).  It  is  important  to  note  that  the 
dots  are  used  simply  to  differentiate  between  the  experimental  and  computed  pro- 
files. The  experimental  data  are  actually  obtained  with  much  higher  resolution  than 
indicated  by  the  symbols.  The  theoretical  profiles  are  projected  onto  the  experimen- 
tal measurement  plane  (as  discussed  in  section  3.4)  to  allow  an  accurate  comparison. 
Displacement  thickness,  momentum  thickness,  and  shape  factor  comparisons  (fig- 
ures 5.25-5.27)  indicate  that  the  computed  profile  is  slightly  decelerated  with  respect 
to  the  measured  boundary  layer  at  all  chord  locations.  This  effect  may  be  caused  by 
any  number  of  factors,  ranging  from  experimental  flow  misalignment  to  the  break- 
down of  the  boundary-layer  equations  near  the  leading  edge. 

Similar  basic-state  differences  were  noted  in  the  previous  ASU  experiments  in- 
volving the  NLF(2)-0415  model  (Dagenhart  1992;  Radeztsky  1994).  In  these  cases, 
an  ad  hoc  adjustment  to  the  angle  of  attack  and/or  Reynolds  number  was  applied  in 
order  to  bring  the  theoretical  basic  state  in  line  with  the  experimental  data.  While 
it  can  be  argued  that  these  adjustments  may  increase  the  validity  of  disturbance- 
amplitude  comparisons,  the  present  work  makes  no  attempt  to  “fix”  the  differences 
in  the  computed  and  measured  basic-state  profiles  for  the  following  reasons. 

1.  Due  to  experimental  constraints,  any  flowfield  adjustments  would  be  based 
solely  on  comparisons  of  the  streamwise  velocity  u.  In  terms  of  the  cross- 
flow  velocity  component,  however,  the  validity  of  these  adjustments  cannot  be 
ascertained. 

2.  The  Cp  measurements  indicate  that  the  experimental  flowfield  agrees  well  with 
the  CFD  design.  The  difference  in  the  level  of  the  Cp  (see  figure  5.6)  could  be 
due  to  a slight  angle-of-attack  misalignment,  however  any  adjustment  to  the 
angle  of  attack  or  sweep  angle  is  unwarranted  in  light  of  the  good  agreement 
between  the  measured  and  computed  pressure  gradient. 
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3.  The  NLF(2)-0415  at  a — —4°  generates  almost  no  lift  when  placed  in  the 
Unsteady  Wind  Tunnel  test  section.  Consequently,  the  upstream  or  forward 
influence  of  the  model  is  minimized.  Flowfield  calculations  indicate  that  the 
test-section  velocity  at  the  location  of  the  Pitot-static  probe  is  within  0.5%  of 
the  velocity  in  the  farfield.  Thus,  a Reynolds-number  adjustment  would  be 
insignificant. 

Simply  stated,  although  comparisons  of  the  measured  and  theoretical  boundary-layer 
profiles  indicate  a slight  disagreement  in  the  basic  state,  there  is  simply  no  experi- 
mental data  available  that  conclusively  justify  a contrived  adjustment  to  the  flowfield 
parameters.  Moreover,  whereas  linear  stability  calculations  are  sensitive  to  small 
changes  in  the  basic  state,  the  strong  nonlinear  characteristics  of  this  problem  ap- 
pear to  make  the  actual  stability  behavior  more  forgiving  of  basic-state  discrepancies. 
This  will  be  highlighted  in  the  following  sections,  where  disturbance  amplitudes  com- 
puted with  nonlinear  parabolized  stability  equations  agree  remarkably  well  with  the 
experimental  results. 

5.4  Natural  Surface  Roughness 

From  the  conditions  of  the  basic-state  scans  discussed  above,  one  needs  only  to  in- 
crease the  Reynolds  number  to  generate  measurable  distortions  of  the  mean  bound- 
ary layer.  Figure  5.28  is  a contour  plot  of  the  normalized  boundary-layer  velocity 
at  x./c  = 0.45  for  Rec  = 3.0  x 10G.  No  artificial  roughness  is  placed  on  the  airfoil 
surface.  The  figure  shows  the  streamwise  velocity  ujUe  in  the  (T,  z)  plane.  (The 
Y axis  is  magnified  by  a factor  of  10  with  respect  to  the  z axis  in  order  to  provide 
more  detail.)  The  flow  is  toward  the  reader  (i.e.,  the  observer  is  looking  upstream 
into  the  oncoming  boundary-layer  flow),  and  the  stationary  vortices  are  turning  in 
the  right-handed  sense.  The  velocity  contours  are  constructed  from  100  mean-flow 


boundary-layer  profiles  acquired  using  the  technique  described  in  section  4.4.1.  It 
is  important  to  note  that  the  wave-like  structure  of  figure  5.28  represents  the  inte- 
grated effect  of  the  weak  stationary  vortices  on  the  streamwise  velocity.  The  vortices 
themselves  are  co-rotating  (V,  w')  disturbances  that  are  too  weak  to  measure  directly. 

In  the  absence  of  artificial  surface  roughness,  the  naturally  occurring  stationary 
crossflow  waves  are  nonuniform  in  span.  Considering  the  near  mirror  finish  of  the 
airfoil  surface  (see  figure  3.4),  this  underscores  the  extreme  sensitivity  to  roughness- 
induced  initial  conditions  provided  by  submicron  surface  irregularities  near  the  leading 
edge. 

Figure  5.28  displays  a dominant  feature  at  a 12  mm  spanwise  spacing,  which  is 
approximately  the  most  amplified  stationary  wavelength  according  to  linear  theory. 
At  the  same  time,  the  richness  in  the  spectral  content  is  evident  and  indicates  non- 
linear interaction  among  many  modes.  This  is  typical  of  all  the  early  experiments 
(Muller  and  Bippes  1989;  Dagenhart  et  al.  1989,  1990;  Bippes  and  Nitschke-Kowsky 
1990;  Bippes  et  al.  1991),  and  leads  to  two  undesirable  side  effects. 

1.  The  spanwise  nonuniformity  creates  some  ambiguity  in  determining  the  dis- 
turbance amplitude.  For  example,  Dagenhart  (1992)  measured  the  disturbance 
amplitude  by  acquiring  data  over  a single  vortex  wavelength.  Clearly  the  growth 
rates  obtained  with  this  method  depend  strongly  on  which  vortex  is  chosen  (i.e., 
which  wave  in  figure  5.28),  as  well  as  the  ability  to  accurately  track  the  same 
vortex  at  all  chord  locations. 

2.  Even  a nonlinear  stability  calculation  that  includes  only  a single  mode  would 
be  inappropriate  to  characterize  the  disturbance  motion.  Instead,  the  compu- 
tations would  need  to  include  the  potentially  infinite  number  of  modes  that  are 
excited  by  the  random  surface  roughness. 


Thus,  the  unknown  natural  roughness  of  the  airfoil  surface  renders  the  comparison 
with  theory  unnecessarily  difficult. 

5.5  Distributed  Roughness 

In  order  to  generate  spanwise-uniform  stationary  crossflow  waves  with  a fixed  spectral 
content,  initial  conditions  are  controlled  with  artificial  surface  roughness  as  outlined 
in  section  3.3.3.  The  roughness  elements  are  distributed  in  full-span  arrays  at  xjc  = 
0.023.  This  location  near  the  first  neutral  point  of  the  stationary  instability— has 
been  shown  to  maximize  the  influence  of  the  roughness  (Radeztsky  et  al.  1993a). 
Applying  the  elements  along  the  entire  span  produces  a disturbance  field  that  is 
invariant  along  lines  of  constant  chord,  i.e.,  the  infinite-span  assumption  is  preserved. 

This  section  presents  the  results  for  the  baseline  configuration  defined  by  data 
set  A in  table  5.1.  The  roughness  height  is  6 /mi  and  the  spamvise  spacing  be- 
tween the  elements  is  12  mm  (i.e.,  [6)  1 2]  roughness).  This  spacing  is  chosen  to  excite 
the  dominant  stationary  wavelength  that  appears  in  the  absence  of  artificial  surface 
roughness.  The  data  are  acquired  at  Rec  = 2.4  x 10fi.  The  roughness  Reynolds 
number  Rek  =0.1  for  this  configuration. 

5.5.1  Vortex  Structure  and  Mean-Flow  Distortion 

Figure  5.29  shows  the  streamwise  velocity  contours  at  xjc  = 0.45  for  the  conditions 
described  above.  The  primary  difference  between  this  plot  and  figure  5.28  is  the  ad- 
dition of  the  artificial  roughness.  (The  only  other  difference  is  the  lower  Reynolds 
number,  which  was  decreased  because  of  the  increase  in  the  disturbance  strength.) 
The  dominance  and  uniformity  of  the  Az  = 12  mm  mode  is  striking.  The  distur- 
bance amplitude  is  well  defined  under  these  conditions  (as  will  be  demonstrated  in 
th(>  following  section),  and  meaningful  comparisons  with  both  linear  and  nonlinear 
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theoretical  predictions  are  possible.  The  dramatic  change  in  the  disturbance  mo- 
tion resulting  from  such  small  roughness  continues  to  highlight  the  sensitivity  of  the 
stationary  waves  to  initial  conditions  provided  by  leading-edge  surface  roughness. 

The  streamwise  velocity  contours  due  to  a single  stationary  vortex  are  isolated 
in  figure  5.30  and  plotted  on  a 1:1  scale.  The  integrated  effect  of  the  combined 
(v'.  w')  motion  of  the  crossflow  vortex  produces  regions  of  upwelling  and  downwelling 
that  transport  low-  and  high-momentum  fluid  within  the  boundary  layer.  The  asym- 
metry of  the  co-rotating  vortices  distorts  this  momentum  transfer  giving  the  ap- 
parent rollover  of  low-momentum  fluid  that  appears  above  high-momentum  fluid. 
The  v'  and  w'  components  are  actually  very  weak,  but  because  the  vortices  are  sta- 
tionary and  nearly  aligned  with  the  potential  flow  direction,  they  act.  on  the  same 
streamwise  fluid  to  produce  the  0(1)  v!  distortions.  This  process  is  described  in 
figures  5.31-5.33,  which  schematically  display  the  {v'.w')  motion  of  two  stationary 
vortices  above  the  measured  streamwise  velocity  contours.  The  measurements  are 
obtained  at  Rec  = 2.4  x 106  with  [6 j 12]  roughness.  One  can  clearly  see  how  the  u' 
distortion  develops  under  the  continuous  presence  of  the  vortices. 

Figure  5.34  shows  the  100  mean  boundary-layer  profiles  from  which  figure  5.29 
is  generated.  These  profiles  are  obtained  at  1 mm  intervals  in  the  swept  span  (z) 
direction  using  the  measurement  technique  described  in  section  4.4.1.  The  dots  mark 
the  spanwise  average  of  the  profiles,  which  accounts  for  the  basic  state  plus  the 
(0,0)  mean-flow  distortion  mode.  Each  individual  profile  contains  approximately 
60  measurement  points.  Each  measurement  point,  in  turn,  represents  a two-second 
average  velocity.  It  is  worth  emphasizing  that  these  time-averaged  profiles  represent  a 
spanwise  modulation  of  the  mean  flow  and  not  an  unsteady  oscillation  in  the  boundary 
layer.  One  can  clearly  see  how  the  momentum  transfer  caused  by  the  stationary 
vortices  has  distorted  the  mean  boundary  layer,  resulting  in  accelerated,  decelerated. 
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and  doubly-inflected  profiles  existing  millimeters  apart.  The  presence  of  the  mean- 
flow  distortion  mode  (indicated  by  the  distortion  of  the  spanwise-averaged  profile) 
shows  the  effect  of  strong  nonlinear  interactions. 

An  important  but  often  overlooked  aspect  of  this  instability  is  dramatically  il- 
lustrated by  figure  5.34.  Quite  simply,  the  mean  boundary  layer  looks  nothing  like 
the  undisturbed  basic  state.  Even  the  spanwise-averaged  profile  is  distorted  and  in- 
flectional due  to  strong  nonlinearities.  This  is  a fundamental  consequence  of  the 
stationary  nature  of  the  disturbance.  Although  the  crossflow  vortex  itself  is  a weak 
(</.  w')  motion,  its  integrated  effect  produces  an  0(1)  u'  disturbance  resulting  in 
a spanwise-varving  mean  flow  subject  to  different  secondary  instability  mechanisms 
depending  on  the  local  characteristics  of  the  boundary  layer.  Consequently,  it  is  phys- 
ically incorrect  to  model  the  stationary  disturbance  as  the  zero-frequency  limit  of  a 
traveling  wave  (although,  of  course,  this  is  mathematically  correct  with  the  frame- 
work of  linear  theory).  By  its  very  definition,  the  traveling  wave  cannot  produce  the 
integrated  effect  that  allows  a weak  disturbance  to  cause  the  strong  distortion  of  the 
streamwise  boundary  layer.  In  light  of  this,  the  continued  failure  of  linear  theory  to 
accurately  model  the  disturbance  growth  (Dagenhart  et  al.  1989,  1990;  Muller  and 
Bippes  1989;  Bippes  et  al.  1991;  Radeztsky  et  al.  1994)  should  not  be  unexpected. 
By  ignoring  the  mean-flow  distortion,  linear  theory  is  computing  the  stability  char- 
acteristics of  boundary  layers  that  do  not  exist.  These  ideas  will  be  revisited  later 
in  this  chapter  when  the  experimental  disturbance  amplitude  is  compared  to  various 
theoretical  predictions. 

5.5.2  Disturbance  Profiles  and  Mode  Shape 

Disturbance  profiles  are  generated  by  subtracting  the  basic  state  plus  mean-flow  dis- 
tortion (i.e.,  the  spanwise-averaged  profile)  from  each  of  the  100  individual  boundary- 
layer  profiles.  These  are  shown  in  figure  5.35  for  the  boundary-layer  profiles  of  fig- 
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ure  5.34.  The  wave-like  nature  of  the  spatial  disturbance  is  indicated  by  the  smooth 
phase  change  between  the  profiles.  The  distortions  reach  a maximum  of  36%  near 
Y — 1 mm  and  are  asymmetric  due  to  the  rollover  phenomenon  displayed  in  fig- 
ure 5.30.  The  disturbance  velocity  is  viewed  as  a contour  plot  in  figure  5.36.  (As 
with  the  total  velocity  contours,  the  Y axis  is  magnified  by  a factor  of  10  to  provide 
more  detail.)  Only  two  contour  levels  are  plotted  to  highlight  the  surplus  and  deficit 
in  the  streamwise  boundary-layer  velocity.  It  should  be  emphasized  that  the  con- 
tours do  not  represent  the  vorticity  of  the  stationary  vortices.  This  data-presentation 
format  accentuates  the  rollover  of  the  low-  and  high-momentum  fluid  and  explains 
the  streaks  in  the  naphthalene  flow  visualization,  but  otherwise  adds  no  information 
that  cannot  be  ascertained  from  the  total  velocity  contours  and  disturbance  profiles. 
Consequently,  disturbance  velocity  contours  will  be  presented  only  for  select  cases. 

The  stationary  wave  mode  shape  is  obtained  by  computing  the  spanwise  rms  of 
the  100  disturbance  profiles  at  each  Y location  in  the  boundary  layer.  Since  the 
stationary  disturbance  creates  a spatially  varying  wave,  this  is  equivalent  to  a time- 
domain  rms  computed  for  a traveling  fluctuation.  Figure  5.37  shows  the  mode  shape 
for  the  disturbance  profiles  of  figure  5.35.  The  rms  amplitude  reaches  a maximum  of 
19%  at  Y = 0.9  mm,  decays  to  a second  local  maximum  of  9%  at  Y — 2.7  mm,  then 
smoothly  vanishes  at  the  boundary-layer  edge.  The  second  lobe  occurring  high  in  the 
boundary  layer  corresponds  to  the  inflectional  distortion  of  the  spanwise-averaged 
profile  in  figure  5.34,  and  indicates  the  presence  of  nonlinear  effects.  This  points  out 
an  important  but  subtle  distinction  between  the  stationary  wave  inode  shape  and  a 
classic  linear  theory  eigenmode.  In  the  presence  of  nonlinearities,  the  disturbance 
mode  profile  will  contain  all  amplified  disturbance  modes  [except  the  (0,0)  mean- 
flow  distortion  mode],  whereas  the  eigenfunction  from  linear  theory  always  contains  a 
single  disturbance  mode.  This  means  the  stationary  wave  mode  shape  as  calculated 
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here  represents  the  amplitude  of  the  total  disturbance. 

5.5.3  Total  Disturbance  Amplitude  Distribution 

The  measurements  described  above  for  x jc  = 0.45  are  obtained  at  5%-chord  intervals 
from  x/c  = 0.05  to  x/c  = 0.45.  Figures  5.38-5.45  show  the  streamwise  velocity 
contours  for  the  additional  measurement  locations.  As  with  the  velocity  contour 
at  x/c  = 0.45  (figure  5.29),  the  Y-to-z  aspect  ratio  of  10:1  is  chosen  to  enhance 
the  detail  in  the  plots.  The  mean  boundary-layer  profiles  from  which  the  contour 
plots  are  constructed  are  shown  in  figures  5.46  5.53.  Figures  5.54  5.61  display  the 
corresponding  disturbance  profiles,  which  are  plotted  on  a consistent  scale  of  ±40% 
to  accentuate  the  disturbance  growth.  The  mode-shape  profiles  are  presented  in 
figure  5.62  for  0.10  < x/c  < 0.45.  This  series  of  plots  clearly  shows  the  development 
of  the  mean  boundary-layer  distortion  caused  by  the  stationary  crossflow  waves.  The 
disturbance  is  too  weak  to  measure  near  the  leading  edge,  but  bv  10%  chord  a very 
small  distortion  of  the  mean  boundary  layer  is  detected  (figure  5.39).  This  distortion 
grows  quickly  up  to  x/c  = 0.25,  where  the  disturbance  profiles  (shown  in  figure  5.58) 
begin  to  exhibit  the  asymmetric  shape  characteristic  of  the  downstream  locations. 
The  beginning  of  the  “second  lobe”  distortion  of  the  mode  shape  is  apparent  at 
x/c  = 0.30  (1'  « 1.6  mm),  confirming  the  presence  of  early  nonlinear  effects.  The 
nonlinear  features  continue  to  develop  to  the  final  measurement  location  of  x/c  = 0.45. 

The  data  at  x/c  — 0.05  highlight  the  accuracy  of  the  measurement  technique 
as  well  as  the  difficulty  in  acquiring  data  in  a thin  boundary  layer.  The  complete 
set  of  100  boundary-layer  scans  (figure  5.46)  contains  approximately  3500  distinct 
measurement  points  covering  99  mm  of  span  within  a 1.2  mm-thick  boundary  layer. 
Each  measurement  point,  in  turn,  represents  the  time  average  of  2000  instantaneous 
analog-to-digital  conversions  for  each  of  five  different  quantities  (test-section  static 
pressure,  temperature,  dynamic  pressure,  and  two  hot  wires).  The  entire  scan  takes 
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three  hours  to  complete.  The  straight-line  extrapolation  to  find  the  airfoil  surface  has 
worked  remarkably  well,  resulting  in  profiles  that  are  consistently  matched  at  the  wall 
and  the  boundary-layer  edge.2  In  fact,  all  100  profiles  are  buried  beneath  the  dots 
of  the  spanwise-averaged  profile.  The  method  is  not  perfect,  however,  and  the  large 
dimensional  du/dY  gradient  of  the  thin  boundary  layer  magnifies  any  small  errors  in 
the  traverse  movement.  Measurements  presented  later  will  show  that  the  stationary- 
wave-induced  distortions  of  the  mean  boundary  layer  are  too  small  to  be  extracted 
from  the  background  noise  under  these  conditions.  Therefore,  the  disturbance  profiles 
and  mode  shape  at  x/c  = 0.05  are  rejected  for  this  data  set  ( [6|  1 2]  roughness  at 
Rec  = 2.4  x 10*’),  and  x/c  = 0.10  is  taken  as  the  initial  measurement  position  for 
which  experimental  data  are  available. 

The  total  stationary  disturbance  amplitude  is  computed  from  the  mode-shape 
profiles  of  figure  5.62.  Three  different  measures  of  the  mode  shape  are  used  to  char- 
acterize the  disturbance  amplitude:  the  maximum,  average,  and  rms  of  the  mode  with 
respect,  to  Y . The  average  of  the  mode  corresponds  to  the  integral  of  | u'  | , and  the 
rms  measure  essentially  gives  the  integral  of  |?/|2.  Figure  5.63  shows  the  disturbance 
amplitude  distribution  for  [6|  12]  roughness  at  Rec  = 2.4  x 10°.  The  dashed  lines  mark 
the  absolute  size  of  the  mode  shapes  as  computed  by  each  of  the  three  measures.  The 
solid  lines  show  the  corresponding  amplification  factor  N.  The  initial  disturbance 
amplitude  at  x/c  = 0.10  is  used  as  the  reference  amplitude  for  each  /V-factor  curve 
[A0  in  equation  (3.16)].  The  ability  of  the  iV-factor  to  collapse  the  data  onto  a single 
curve  is  typical.  By  equation  (3.14),  the  roughly  constant  slope  of  the  /V-factor  curves 
indicates  that  the  disturbance  growth  is  approximately  exponential  for  x/c  < 0.25. 
For  xj  c > 0.30,  the  local  spatial  growth  rate  decreases  and  the  disturbance  ampli- 

2 The  straight-line  extrapolation  is  strictly  valid  only  for  a zero- pressure-gradient  flow  with  no 
disturbances.  However,  the  profile  curvature  at  the  wall  is  negligible  for  this  experiment  and  the 
technique  is  justified  a posteriori  by  figure  5.10. 
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tilde  saturates  due  to  nonlinear  effects.  This  nonlinear  amplitude  saturation  occurs 
well  before  the  transition  location  of  ( x/c)iT  = 0.52.  It  is  important  to  note  that  the 
quantitative  value  of  N is  not  relevant  because  of  the  arbitrariness  in  choosing  the 
normalizing  amplitude.  Only  the  slope  is  important  as  it  provides  the  local  spatial 
growth  rate. 

Figure  5.64  compares  the  experimental  iV-factor  (as  computed  from  the  maximum 
of  the  mode  shape  profiles)  with  various  theoretical  predictions  of  Haynes  (Haynes 
and  Reed  1996).  The  nonlinear  parabolized  stability  equations  (PSE)  results  are  com- 
puted using  initial  amplitudes  provided  by  the  experiment,  and  great  care  is  taken  to 
ensure  the  computations  are  presented  in  the  coordinate  system  to  which  the  mea- 
surements are  restricted  (see  section  3.4).  The  agreement  is  excellent,  especially  in 
predicting  the  amplitude  saturation.  (At  this  time,  the  Haynes  and  Reed  formulation 
does  not  contain  curvature,  which  is  known  to  be  stabilizing  and  may  account  for  the 
small  differences  in  the  disturbance  growth.)  In  contrast,  the  Orr-Sommerfeld  and 
linear  PSE  results  fail  to  predict  the  details  of  the  disturbance  growth.  The  early 
qualitative  agreement  with  the  linear  PSE  computations  indicates  that  the  nonlinear 
effects  are  initially  weak  up  to  x/c  = 0.25,  at  which  time  the  growth  rates  depart 
from  linear  behavior.  The  saturation  of  the  disturbance  amplitude  is  dramatic  for 
x/c  > 0.30.  In  this  region  the  (0,0)  mean-flow-distortion  mode  is  observed  in  the 
spanwise-averaged  boundary-layer  profile,  and  the  mode  shapes  exhibit  the  charac- 
teristic second  lobe  indicative  of  strong  nonlinear  effects.  Figure  5.64  absolutely  and 
conclusively  removes  all  doubt  about  the  importance  of  nonlinear  effects  for  the  sta- 
tionary crossflow  instability. 

5.5.4  Mode-Shape  Comparisons 

The  ability  of  the  nonlinear  PSE  to  capture  the  details  of  the  disturbance  growth 
can  be  assessed  by  comparing  the  measured  and  computed  mode  shapes.  These 
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comparisons  are  shown  in  figures  5.65  -5.72  for  0.10  < x/c  < 0.45.  The  theoretical 
mode  shapes  provided  by  Haynes  (Haynes  and  Reed  1996)  are  appropriately  projected 
onto  the  experimental  measurement  plane  as  outlined  in  section  3.4.  The  dots  simply 
differentiate  between  the  experimental  and  computed  results  and  do  not  indicate 
actual  measurement  points.  The  initial  conditions  for  the  nonlinear  PSE  are  chosen 
so  the  calculations  match  the  experiment  at  the  first  location  for  which  measurements 
are  available.  This  means  the  disturbance  amplitudes  are  equal  by  definition  at 
x/c  = 0.10  (figure  5.65).  The  mean  boundary-layer  distortions  are  initially  very 
weak,  which  explains  the  noise  in  the  experimental  data  at  10%  chord. 


The  nonlinear  PSE  result  overpredicts  the  maximum  amplitude  by  0.3%  at  x/c  = 
0.15  (figure  5.66),  and  their  is  evidence  of  some  stretching  in  the  Y direction.  These 
trends  continue  downstream,  and  explain  the  slight  difference  between  the  experi- 
mental and  nonlinear  PSE  iV- factor  curves  in  figure  5.64.  At  x/c  = 0.30,  both  the 
measured  and  computed  mode  shapes  begin  to  exhibit  the  second  lobe  indicating  the 
presence  of  strong  nonlinearities  (figure  5.69).  The  nonlinear  features  continue  to 
develop  until  the  final  measurement  location  at  x/c  = 0.45. 

It  is  important  to  focus  on  the  “shape”  of  the  mode  as  well  its  “size”  when  com- 
paring the  experimental  and  theoretical  data.  Although  the  computations  slightly 
overestimate  the  measured  data,  all  of  the  essential  features  of  the  stationary  dis- 
turbance are  captured  by  the  nonlinear  PSE.  In  particular,  the  development  of  the 
nonlinear  second  lobe  is  predicted  remarkably  well.  The  success  of  the  nonlinear 
PSE  and  its  agreement  with  the  experiment  underscores  the  significance  of  the  mean 
boundary-layer  distortion  as  the  dominant  aspect  of  the  stationary  crossflow  insta- 
bility. 
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5.5.5  Wavelength  Separation 

Spamvise  Scans  at  Constant  Y 

As  pointed  out  previously,  disturbance  amplitudes  computed  from  mode-shape 
profiles  contain  all  amplified  stationary  modes.  However,  the  mode  shapes  provide 
no  information  concerning  the  wavenumber  content  of  the  disturbance.  These  data  are 
important  in  light  of  the  sustained  region  of  nonlinear  interaction  shown  in  figure  5.64, 
hence  the  spamvise  scan  at  constant  Y technique  (section  4.4.2)  is  used  to  obtain 
individual-mode  growth  rates.  At  each  chord  position,  the  scan  is  performed  at  the 
height  corresponding  to  the  peak  of  the  total  disturbance  mode  shape.  This  allows  the 
amplitude  of  the  individual  modes  to  be  directly  compared  to  the  total  disturbance 
amplitude  computed  from  the  maximum  of  the  mode-shape  profile. 

Figures  5.73  5.81  show  the  hot-wire  surveys  for  0.05  < x/c  < 0.45.  The  data  are 
acquired  for  lie,  — 2.4  x 10(>  with  [6|12]  roughness.  These  plots  show  the  spamvise 
distribution  of  the  normalized  streamwise  velocity  at  a fixed  height  in  the  boundary 
layer.  The  velocity  trace  is  equivalent  to  that  which  would  be  obtained  by  taking 
a spamvise  slice  across  a set  of  wall-normal  boundary-layer  profiles  at  a constant 
height  above  the  airfoil  surface.  The  primary  difference  is  that  the  data  cover  the 
entire  spamvise  extent  of  the  available  measurement  region  (240  mm)  at  better  than 
four  times  the  spamvise  resolution  of  the  boundary- layer  scans.  This  allows  accurate 
decomposition  in  the  spectral  domain. 

The  development  of  the  stationary-wave-induced  distortion  of  the  mean  bound- 
ary layer  is  clearly  evident  for  x/c  > 0.15.  At  x/c  = 0.45  (figure  5.81),  the  velocity 
fluctuations  range  from  u/Ue  = 0.30  to  u/Ue  = 0.88,  and  the  rms  of  the  signal  (com- 
puted by  first  subtracting  the  mean  velocity  of  u/Ue  = 0.67)  is  18.5%.  This  value 
agrees  well  with  the  maximum  of  the  mode-shape  profile  (figure  5.37),  indicating  the 
scan  proceeded  at  the  desired  height  in  the  boundary  layer.  The  lack  of  any  overrid- 


ing  DC  trend  confirms  that  the  traverse-alignment  technique  sufficiently  models  the 
misalignment  between  the  airfoil  surface  and  the  hot-wire  traverse.  Moreover,  the 
uniformity  of  the  disturbance  over  the  entire  span  verifies  the  spanwise  independence 
of  the  instability,  and  once  again  validates  the  infinite-span  assumption. 

The  velocity  profiles  are  not  as  clean  near  the  leading  edge.  The  nominal  u/Ue  at 
x/c  — 0.10  drops  slightly  as  the  scan  progresses  (from  left  to  right  in  figure  5.74),  and 
the  data  appear  contaminated  by  background  noise.  These  features  are  exaggerated 
at  x/c  = 0.05  (figure  5.73).  This  is  not  caused  by  a failure  of  the  traverse-alignment 
technique,  but  instead  results  from  the  magnification  of  small  traverse  errors  by  the 
increasingly  large  du/dY  gradient  of  the  thin  boundary  layer.  (Recall  that  this  effect 
was  observed  in  the  total  disturbance  profiles  discussed  previously.)  Fortunately, 
the  disturbance  amplitude  is  not  a strong  function  of  V'  near  the  leading  edge,  as 
indicated  bv  the  relative  “flatness”  of  the  mode-shape  profile  for  small  x/c.  (shown 
in  figure  5.62).  Therefore,  the  “drift”  in  the  scan  can  be  tolerated  without  rendering 
the  data  useless. 

Spectral  Analysis 

The  spatial  power  spectra  for  the  spanwise  scans  are  shown  in  figures  5.82  5.90. 
These  plots  show  the  power  spectral  density  (PSD)  as  a function  of  spanwise  wave- 
length A*  (note  the  equivalent  label  A.,  in  the  figures).  The  power  estimates  are 
generated  with  the  FFT-bascd  power  spectrum  technique  described  in  section  4.4.2. 
With  64  times  spectral  smoothing  through  zero  padding,  the  method  accurately  ex- 
tracts the  peaky  spectra  with  sufficient  resolution  to  allow  disturbance  amplitude 
calculations  by  integration  of  the  PSD. 

The  data  for  0.20  < x/c  < 0.45  are  plotted  on  the  same  scale  to  highlight  the 
disturbance  growth.  The  spectra  are  dominated  by  a 12  mm  component  corresponding 
to  the  spacing  of  the  roughness  elements.  However,  a (dear  contribution  from  the  (0,  2) 
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mode  at  A2  = 6 mm  appears  for  xjc  > 0.35.  Careful  examination  of  these  data  reveals 
that  the  6 mm  component  first  becomes  measurable  at  xjc  — 0.30.  Moreover,  a very 
weak  4 mm  component  is  detected  for  xjc  > 0.35.  (These  features  are  too  small  to 
appear  on  the  scale  of  the  plots.)  This  wave  doubling  was  observed  by  Saric  and 
Yeates  (1985)  and  predicted  by  Reed  (1988).  Thus,  the  growth  of  the  (0,2)  and  (0,3) 
harmonics  and  their  nonlinear  interaction  with  the  fundamental  mode  leads  to  the 
amplitude  saturation  observed  in  figure  5.64. 

The  power  axis  is  magnified  to  show  the  small  disturbance  levels  for  xjc  < 0.15 
(figures  5.82  5.84).  The  12  mm  component  continues  to  dominate  the  PSD  at  xjc  ~ 

0.15,  however  a small  peak  appears  at  A2  = 7.2  mm  and  a tiny  bump  can  be  seen  at 
A2  = 6 mm.  Although  the  6 mm  feature  is  arguably  the  (0,  2)  mode,  it  is  rejected  since 
its  amplitude  is  dangerously  close  to  the  noise  floor  of  the  spectrum.  At  xjc  — 0.10 
the  fundamental  mode  at  Xz  — 12  mm  is  smaller  (but  still  measurable),  the  6 mm 
feature  has  disappeared,  and  the  7.2  mm  component  is  larger.  At  xjc  = 0.05  the 
disturbance  energy  is  almost  entirely  contained  in  the  7.2  mm  wavelength  and  its 
harmonic  at  Xz  — 3.6  mm.  In  particular,  there  is  no  evidence  of  the  fundamental 
crossflow  mode  (Xz  = 12  mm)  or  its  harmonics.  These  trends  are  caused  by  two 
separate  phenomena: 

1.  The  lack  of  any  measurable  crossflow  modes  for  small  xjc  is  simply  due  to 
the  infinitesimal  amplitude  of  the  stationary  disturbance  near  the  leading  edge. 
This  underscores  the  ability  of  the  k = 6 /mi  roughness  to  provide  uniform 
stationary  waves  without  excessive  initial  amplitudes. 

2.  The  emergence  of  the  7.2  mm  component,  as  well  as  the  increased  noise  in  the 
PSD,  is  caused  by  the  magnification  of  small  traverse  errors  in  the  thin  boundary 
layer  near  the  leading  edge.  This  has  been  discussed  previously,  and  results 
from  the  increase  in  the  dimensional  wall-normal  velocity  gradient  du/dY . The 
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7.2  mm  wavelength  itself  comes  from  a small,  once-per-revolution  eccentricity 
of  the  traverse  vertical  motion  (Z)  lead  screws  (Radeztsky  1994). 

It  should  be  emphasized  that  the  traverse  motion  is  highly  accurate  (as  shown  in 
table  2.2),  and  resolution  issues  are  a concern  only  under  the  influence  of  the  large 
velocity  gradient  of  a thin  boundary  layer.  Fortunately,  the  accuracy  of  the  spectral 
decomposition  technique  allows  the  crossflow  modes  to  be  separated  from  the  noise 
even  when  the  disturbance  amplitude  is  very  small. 

Individual- Mode  Disturbance  Amplitudes 

By  Parseval’s  theorem,  the  square  root  of  the  integral  of  the  PSD  equals  the  rms 
amplitude  of  the  signal.  Thus,  the  disturbance  energy  for  a single  mode  is  computed 
by  integrating  the  corresponding  peak  in  the  spectrum.'5  The  extent  of  the  peak  is 
defined  by  the  first  local  minimum  on  each  side. 

Figure  5.91  shows  the  disturbance  amplitude  as  a function  of  x/c  for  Rec  = 2.4  x 
106  with  [6 1 12]  roughness.  As  with  the  total  disturbance  amplitude  (figure  5.63), 
the  dashed  lines  represent  the  rms  amplitude  of  the  individual  modes,  while  the 
amplification  factor  N is  plotted  with  solid  lines.  The  iV-factors  are  relative  to  the 
point  at  which  the  mode  is  first  detected.  The  fundamental  (A2  = 12  mm)  mode 
contains  most  of  the  disturbance  energy,  however  the  (0,2)  and  (0,3)  modes  are  also 
amplified  in  the  region  of  nonlinear  amplitude  saturation. 

Recall  that  at  each  chord  location  the  spanwise  scan  is  taken  at  the  boundary-layer 
height  corresponding  to  the  maximum  of  the  mode-shape  profile.  This  permits  us  to 
compare  the  amplitude  of  the  individual  modes  with  that  of  the  total  disturbance 
computed  from  the  maximum  of  the  mode  shape.  Figure  5.92  shows  this  comparison. 
It  is  a remarkable  verification  of  the  two  vastly  different  measurement  techniques 

■'  Actually,  the  PSD  is  converted  from  mean-square  to  sum-square  power,  and  the  power  estimates 
in  each  peak  are  summed  to  provide  the  disturbance  energy.  This  avoids  introducing  error  with  an 
approximate  integration  technique. 
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that  the  amplitude  of  the  total  and  fundamental  (\z  = 12  mm)  disturbances  agree 
for  x j c ■ < 0.30.  In  particular,  this  validates  the  spectral  decomposition  technique  and 
highlights  the  ability  of  the  FFT — when  correctly  implemented  to  extract  peaky 
spectra  from  a limited  data  series.  Both  the  total  and  fundamental  disturbances 
exhibit  the  characteristic  nonlinear  amplitude  saturation  at  xfc  = 0.30.  At  this 
location,  the  amplitude  of  the  fundamental  diverges  from  the  total  disturbance  and 
the  (0,  2)  harmonic  at  Xz  = 6 mm  becomes  measurable.  The  6 mm  component 
contains  significant  energy  and  the  (0,3)  mode  (A-  = 4 mm)  becomes  unstable  in 
the  region  of  strong  nonlinear  interaction  ( x/c  > 0.35).  It  is  interesting  to  note  that 
the  amplitude  saturation  is  not  limited  to  the  fundamental  disturbance,  but  is  also 
observed  in  the  6 mm  mode.  Clearly,  the  amplitude  saturation  phenomenon  is  caused 
by  the  nonlinear  interaction  between  the  fundamental  disturbance  and  its  harmonics. 

5.6  Summary 

The  data  presented  in  this  chapter  provide  detailed  experimental  measurements  and 
theoretical  calculations  for  the  baseline  configuration  defined  as  Rec  = 2.4  x 106  with 
[6|  12]  roughness.  These  results  provide  important  insight  into  the  fundamental  nature 
of  the  stationary  crossflow  instability.  When  the  initial  conditions  are  governed  by 
“natural’'  surface  roughness,  the  stationary  waves  are  nonuniform  in  span  and  con- 
tain many  disturbance  modes.  Micron-sized  artificial  roughness  elements  placed  near 
the  leading  edge  generate  uniform  stationary  waves  with  a controlled  wavenumber 
content.  Detailed  hot-wire  measurements  show  early  departure  from  linear  behavior 
and  saturation  of  the  disturbance  amplitude.  Spectral  analysis  shows  that  this  satu- 
ration phenomenon  is  caused  by  the  nonlinear  growth  of  harmonic  disturbances  and 
their  interaction  with  the  fundamental  mode.  The  nonlinear  effects  are  observed  well 
before  transition.  PSE  computations  agree  remarkably  well  with  the  experiment  and 
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confirm  the  importance  of  nonlinear  effects. 

These  measurements  continue  to  highlight  the  sensitivity  to  roughness-induced 
initial  conditions  and  the  importance  of  the  receptivity  process  for  this  instability 
mechanism.  With  this  in  mind,  the  following  chapter  presents  a parametric  study  in 
which  the  disturbance  growth  is  investigated  under  varying  freestream  and  surface- 
roughness  conditions. 


CHAPTER  6 


Results — Part  2:  Parametric  Study 


This  chapter  discusses  the  experimental  results  for  data  sets  B through  T defined  in 
table  5.1,  which  is  repeated  on  the  following  page  for  convenience.  Along  with  data 
set  A , data  sets  B and  C examine  the  effects  of  varying  the  chord  Reynolds  number. 
Multiple-mode  crossflow  waves  are  produced  with  data  set  T>.  Finally,  data  sets  £ 
and  T combine  with  data  set  A to  investigate  the  dependence  on  initial  disturbance 
amplitude. 

6.1  Reynolds  Number  Variation 

Measurements  for  data  sets  B and  C are  presented  in  this  section.  The  roughness  con- 
figuration remains  at  k = 6 //,rn  with  12  mm  spanwise  spacing  (i.e.,  [6 1 1 2]  roughness) 
for  these  investigations.  Transition  occurs  somewhere  past  the  pressure  minimum  for 
data  set  B {Re,.  = 1.6  x 10fl),  and  a complete  set  of  boundary-layer  profile  scans 
are  taken  from  x/c  — 0.05  to  x/c  = 0.60  in  5%-chord  increments.  For  data  set  C 
(Re,.  = 3.2  x 106),  transition  occurs  at  x/c  — 0.32  with  some  turbulent  wedges  near 
x/c  = 0.30.  In  this  case,  boundary-layer  profile  scans  and  spanwise  scans  at  constant 
V'  are  obtained  for  0.05  < x/c  < 0.29.  The  last  measurement  location  is  chosen  to  en- 
sure that  no  data  are  acquired  in  a turbulent  boundary  layer.  The  roughness  Reynolds 
numbers  (Re/,)  are  0.061  and  0.17  for  Re,  = 1.6  x 106  and  3.2  x 106,  respectively. 
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Table  5.1:  Experimental  data  set  configuration. 


Data 

Set 

ReJ  106 

Roughness 

Scan  Type 

k [/mi 

Xz  [mm 

BL 

Span 

A 

2.4 

6 

12 

• 

• 

B 

1.6 

6 

12 

• 

C 

3.2 

6 

12 

• 

• 

V 

2.4 

6 

36 

• 

• 

£ 

2.4 

18 

12 

• 

• 

T 

2.4 

48 

12 

• 

• 

6.1.1  Effect  of  Decreasing  Rec 

The  streamwise  velocity  contours,  boundary-layer  profiles,  and  disturbance  profiles  for 
Rec  = 1.6  x 10°  are  shown  in  figures  6.1  6.12,  6.13-6.24,  and  6.25  6.36,  respectively. 
As  with  the  earlier  contour  plots,  the  10:1  V-to-2  aspect  ratio  is  chosen  to  enhance 
the  detail.  (From  this  point  forward,  the  10:1  aspect  ratio  will  be  assumed  unless 
otherwise  noted.)  The  disturbance  profiles  are  plotted  on  a consistent  scale  of  ±40% 
to  facilitate  direct  comparisons  with  data  set  A.  The  data  look  qualitatively  similar 
to  the  results  for  Re,  = 2.4  x 106,  except  the  disturbance  amplitudes  are  smaller  as 
expected.  This  is  a consequence  of  the  reduction  in  Re./,  and  growth  rate  with  the 
lower  chord  Reynolds  number. 

A careful  analysis  of  the  disturbance  motion  shows  that  stationary  crossflow  waves 
are  not  measurable  for  x/c  < 0.15.  This  is  illustrated  by  comparing  the  disturbance 
velocity  contours  at  x/c  — 0.15  and  x/c  = 0.20  (figures  6.37  and  6.38).  The  periodicity 
of  the  velocity  surplus  and  deficit  clearly  indicates  a stationary  wave  at  x/c  = 0.20. 
However,  at  x/c  = 0.15  no  definite  structure  stands  out  above  the  background  noise 
of  the  measurements.  Thus,  the  data  for  x/c  < 0.15  are  not  considered. 

The  mode-shape  profiles  for  0.20  < x/c  < 0.60  are  shown  in  figure  6.39.  It  is 
worth  emphasizing  that  the  symbols  simply  identify  the  curves  and  do  not  represent 
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measurement  points.  The  strong  distortion  and  second  lobe  that  dominate  the  Rec  = 
2.4  x 106  data  are  missing,  although  the  mode  shapes  at  x/c  = 0.55  and  x/c  — 0.60 
are  beginning  to  “push  out”  near  Y = 2.3  mm  in  an  early  sign  of  nonlinear  effects. 
Tracking  the  maximum,  average,  and  rms  of  the  mode  shapes  yields  the  amplitude 
distribution  shown  in  figure  6.40.  As  noted  earlier,  the  amplification  factor  N collapses 
the  three  measures  to  a single  curve.  It  is  interesting  that  a slight  amplitude  saturation 
is  observed  for  x/c  > 0.50  even  though  the  disturbance  amplitude  is  relatively  small. 
This  reaffirms  the  importance  of  nonlinear  effects  for  boundary  layers  dominated  by 
the  stationary  crossflow  instability. 

6.1.2  Effect  of  Increasing  Rec 

Increasing  the  chord  Reynolds  number  to  3.2  x 10f>  provides  a third  data  set.  with 
[6 1 1 2]  roughness.  Figures  6.41-6.46,  6.47-6.52,  and  6.53-6.58  show  the  streamwise 
velocity  contours,  boundary-layer  profiles,  and  disturbance  profiles,  respectively.  The 
data  look  very  similar  to  those  obtained  for  Rec  = 2.4  x 106.  The  distortion  of  the 
mean  boundary  layer,  however,  develops  much  more  quickly  due  to  the  increased 
growth  rates. 

The  disturbance  remains  too  weak  to  measure  at  x/c  = 0.05  (figures  6.41,  6.47, 
and  6.53)  even  though  the  roughness  Reynolds  number  (Rek)  is  50%  larger  under 
these  conditions  than  for  Rec  = 2.4  x 10°.  (This  will  be  confirmed  by  the  spanwise 
scans  presented  below.)  Thus,  the  disturbance  amplitudes  are  computed  for  x/c  > 
0.10.  Figure  6.59  shows  the  mode-shape  profiles  for  these  chord  locations.  The  mode 
shape  first  begins  to  distort  at  x/c  - 0.20,  and  along  with  the  asymmetry  of  the 
disturbance  profiles  at  this  location  (figure  6.56),  indicates  that  nonlinearities  are 
present.  From  x/c  = 0.20  to  x/c  = 0.29  the  mode  shape  quickly  develops  the  second 
lobe  characteristic  of  strong  nonlinear  effects.  Figure  6.60  shows  the  total  disturbance 
amplitude  and  amplification  factor  N computed  from  the  mode  shapes.  Nonlinear 
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behavior  is  observed  for  x/c  > 0.20,  with  strong  amplitude  saturation  occurring  for 
x/c  > 0.25. 

Data  set  C also  contains  a full  set  of  spanwise  scans  suitable  for  spectral  analysis. 
Figures  6.61  6.66  show  the  velocity  profiles  for  0.05  < x/c  < 0.29.  The  corresponding 
spectra  are  presented  in  figures  6.67-6.72.  The  scans  are  contaminated  by  noise  for 
small  x/c  due  to  the  magnification  of  small  traverse  errors  by  the  large  wall-normal 
gradient  du/dY.  This  effect  was  observed  for  Rec  = 2.4  x 106,  but  is  intensified 
for  Rt  = 3.2  x 10f)  because  of  the  increase  in  freestream  velocity  and  decrease  in 
boundary-layer  thickness.  Consequently,  the  hot-wire  scan  at  x/c  = 0.05  is  over- 
whelmed by  the  A2  = 7.2  mm  traverse  eccentricity  discussed  in  section  5.5.5,  and  no 
crossflow  modes  can  be  extracted  (figure  6.67). 

The  background  noise  level  is  smaller  at  x/c  — 0.10,  and  figure  6.68  shows  that 
there  is  measurable  energy  in  the  fundamental  mode  (A,  = 12  mm).  The  (0,2) 
harmonic  at  A;  = 6 mm  also  appears  to  be  amplified,  but  we  are  forced  to  ignore 
it  since  it  is  arguably  the  same  level  as  the  background  noise.  By  x/c  = 0.15,  the 
fundamental  crossflow  wavelength  dominates  and  the  (0,  2)  mode  contains  significant 
energy  (figure  6.69).  In  fact,  the  A,  = 6 mm  component  can  be  clearly  seen  in 
the  velocity  profile  (figure  6.63).  This  trend  continues  downstream  until  the  last 
measurement  location  at  x/c  = 0.29.  The  (0,3)  mode  (A;  = 4 min)  becomes  unstable 
at  x/c  = 0.25  and  undergoes  strong  growth  until  x/c  = 0.29.  Unlike  the  Rec  = 
2.4  x 10<J  case,  the  data  at  x/c  = 0.29  show  evidence  of  a weak  (0,  4)  mode  = 3 mm) 
just  before  transition.  Its  amplitude,  however,  is  an  order  of  magnitude  below  that 
of  the  (0.  3)  mode  and  is  barely  detectable  above  the  background  noise  level  of  the 
spectrum. 


These  features  are  encapsulated  in  figure  6.73,  which  shows  the  disturbance  ampli- 
tude and  corresponding  amplification  factor  N for  the  dominant  stationary  crossflow 


modes.  These  data  are  computed  by  integrating  the  appropriate  peaks  in  the  spectra. 
As  with  the  results  for  Rec  = 2.4  x 10b  (figure  5.91),  the  amplitude  of  the  fundamen- 
tal mode  saturates  in  the  region  where  the  harmonics  contain  significant  disturbance 
energy.  However,  in  contrast  to  figure  5.91  the  amplitude  of  the  fundamental  actually 
decreases  at  the  last  measurement  location  and  the  (0,  2)  and  (0,  3)  modes  show  no 
evidence  of  amplitude  saturation. 

Figure  6.74  compares  the  amplitude  of  the  individual  modes  to  that  of  the  total 
disturbance  computed  from  the  maximum  of  the  mode-shape  profiles.  (Recall  that 
the  spanwise  scans  proceed  at  the  height  corresponding  to  the  maximum  of  the  mode 
shape.)  This  plot  displays  features  similar  to  the  results  for  Rec  = 2.4  x 106  (shown 
in  figure  5.92).  In  particular,  the  nonlinear  amplitude  saturation  is  characterized  bv 
the  decrease  in  the  amplitude  of  the  fundamental  mode  (with  respect  to  the  total 
disturbance)  and  the  corresponding  growth  of  the  harmonics.  The  accuracy  of  the 
two  measurement  techniques  is  indicated  by  the  excellent  agreement,  between  the  total 
and  fundamental  disturbances  for  x/c  < 0.20. 

6.1.3  Disturbance  Amplitude  Comparison 

The  data  presented  above  are  combined  with  the  results  for  Rec  = 2.4  x 10(‘  (presented 
in  section  5.5)  to  investigate  the  effect  of  chord  Reynolds  number  on  the  stationary 
crossflow  instability. 

Figure  6.75  shows  the  total  disturbance  amplitude  distribution  for  Rec  = 1.6  x 106, 
2.4  x 10fi,  and  3.2  x 10(i  with  [6|  12]  roughness.  These  curves  are  generated  from 
the  maximum  of  the  mode-shape  profiles.  As  before,  the  dashed  lines  show  the 
absolute  disturbance  amplitude  and  the  solid  lines  represent  the  amplification  factor 
N.  Amplitude  saturation  is  evident  for  all  Reynolds  numbers.  Recall  that  transition 
is  occurring  at  x/c  = 0.32  for  Rev  — 3.2  x 106,  at  x/c  = 0.52  for  Rec  = 2.4  x 10(>,  and 
somewhere  past  the  pressure  minimum  for  R.ec  = 1.6  x 10fi.  It  is  interesting  to  note 
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that  the  amplitude  of  the  stationary  crossflow  disturbance  is  approximately  equal  at 
transition  for  the  two  largest  Reynolds  number  cases.  The  jV- factors  are  different, 
but  this  is  due  to  the  difference  in  the  reference  amplitude  at  x/c  = 0.10.  It  appears 
there  may  exist  an  equilibrium  amplitude  toward  which  the  disturbance  grows.  This 
idea  will  be  revisited  later  in  this  chapter. 

6.2  Multiple-Mode  Crossflow  Waves 

Flow  visualizations  presented  in  section  5.2  indicate  that  increasing  the  spanwise 
spacing  of  the  roughness  elements  generates  stationary  crossflow  waves  containing 
multiple  disturbance  modes.  Moreover,  the  characteristic  saw-tooth  pattern  of  the 
transition  front  becomes  more  regular  and  the  average  transition  location  moves  for- 
ward when  the  roughness  spacing  is  increased.  Clearly,  changing  the  initial  conditions 
bv  altering  the  surface  roughness  distribution  affects  the  stability  characteristics  and 
warrants  further  investigation. 

This  section  describes  the  measurements  for  data  set  V in  table  5.1.  The  roughness 
height  remains  at  G /an,  but  the  spanwise  spacing  is  increased  to  36  mm.  This 
is  denoted  as  [6 1 36]  roughness  in  the  shorthand  notation  described  in  section  5.1. 
The  chord  Reynolds  number  is  2.4  x 106.  The  disturbance  structure  is  mapped  for 
0.05  <:C  x/c  ^ 0.45  using  the  boundary-layer  profile  technique.  Spanwise  scans  at 
constant  Y provide  individual-mode  growth  rates  for  0.10  < x/c  < 0.45. 

6.2.1  Disturbance  Structure 

The  streamwise  velocity  contours  for  0.05  < x/c  < 0.45  are  plotted  in  figures  6.76 
6.84.  Figures  6.85-6.93  show  the  boundary-layer  profiles  from  which  the  contours  are 
generated.  The  corresponding  disturbance  profiles  are  presented  in  figures  6.94-6.102. 
As  with  the  previous  data  sets,  no  stationary-wave-induced  distortions  of  the  mean 
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boundary  layer  are  measurable  at  xjc  = 0.05  (figure  G.76).  In  addition,  figure  6.77 
shows  no  measurable  stationary  disturbance  at  xjc.  — 0.10.  This  is  in  contrast  to 
the  data  for  [6|12]  roughness  (section  5.5),  which  contain  measurable  energy  in  the 
fundamental  crossflow  mode  at  10%  chord. 

The  data  at  x/c  = 0.15  (figures  6.78,  6.87,  and  6.96)  show  the  first  indication  of  un- 
stable crossflow  modes.  Figure  6.79  clearly  exhibits  the  distortion  of  the  mean  bound- 
ary layer  at  xjc  = 0.20,  and  there  is  an  apparent  36  mm  feature  corresponding  to  the 
roughness  spacing.  The  disturbance  profiles  become  asymmetric  at  xjc  = 0.25  due 
to  the  nonlinear  interaction  among  the  various  modes  (figure  6.98).  By  xjc  = 0.30, 
the  velocity  contour  (figure  6.81)  displays  a more  complicated  structure  indicating 
the  presence  of  multiple  stationary  modes.  This  pattern  develops  downstream  into 
the  quite  dramatic  distortion  of  the  mean  boundary  layer  at  xjc  = 0.45  (figure  6.84). 
The  primary  features  are  36  mm  apart  and  represent  the  fundamental  crossflow  mode. 
Unlike  the  12  mm  forcing,  however,  there  is  clear  evidence  of  significant  energy  in 
the  harmonic  wavelengths.  These  features  are  accentuated  in  the  contour  plot  of  the 
velocity  surplus/deficit  shown  in  figure  6.103.  The  ability  of  the  artificial  roughness 
to  provide  uniform  initial  conditions  even  for  large  spanwise  spacing  is  indicated  by 
the  periodicity  of  the  boundary-layer  distortions. 

The  mode-shape  profiles  for  0.15  < xjc  < 0.45  are  plotted  in  figure  6.104.  Once 
again,  the  symbols  are  used  to  identify  the  curves  but  do  not  indicate  measurement 
points.  The  nonlinear  distortion  appears  at  xjc  = 0.35,  after  which  the  mode  shape 
quickly  develops  the  now-familiar  second  lobe  indicating  the  presence  of  st  rong  nonlin- 
ear effects.  As  a measure  of  the  total  disturbance  energy,  the  profiles  look  qualitatively 
similar  to  those  obtained  with  [6|  12]  roughness.  However,  one  should  be  cautioned 
against  placing  too  much  emphasis  on  this  comparison.  The  velocity  contours  clearly 
show  the  dramatic  difference  in  the  disturbance  field  for  the  two  roughness  spacings. 
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Figure  6.105  shows  the  chordwise  distribution  of  the  total  disturbance  amplitude 
and  amplification  factor  N.  These  curves  are  computed  from  the  maximum,  average, 
and  rms  of  the  mode-shape  profiles.  Amplitude  saturation  occurs  for  x/c  > 0.35. 

6.2.2  Wavelength  Separation 

The  rich  modal  content  of  the  disturbance  structure  for  [6|36]  roughness  is  illuminated 
by  spanwdse  scans  and  spectral  analysis.  The  scans  are  performed  at  the  boundary- 
layer  height  corresponding  to  the  maximum  of  the  mode  shapes  for  0.10  < x/c  < 0.45. 
The  hot-wire  surveys  are  presented  in  figures  6.106-6.113.  Figures  6.114  6.122  show 
the  power  spectral  density  for  these  cases.  Unstable  crossflow  modes  are  first  detected 
at  x/c  — 0.15.  The  spectrum  at  this  location  (figure  6.115)  shows  energy  in  the  (0,  3). 
(0,4),  and  (0,5)  modes  (Az  = 12  mm,  9 mm,  and  7.2  mm).  The  (0,5)  mode  must 
be  ignored  since  7.2  mm  happens  (quite  coincidentally)  to  be  the  exact  wavelength 
of  the  traverse-related  disturbance  discussed  in  section  5.5.5.  A small  amount  of 
energy  appears  in  the  (0,2)  mode  (A*  = 18  mm),  however  this  is  disregarded  as  it  is 
essentially  within  the  noise  of  the  spectrum.  It  is  interesting  that  the  fundamental 
mode  (A.  = 36  mm)  is  too  weak  to  measure  at  x/c  = 0.15  even  though  several 
harmonics  are  amplified. 

The  velocity  profile  at  x/c  = 0.20  (figure  6.108)  shows  the  first  sign  of  a 36  mm 
feature  corresponding  to  the  roughness  spacing.  However,  the  spectrum  (shown  in 
figure  6.116)  indicates  the  energy  in  the  fundamental  mode  is  infinitesimal  compared 
to  the  (0,2) -(0,5)  modes.  In  fact,  based  on  a consistent  definition  of  the  spectral 
noise,  the  36  mm  “peak”  must  be  ignored  even  though  its  effect  is  evident  in  the 
velocity  profile.  It  appears  that  the  fundamental  mode,  although  tiny  in  amplitude, 
is  providing  enough  energy  to  “pump”  the  more  unstable  wavelengths  on  a 36  mm 
scale.  The  fundamental  mode  itself  becomes  measurable  at  x/c  — 0.25  (figure  6.117), 
although  the  harmonics — especially  at  Az  = 12  mm  and  9 mm  contain  most  of  the 
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disturbance  energy. 

The  higher  harmonics  become  unstable  at  x/c  = 0.30.  The  spectrum  at  this 
location,  plotted  in  figure  6.118,  shows  the  emergence  of  the  (0,  6)-(0,  8)  modes 
(A*  = 6 mm,  5.1  mm,  and  4.5  mm).  These  modes  continue  to  grow  downstream 
until  the  last  measurement  position  at  x/c  = 0.45.  At  this  location  the  velocity 
profile  (figure  6.113)  clearly  shows  the  effect  of  multiple  unstable  modes,  which  are 
well-defined  in  the  spectrum  (figure  6.121).  It  is  important  to  note  that  the  broad- 
ening of  the  PSD  with  increasing  A is  an  artifact  of  plotting  the  spectrum  in  the 
wavelength  versus  the  wavenumber  domain.  Focusing  on  the  small-wavelength  dis- 
turbances (figure  6.122)  shows  unstable  harmonics  up  to  and  including  the  (0, 9)  mode 
(A*  = 4 mm).  There  appear  to  be  no  amplified  subharmonics  (wavelength  doubling) 
of  the  roughness  spacing. 

Integrating  the  spectral  power  in  each  peak  provides  the  disturbance  amplitude  for 
each  unstable  mode.  Figure  6.123  shows  the  amplification  factor  N for  all  wavelengths 
amplified  by  the  36  mm  forcing.  The  (0, 5)  mode  (A.  = 7.2  mm)  is  not  plotted  because 
of  its  contamination  by  the  vertical  traverse  motion.  The  fundamental  disturbance 
is  not  detected  until  x/c.  — 0.25,  but  has  the  largest  growth  rate.  The  harmonics  at 
A2  = 18  mm,  12  mm,  and  9 mm  are  detected  before  the  fundamental,  yet  display 
smaller  growth  rates.  Amplitude  saturation  occurs  for  the  Az  = 12  mm,  9 mm, 
6 mm,  and  5.1  mm  disturbances.  Although  the  amplitude  of  the  short  wavelength 
disturbances  (Az  = 4.5  mm  and  4 mm)  is  very  small,  they  are  growing  as  we  move 
toward  transition. 

Since  the  individual  modes  are  detected  at  different  locations  with  different  initial 
amplitudes,  the  reference  amplitude  for  the  iV-factor  calculations  is  different  for  each 
wavelength.  Consequently,  direct  comparisons  between  the  value  of  N are  not  mean- 
ingful. In  this  case  it  is  instructive  to  consider  the  disturbance  amplitude  distribution 


102 


shown  in  figure  G.124.  This  plot  shows  the  rms  amplitude  of  the  amplified  distur- 
bances as  a function  of  x/c,  and  can  be  used  to  compare  the  relative  disturbance 
energy  in  individual  modes. 

Based  on  preliminary  nonlinear  PSE  calculations,  it  appears  the  36  mm  forcing  is 
providing  initial  energy  for  several  harmonics  as  well  as  the  fundamental  disturbance. 
This  explains  the  early  growth  of  the  (0, 3)  and  (0,  4)  modes  (Az  = 12  mm  and  9 mm), 
which  are  near  the  most  unstable  wavelength  for  this  experiment.  The  importance  of 
the  spectral  content  of  the  roughness  distribution  may  also  indicate  that  the  roughness 
excite  the  stationary  instability  through  a wave-like  pressure  disturbance  and  not  a 
local  injection  of  vorticity.  However,  one  should  be  cautioned  that  a more  detailed 
receptivity  study  is  required  to  fully  address  this  topic. 


6.3  Initial  Amplitude  Variation 


The  height  of  the  artificial  roughness  elements  has  been  fixed  at  6 pm  in  all  the 
results  presented  up  to  this  point.  In  this  section,  the  effect  of  increasing  the  initial 
disturbance  amplitude  by  changing  the  roughness  height  is  investigated. 

Data  sets  £ and  T defined  in  table  5.1  contain  measurements  for  two  additional 
roughness  heights:  18  pm  and  48  pm.  The  details  of  the  elements  are  described  in 
section  3.3.3.  In  both  cases,  the  spanwise  spacing  of  the  elements  is  12  mm  and  the 
data  are  obtained  for  Rec  - 2.4  x 10°.  Under  these  conditions,  the  roughness  Reynolds 
numbers  are  1.0  and  7.0  for  k — 18  pm  and  48  pm,  respectively.  The  measurements 
presented  in  this  section  are  combined  with  the  results  for  [6|  12]  roughness  at  Rec  = 
2.4  x l()fi  (data  set  A)  to  provide  disturbance  amplitudes  covering  a nearly  two-order- 
of-magnitude  change  in  Re*,. 
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6.3.1  Data  Set  S 

Figures  6.125-6.133  show  the  streamwise  velocity  contours  for  [18|  12]  roughness  at 
Rec  — 2.4  x 106.  The  mean  boundary-layer  profiles  used  to  construct  the  contour  plots 
are  presented  in  figures  6.134-6.142,  and  figures  6.143  6.151  show  the  corresponding 
disturbance  profiles.  In  contrast  to  the  [6|12]  roughness  at  the  same  Reynolds  num- 
ber, measurable  stationary  crossflow  waves  are  detected  at  x/c  = 0.05  (this  will  be 
confirmed  by  spectral  analysis).  The  distortion  of  the  mean  boundary  layer  is  clearly 
evident  in  the  velocity  contour  at  x/c  = 0.10  (figure  6.126).  The  asymmetry  of  the 
disturbance  profiles  at  x/c  = 0.20,  shown  in  figure  6.146,  indicates  the  early  “rollover” 
of  the  streamwise  velocity  component  due  to  the  distortions  of  the  mean  boundary 
layer. 

The  stationary  wave  mode-shape  profiles  are  plotted  in  figure  6.152  for  0.05  < 
x/c  < 0.45.  The  nonlinear  distortion  leading  to  the  second  lobe  occurs  by  x/c  = 0.25 
and  continues  to  develop  downstream.  Figure  6.153  shows  the  total  disturbance 
amplitude  computed  from  the  maximum,  average,  and  rms  of  the  mode  shapes.  Al- 
though qualitatively  similar  to  the  [6 1 12]  roughness  case  (figure  5.63),  there  are  several 
important  differences: 

1.  The  initial  disturbance  amplitude  is  larger  and  can  be  measured  at  x/c  — 0.05. 

2.  The  amplitude  saturation  is  more  pronounced  and  occurs  earlier,  at  x/c  ss  0.25 
(saturation  occurs  at  x/c  ~ 0.3  for  the  k — 6 pm  roughness). 

3.  The  saturation  amplitude  at  x/c  = 0.45  is  comparable  to  the  6 pm  roughness 
case.  (This  will  be  discussed  in  detail  later.) 

4.  The  disturbance  amplitude  shows  a “dip”  at  x/c  = 0.35  which  is  not  present 
for  the  [6|  12]  roughness.  In  fact,  the  amplitude  measures  show  a “two-stage” 


104 

saturation  process  described  by  an  initial  amplitude  saturation  at  x/c  — 0.35, 

followed  by  second  growth  to  x/c  = 0.40  and  a final  saturation  at  x/c  — 0.45. 

Spanwise  scans  and  spectral  decomposition  provide  important  information  re- 
garding the  two-stage  amplitude  saturation.  The  velocity  profiles  and  corresponding 
powder  spectra  are  shown  in  figures  6.154-6.162  and  6.163-6.171.  The  scans  are  per- 
formed at  the  boundary-layer  height  corresponding  to  the  maximum  of  the  mode 
shapes.  The  spectrum  at  x/c  = 0.05  (figure  6.163)  displays  the  familiar  Xz  = 7.2  mm 
traverse-induced  noise,  however  the  fundamental  crossflow  mode  (Az  = 12  mm)  is 
clearly  amplified  at  this  location.  Figure  6.164  indicates  that  the  (0,2)  harmonic 
(Az  = 6 mm)  becomes  unstable  by  x/c  = 0.10.  The  (0,3)  mode  (\z  = 4 mm)  is 
detected  at  x/c  — 0.25. 

The  disturbance  amplitude  and  amplification  factor  N for  the  unstable  crossflow 
modes  are  plotted  in  figure  6.172.  Figure  6.173  shows  the  same  data  combined  with 
the  total  disturbance  amplitude  computed  from  the  maximum  of  the  mode-shape 
profiles.  The  fundamental  mode  exhibits  the  same  two-stage  amplitude  saturation 
that  characterizes  the  total  disturbance.  Moreover,  this  behavior  is  also  observed  in 
the  harmonic  wavelengths.  The  Az  = 6 mm  mode  show's  early  growth,  saturates  at 
x/c  = 0.25,  then  grows  again  only  to  saturate  a second  time  at  x/c  — 0.45.  Even  the 
(0,  3)  mode  (A,  = 4 mm)  shows  the  same  general  trend.  It  is  no  surprise  that  the 
strong  nonlinear  effects  beginning  at  x/c  = 0.25  correlate  with  the  resurgence  of  the 
6 mm  component  and  emergence  of  the  4 mm  wavelength. 

6.3.2  Data  Set  JF 

The  measurements  presented  above  are  repeated  for  [48|  12]  roughness.  The  stream- 
wise  velocity  contours,  mean  boundary-layer  profiles,  and  disturbance  profiles  for 
0.05  < x/c  < 0.45  are  showm  in  figures  6.174-6.182,  6.183-6.191,  and  6.192  6.200, 
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respectively.  The  total  disturbance  mode-shape  profiles  are  plotted  in  figure  6.201. 
As  with  the  corresponding  plots  for  the  previous  data  sets,  the  symbols  do  not  in- 
dicate measurement  points,  but  are  simply  used  as  curve  identifiers.  The  data  are 
very  similar  to  those  obtained  for  [18|  12]  roughness.  However,  there  are  two  notable 
differences.  First,  the  distortion  of  the  mode  shape  and  development  of  the  second 
lobe  begins  at  x/c  = 0.20,  compared  to  xjc  = 0.25  for  the  k = 18  /mi  roughness. 
This  indicates  that  nonlinear  effects  develop  more  quickly  with  increased  initial  dis- 
turbance amplitude,  as  expected.  Second,  the  disturbance-layer  thickness — defined 
as  the  height  at  which  the  mean  boundary-layer  distortions  vanish — has  increased  by 
approximately  0.5  mm  at  x/c  = 0.45  (about  a 1 mm  increase  with  respect  to  the 
[6 1 12]  roughness). 

The  total  disturbance  amplitude  and  amplification  factor  N determined  from  the 
mode-shape  profiles  are  plotted  in  figure  6.202.  The  trends  noted  for  the  k = 18  /mi 
roughness  continue  to  develop.  The  initial  amplitude  is  larger  and  saturation  occurs 
earlier,  now  at  x/c  « 0.20.  The  two-stage  saturation  phenomenon  still  appears  at 
x/c  = 0.35,  and  the  saturation  amplitude  remains  relatively  unchanged  from  both 
the  [6|  12]  and  [18|  12]  roughness  cases.  In  contrast  to  these  previous  measurements 
(figures  5.63  and  6.153),  the  disturbance  amplitude  actually  decreases  from  x/c  = 0.40 
to  x/c  = 0.45. 

Spanwise  scans  performed  at  the  boundary-layer  height  corresponding  to  the  maxi- 
mum of  the  mode-shape  profiles  give  the  disturbance  amplitude  in  individual  crossflow 
modes.  The  velocity  profiles  are  shown  in  figures  6.203-6.211.  The  corresponding 
spectra  are  plotted  in  figures  6.212  6.220.  Integrating  the  spectral  peaks  provides 
the  individual-mode  rms  amplitude  and  jV-factor  curves  shown  in  figure  6.221.  The 
fundamental  crossflow  mode  exhibits  the  two-stage  saturation  phenomenon,  however 
the  (0, 2)  harmonic  (A*  - 6 mm)  no  longer  displays  this  feature.  Moreover,  the 
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6 mm- wavelength  disturbance  is  not  detected  until  x/c  = 0.15,  whereas  it  contains 
measurable  energy  at  x/c  = 0.10  for  [18)12]  roughness.  On  the  other  hand,  the  (0,3) 
mode  (A,  = 4 mm)  becomes  unstable  earlier  but  otherwise  behaves  similarly  to  that 
for  A:  = 18  /mi  roughness.  Comparing  these  data  to  the  total  disturbance  amplitude 
(figure  6.222)  shows  the  typical  trends.  The  region  of  nonlinear  amplitude  saturation 
is  characterized  by  a decrease  in  the  amplitude  of  the  fundamental  mode  (with  respect 
to  the  total  disturbance)  and  the  corresponding  growth  of  the  harmonic  disturbances. 

6.3.3  Disturbance  Amplitude  Comparison 

The  total  disturbance  mode-shape  profiles  offer  the  opportunity  to  quantify  the  effects 
of  roughness  height  on  the  growth  and  structure  of  the  stationary  waves.  Figure  6.223 
shows  the  mode  shapes  at  x/c  = 0.10  and  x/c  = 0.45  for  [6)12],  [18|  12],  and  [48)12] 
roughness.  The  larger  roughness  generates  a much  larger  initial  amplitude,  however 
the  disturbance  profiles  relax  downstream  to  a similar  mode  shape.  The  disturbance 
energy  is  redistributed  higher  in  the  boundary  layer  with  increasing  initial  amplitude, 
actually  increasing  the  disturbance-layer  thickness  as  noted  earlier. 

Figure  6.224  shows  the  total  disturbance  amplitude  distributions  for  the  three 
roughness  heights.  These  curves  are  simply  the  maximum  measures  of  the  mode- 
shapes  extracted  from  figures  5.63,  6.153,  and  6.202.  As  indicated  above,  the  initial 
amplitude  and  growth  rate  increase  with  larger  roughness.  However,  the  amplitude 
prior  to  transition  remains  essentially  unchanged.  It  appears  that  once  the  nonlinear 
('fleets  are  strong  enough  to  cause  saturation,  the  saturation  amplitude  is  fixed. 

Figure  6.225  shows  the  amplification  factor  N for  the  amplitude  distributions  of 
figure  6.224.  The  iV-factors  are  computed  relative  to  x/c  = 0.10  since  this  is  the  first 
location  at  which  measurements  are  available  for  all  roughness  heights.  Because  of 
the  larger  initial  amplitude  and  growth  rate,  the  disturbance  amplitude  reaches  the 
saturation  level  earlier  for  the  larger  roughness.  Consequently,  the  relative  growth 


107 


decreases  with  increased  roughness  height  and  the  amplitude  saturation  becomes  more 
pronounced.  When  analyzing  figure  6.225,  it  is  important  to  bear  in  mind  that  the 
numerical  value  of  N is  irrelevant  because  the  initial  conditions  are  different  for  each 
roughness  configuration.  Only  the  slope  is  important  as  it  provides  the  local  spatial 
growth  rate. 

Recall  that  transition  occurs  at  xjc  k,  0.50  for  all  roughness  heights  considered 
(section  5.2).  In  light  of  this,  figure  6.225  shows  that  for  distributed  surface  rough- 
ness, nonlinear  effects  can  dominate  the  stationary  crossflow  instability  well  before 
transition.  In  the  case  of  [48|  12]  roughness,  for  example,  strong  nonlinear  interactions 
and  amplitude  saturation  are  observed  over  60%  of  the  extent  of  the  laminar  bound- 
ary layer.  The  insensitivity  of  the  transition  location  to  the  height  of  distributed 
roughness  is  in  contrast  to  Radeztsky  et  al.  (1993a),  who  showed  that  a single  rough- 
ness element  induces  an  isolated  transition  w'edge  whose  location  is  a strong  function 
of  the  roughness  height.  Moreover,  when  compared  with  Radeztsky  et  al.  (1993a), 
the  transition  location  in  the  case  of  distributed  roughness  is  at  a larger  xfc  than 
for  isolated  roughness.  Clearly  there  are  many  receptivity  issues  that,  warrant  further 
investigation. 


6.4  Special  Considerations 

6.4.1  Which  Wavelength  Is  Conserved? 

In  section  3.5.2  it  was  noted  that  linear,  parallel  stability  analyses  require  certain 
constraints  on  one  or  more  of  the  parameters  of  equation  (3.13)  in  order  to  obtain 
solutions  to  the  Orr-Soinmerfeld  equation.  Many  so-called  integration  strategies  have 
been  proposed,  each  attempting  to  choose  appropriate  constraints  based  on  some 
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assumed  physical  behavior  of  the  instability.1  Arnal  (1994)  reviews  the  common 
st  rategies,  which  typically  involve  an  assumption  concerning  the  direction  of  growth, 
the  direction  of  propagation,  and/or  the  wavelength  of  the  unstable  disturbance. 

No  single  integration  strategy  is  universally  accepted  as  providing  the  most  ap- 
plicable results  for  the  general  three-dimensional  flowfield.  However,  it  is  possible 
to  choose  physically  meaningful  constraints  for  infinite-span  swept  wings.  In  this 
case  the  disturbance  field  must  be  the  same  at  each  spanwise  location,  from  which 
it  follows  that  the  spanwise  growth  rate  //  = 0 (Mack  1988).  This  means  the  dis- 
turbance grows  only  in  the  chordwise  direction,  hence  the  need  to  impose  an  ad  hoc 
constraint  on  the  growth  direction  is  removed.  The  experimental  evidence  support- 
ing lii  = 0 is  provided  by  any  of  the  streamwise  velocity  contours  and  spanwise  scans 
presented  earlier  in  this  chapter.  For  example,  figures  5.29  and  5.81  clearly  show  a 
spanwise-invariant  disturbance  field. 

The  appropriate  constraint  concerning  the  crossflow  wavelength  (or  equivalently 
the  direction  of  propagation  of  the  disturbance)  is  not  so  obvious.  In  early  exper- 
iments at  ASU,  Dagenhart  (1992)  measured  the  stationary  disturbance  wavelength 
by  counting  the  number  of  light  and  dark  streaks  that  appear  in  a naphthalene  flow 
visualization.  To  within  the  resolution  of  this  technique,  the  crossflow  wavelength 
normal  to  the  vortex  axis  (A(:f)  appeared  to  remain  constant  over  the  model.  In  addi- 
tion, no  “drop-outs”  or  other  adjustments  to  the  vortex  spacing  were  observed.  Thus. 
Dagenhart  concluded  that  the  appropriate  wavelength  constraint  for  the  linear,  par- 
allel stability  calculations  on  infinite-span  wings  is  constant  Arf.  These  early  results, 
however,  were  restricted  by  certain  fundamental  limitations  of  the  technique: 

1.  With  naphthalene  How  visualization,  the  edges  of  the  vortex  streaks  are  not 

1 The  term  “integration  strategies”  comes  from  the  ultimate  goal  of  computing  the  amplification 
factor  N via  equation  (3.15). 
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always  well  defined  and  can  be  obscured  by  inconsistency  in  the  concentration 
and  application  of  the  naphthalene.  Moreover,  the  wavelength  measurements 
were  obtained  under  the  conditions  of  natural  surface  roughness,  which  has 
been  shown  to  generate  a complicated  stationary  structure  containing  many 
disturbance  modes  (figure  5.28).  Consequently,  measurements  of  the  separation 
between  the  vortex  streaks  are  rather  subjective  and  open  to  the  interpretation 
of  the  observer. 

2.  As  discussed  in  section  4.5,  the  naphthalene  sublimes  quickly  near  the  leading 
edge  due  to  the  high  shear  of  the  thin,  laminar  boundary  layer.  As  a result, 
Dagenhart’s  measurements  were  confined  to  the  mid-chord  region  of  the  model, 
in  which  the  direction  of  propagation  of  the  stationary  wave  is  nearly  constant. 
In  light  of  this,  it  is  no  surprise  that  vortex  drops-outs  or  adjustments  to  the 
spanwise  spacing  were  not  observed. 

With  a high-resolution  instrumentation  traverse  and  improved  measurement  tech- 
niques, the  present  experiment  is  much  better  equipped  to  objectively  analyze  the 
wavenumber  distribution  for  the  stationary  crossflow  instability.  In  fact,  the  power 
spectra  for  the  spanwise  scans  contain  all  the  necessary  information.  By  virtue  of  the 
measurement  technique,  the  location  of  each  peak  in  the  spectrum  gives  the  spanwise 
wavelength  A,  = 2tt / [Jr  of  the  corresponding  unstable  crossflow  mode.  Tracking  this 
spectral  peak  at  each  chord  location  gives  the  wavelength  distributions  shown  in  fig- 
ure 6.226.  This  plot  shows  the  spanwise  wavelength  (A, ) of  the  fundamental  crossflow 
mode  [i.e.,  the  (0,1)  mode]  as  a function  of  x/c  for  all  experimental  data  sets  with 
12  mm  roughness  spacing.  Also  plotted  is  the  variation  of  Xz  predicted  by  linear  the- 
ory (data  provided  by  Haynes)  using  the  constant  crossflow  wavelength  constraint. 
The  value  of  ACf  - 8.5  mm  is  chosen  because  it  gives  \z  = 12  mm  at  x/c  — 0.05, 
corresponding  to  the  initial  roughness  forcing.  Clearly,  fixing  Acf  fails  to  produce  the 
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correct  spanwise  wavelength.  Instead,  the  spanwise  wavelength  remains  constant  as 
proposed  by  Mack  (1988).  It  should  be  noted,  however,  that  Xz  does  not  dramati- 
cally change  when  ACf  is  fixed,  which  explains  why  the  qualitative  measurements  of 
Dagenhart  produced  the  incorrect  conclusion. 

In  summary,  the  spanwise  growth  rate  /?f  is  zero  and  the  spanwise  wavenumber  pr 
is  constant  for  infinite-span  swept  wings. 

6.4.2  Effect  of  Acoustic  Forcing 

As  discussed  in  section  1.3,  Muller  and  Bippes  (1989),  Bippes  (1990,  1991),  and 
Bippes  et  al.  (1991)  have  shown  that  the  growth  of  traveling  crossflow  waves  depends 
strongly  on  the  freestream  turbulence  level.  However,  in  these  experiments  the  initial 
conditions  for  the  stationary  waves  came  from  the  unknown  natural  roughness  of  the 
model  surface.  What  is  not  well  understood  is  the  importance  of  freestream  distur- 
bance's (acoustic  and  vortical)  in  relation  to  surface  roughness  and  their  combined 
effect  on  transition.  At  issue  here  is  the  receptivity  of  the  boundary  layer  to  various 
types  of  disturbances,  both  in  the  freestream  and  on  the  airfoil.  While  this  problem  is 
far  from  being  solved  and  continues  to  evade  a robust  theory,  the  experimentalist  can 
provide  valuable  information  concerning  the  response  of  the  boundary  layer  under 
varying  environmental  conditions. 

With  this  in  mind,  the  present  work  also  investigates  the  interaction  between 
freestream  acoustic  disturbances  and  artificial  surface  roughness  in  relation  to  their 
effect,  on  transition  in  crossflow-dominated  flows.  The  study  ignores  the  details  of 
the  instability  growth,  concentrating  instead  on  the  transition  location  as  a global 
measure  of  the  stability  of  the  boundary  layer.  The  transition  location  is  determined 
from  naphthalene  flow  visualization  using  the  technique  described  in  section  4.5. 

The  tests  are  performed  with  the  [48|  12]  surface  roughness  distribution  (i.e.,  k = 
48  //m  roughness  with  a 12  mm  spanwise  spacing  applied  at  x/c  — 0.023).  Freestream 
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Table  6.1:  Effect  of  freestream  acoustic  forcing  on  transition  location  for  Rec  = 
2.4  x 106  and  [48|  12]  roughness. 


Frequency 
Range  [Hz] 

Transition 
Location  [ x/c, ] 

No  Sound 

0.49 

30-200 

0.49 

200-lk 

0.49 

lk-5k 

0.49 

acoustic  disturbances  are  generated  with  an  array  of  nine  speakers  mounted  in  the 
plenum  wall.  The  source  signal  contains  bandwidth-limited  random  noise.  Three  fre- 
quency ranges  are  chosen:  30  -200  Hz,  200  1000  Hz,  and  1-5  kHz  corresponding  to  the 
unstable  frequencies  for  traveling  crossflow  waves,  Tollmien-Schlichting  disturbances, 
and  secondary  instabilities,  respectively.  The  sound  pressure  level  is  held  constant 
at  100  dB  (20  /rPa  reference).  Table  6.1  displays  the  average  transition  location  for 
Rec  = 2.4  x 10fa  under  these  conditions.  The  acoustic  forcing  has  absolutely  no  effect 
on  the  transition  location.  Moreover,  no  changes  in  any  detail  of  the  characteristic 
saw-tooth  transition  pattern  are  noted.  These  results  confirm  the  observations  of 
Takagi  et  al.  (1991)  and  Radeztsky  et  al.  (1993a),  as  well  as  the  more  recent  findings 
of  Devhle  and  Bippes  (1996),  and  indicate  that  acoustic  disturbances  play  only  a 
passive  role  in  the  crossflow  instability. 


CHAPTER  7 


Conclusions 


7*1  Summary  of  Key  Results 

The  important  results  of  this  experimental  investigation  are  summarized  in  the  fol- 
lowing specific  conclusions. 

1.  Stationary  crossflow  waves  dominate  the  transition  process  on  swept  wings  in  a 
low-disturbance  environment.  Even  under  the  conditions  of  0.25  fim  rms  natural 
surface  roughness,  the  stationary  waves  cause  strong  nonlinear  distortions  of 
the  mean  streamwise  boundary  layer  and  the  jagged  transition  front  (due  to  a 
local  breakdown  caused  by  secondary  instabilities)  indicates  sensitivity  to  the 
randomly  distributed  roughness. 

2.  Leading-edge  roughness  provides  the  all-important  initial  conditions  for  the  sta- 
tionary waves.  In  the  absence  of  art  ificial  roughness,  submicron  irregularities  in 
the  natural  surface  finish  generate  nonuniform  disturbance  motion  containing 
many  modes.  Systematic  introduction  of  equally-spaced  roughness  elements 
produces  an  ideal  fundamental  mode  at  that  spacing.  The  micron-sized  ele- 
ments introduce  known  initial  conditions  without  saturating  the  initial  distur- 
bance amplitude,  thus  providing  the  necessary  database  for  comparisons  with 
both  linear  and  nonlinear  computations. 
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3.  In  contrast  to  the  results  for  isolated  roughness  elements  (Radoztsky  et  al. 
1993a),  increasing  the  height  of  the  distributed  roughness  has  surprisingly  little 
affect  on  the  transition  location.  Moreover,  transition  occurs  at  larger  x/c.  for 
distributed  roughness  than  for  an  isolated  element  of  the  same  height.  This 
behavior  is  also  observed  when  the  spacing  of  the  distributed  roughness  is  in- 
creased from  12  mm  to  36  mm,  in  which  case  transition  moves  forward  and 
the  characteristic  saw-tooth  transition  front  becomes  dramatically  more  regular 
in  span.  This  indicates  that,  through  enhanced  nonlinear  interactions  among 
multiple  modes,  the  wavenumber  content  of  the  stationary  disturbance  can  sig- 
nificantly influence  the  local  transition  characteristics  of  the  boundary  layer. 

4.  Extensive  hot-wire  measurements  document  the  detailed  structure  of  the  distor- 
tion of  the  mean  boundary  layer  caused  by  the  integrated  effect  of  the  stationary 
waves.  Boundary-layer  profiles  clearly  show  the  growth  of  the  mean-flow  distor- 
tion [(0,  0)  mode]  and  the  corresponding  development  of  the  nonlinear  second 
lobe  in  the  total  disturbance  inode-shape  profile.  This  early  nonlinear  mode 
interaction  causes  the  total  disturbance  amplitude  to  saturate  well  before  tran- 
sition. 

5.  Accurate  spectral  decomposition  provides  individual-mode  growth  rates  for  the 
fundamental  crossflow  mode  and  all  amplified  harmonics.  These  data  show 
that  the  amplitude  saturation  phenomenon  is  caused  by  the  nonlinear  growth 
of  the  harmonic  disturbances  and  their  interaction  with  the  fundamental  mode. 
The  amplitude  saturation  is  not  limited  to  the  fundamental  mode,  but  is  also 
observed  in  the  harmonic  wavelengths. 

C.  Uniformly  distributed  roughness  at  the  most  unstable  wavelength  (A^  = 12  mm) 
produces  an  ideal  fundamental  mode  and  poor  excitation  of  other  modes.  In 
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fact,  harmonic  disturbances  are  amplified  only  through  nonlinear  coupling.  In 
this  sense,  forcing  the  most  unstable  wavelength  provides  a minimum  destabi- 
lization of  the  total  disturbance. 

7.  Increasing  the  spanwise  spacing  of  the  roughness  elements  generates  multiple- 
mode stationary  waves  that  produce  an  intricate  distortion  of  the  mean  bound- 
ary layer  leading  to  earlier  transition.  Spectral  analysis  shows  that  harmonic 
disturbances  can  be  amplified  before  the  fundamental  mode,  indicating  that 
the  Fourier  decomposition  of  the  roughness  distribution  is  an  important  aspect 
of  the  receptivity  process.  This  suggests  that  the  roughness  may  excite  the 
stationary  instability  through  a wave-like  pressure  disturbance  and  not  a local 
injection  of  vorticitv. 

8.  Increasing  the  height  of  the  distributed  roughness  produces  a larger  initial  dis- 
turbance and  enhanced  nonlinear  interaction  leading  to  earlier  amplitude  satu- 
ration. In  addition,  a “two-stage”  saturation  is  observed  for  the  larger  roughness 
heights.  However,  the  total  disturbance  relaxes  downstream  to  a similar  mode 
shape  profile  and  the  transition  location  remains  largely  unaffected.  This  sug- 
gests that  once  the  nonlinear  effects  are  strong  enough  to  cause  saturation,  the 
saturation  amplitude  is  fixed. 

9.  No  wavelength  doubling  (i.e.,  modes  with  wavelengths  larger  than  the  funda- 
mental) is  observed  for  any  roughness  configuration.  On  the  other  hand,  both 
the  12  mm  and  36  mm  forcing  produce  amplified  harmonics  up  to  the  Xz  — 4 mm 
wavelength.  It  is  interesting  to  note  that  the  diameter  of  the  roughness  elements 
is  3.7  mm.  Thus,  the  measurements  show  no  evidence  of  unstable  harmonics 
with  wavelengths  smaller  than  the  roughness  diameter  (with  the  exception  of 
the  infinitesimal  A2  = 3 mm  disturbance  immediately  before  transition  for  [6 1 1 2] 
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roughness  at  Rec  = 3.2  x 106).  This  result  is  unexplained  and  illuminates  the 
importance  of  the  receptivity  process  for  stationary  crossflow  waves. 

10.  In  contrast  to  Radeztsky  et  al.  (1994),  the  early  disturbance  growth  shows 
qualitative  agreement  with  linear  theory  even  though  the  present  distortion  of 
the  mean  boundary  layer  is  an  order  of  magnitude  larger.  This  shows  that 
there  is,  in  fact,  a linear  region  for  the  stationary  crossflow  instability,  and  may 
indicate  that  the  large  roughness  used  in  the  previous  experiments  caused  the 
linear  receptivity  regime  to  be  bypassed. 

11.  Although  a region  of  linear  growth  is  observed,  the  importance  of  nonlinear 
effects  as  the  dominant  aspect  of  the  stationary  crossflow  instability  is  firmly 
established  and  incontrovertible.  This  is  reinforced  by  the  outstanding  nonlin- 
ear PSE  calculations  that  agree  remarkably  well  with  the  experimental  data. 
Together  these  results  show  that  the  important  physical  mechanism  is  not  the 
weak  ( />',  »/)  motion  of  the  stationary  wave  itself,  but  rather  its  ability  to  induce 
0(1)  //'  distortions  of  the  mean  boundary-layer  flow.  In  light  of  this,  the  failure 
of  linear  theory  to  capture  the  details  of  the  disturbance  growth  is  perhaps  no 
surprise.  Linear  theory  simply  cannot  account  for  the  distortion  of  the  basic 
state  resulting  from  the  integrated  effect  of  the  stationary  wave. 

12.  The  secondary  instability  that  leads  to  transition  appears  to  be  most  effective  in 
the  presence  of  multiple-mode  disturbances.  This  is  suggested  by  the  uniformity 
of  the  saw-tooth  transition  front  and  the  decrease  in  the  transition  Reynolds 
number  with  increased  roughness  spacing.  Because  of  these  effects,  even  linear 
stability  analyses  of  the  distorted  mean  flow  are  insufficient  to  describe  the 
transition  process. 
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13.  Traditional  transition-correlation  techniques  such  as  the  classic  cN  method  are 
ineffective  under  these  conditions.  The  nonlinear  amplitude  saturation  occurs 
well  before  transition,  thus  the  disturbance  amplitude  is  relatively  constant  over 
a large  extent  of  the  laminar  boundary  layer.  Moreover,  increasing  the  rough- 
ness height  moves  the  saturation  point  forward,  while  the  saturation  amplitude 
and  transition  location  remain  unchanged.  On  the  other  hand,  increasing  the 
roughness  height  from  6 /on  to  48  ^m  causes  a nearly  50%  decrease  in  the  distur- 
bance amplification  factor  N at  transition.  Thus,  transition  correlations  based 
on  iV-factor  criteria  are  not  possible.  In  fact,  the  strong  nonlinear  saturation 
and  its  dependence  on  initial  conditions  show  that,  in  a stationary-crossflow- 
dorninated  boundary  layer,  no  transition  prediction  is  possible  without  proper 
account  for  the  receptivity  process. 

7.2  Closing  Remarks 

Whereas  transition  correlation  is  often  based  on  semi-empirical  observations,  accurate 
transition  prediction  requires  an  intimate  understanding  of  the  physics  involved.  The 
present  investigation  provides  this  information  for  swept-wing  boundary  layers  dom- 
inated by  the  stationary  crossflow  instability.  The  importance  of  nonlinear  effects  is 
firmly  established,  and  the  ability  of  micron-sized  roughness  to  alter  the  disturbance 
growth  confirms  that  receptivity  issues  cannot  be  ignored.  This  experiment  provides 
a detailed  database  for  the  development  and  validation  of  computational  methods 
concerning  three-dimensional  boundary-layer  stability  and  transition  prediction. 
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Figure  2.2:  Side  view  of  traverse  carriage.  All  dimensions  in  millimeters 


ure  2.3:  Front  view  of  traverse  carriage.  All  dimensions  in  millimeters 
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Figure  3.7:  Schematic  of  swept  wing  boundary- layer  profile  showing  the  tangential,  crossflow,  and 
total  velocity  profile. 
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Figure  4.1:  Typical  hot-wip 
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Figure  4.2:  Typical  hot-wire  temperature  compensation.  The  scatter  in  the  data  is  caused  by  short 
sampling  times. 
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Figure  4.3:  Comparison  of  inferred  hot-wire  velocity  with  and  without  temperature  compensation. 
The  uncompensated  velocity  readings  suffer  from  errors  exceeding  15%.  With  compensation  the  error 
is  less  than  1%  over  the  entire  speed  range. 
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Figure  4.4:  Wind  tunnel  “cruise  control”  program.  The  test-section  speed,  Reynolds  number,  and 
nondimensional  frequency  are  represented  with  dual  sliders.  For  each  set,  the  lower  slider  is  the  target 
and  the  upper  slider  is  the  measured  value. 


Figure  4.7:  Typical  spanwise  traverse-alignment  scan.  The  misalignment  between  the  hot-wire  traverse 
and  the  model  is  indicated  by  the  apparent  ‘‘drift”  in  the  traverse  motion.  Scan  is  at  x/c  = 0.30. 
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Figure  5.1:  Naphthalene  flow  visualization  for  Rec  — 2.2  x 1 06,  [6|  12]  roughness. 
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Figure  5.2:  Naphthalene  flow  visualization  for  Rec  = 3.2  x 106,  [6|  12]  roughness. 


153 


Figure  5.3:  Naphthalene  flow  visualization  for  Rec  = 2.2  x 106,  [6 1 36]  roughness. 
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Figure  5.4:  Naphthalene  flow  visualization  for  Rec  — 3.2  x 106,  [6|36]  roughness. 


Figure  5.6:  Comparison  of  experimental  and  theoretical  pressure  distribution  for  Re 
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Figure  5.7:  Comparison  of  experimental  and  theoretical  pressure  distribution  for  Re, 
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Figure  5.19:  Comparison  of  experimental  and  theoretical  basic-state  boundary-layer  profile  at  x/c, 
0.35.  Rec  — 1.6  x 106.  The  dots  are  simply  identifiers  and  do  not  represent  measurement  points. 


Figure  5.21:  Comparison  of  experimental  and  theoretical  basic-state  boundary-layer  profile  at  xt 
0.45.  Rec  = 1.6  x 106.  The  dots  are  simply  identifiers  and  do  not  represent  measurement  points. 


Figure  5.23:  Comparison  of  experimental  and  theoretical  basic-state  boundary-layer  profile  at  xj 
0.55.  Rec  = 1.6  x 106.  The  dots  are  simply  identifiers  and  do  not  represent  measurement  points. 
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igure  5.31:  Effect  of  combined  (v  ,w  ) disturbance  motion  on  the  streamwise  boundary  layer  velocity 
t xjc  — 0.20.  Rec  = 2.4  x 106,  [6|  12]  roughness. 
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Figure  5.33:  Effect  of  combined  (v  . w ) disturbance  motion  on  the  streamwise  boundary  layer  velocity 
at  x/c,  = 0.45.  Rec  — 2.4  x 106,  [6|  12]  roughness. 
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5.38:  Streamwise  velocity  contours  at  x/c  — 0.05.  Re 
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Figure  5.65:  Comparison  of  experimental  and  theoretical  crossflow  mode  shape  at  x/c  = 0.10.  Rec  = 
2.4  x 106,  [6|  12]  roughness.  The  dots  are  simply  identifiers  and  do  not  represent  measurement  points 


Figure  5.66:  Comparison  of  experimental  and  theoretical  crossflow  mode  shape  at  x/c  — 0.15.  Rec  = 
2.4  x 1()6.  [6|  12]  roughness.  The  dots  arc  simply  identifiers  and  do  not  represent  measurement  points 
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Figure  5.68:  Comparison  of  experimental  and  theoretical  crossflow  mode  shape  at  x/c  = 0.25.  Rec  - 
2.4  x 106,  [6(12]  roughness.  The  dots  are  simply  identifiers  and  do  not  represent  measurement  points 


Figure  5.69:  Comparison  of  experimental  and  theoretical  crossflow  mode  shape  at  x/c  = 0.30.  Rec  = 
2.4  x 106,  [6)12]  roughness.  The  dots  are  simply  identifiers  and  do  not  represent  measurement  points 


Figure  5.70:  Comparison  of  experimental  and  theoretical  crossflow  mode  shape  at  x/c,  = 0.35.  Rec  = 
2.4  x 10°,  [6 1 1 2]  roughness.  The  dots  are  simply  identifiers  and  do  not  represent  measurement  points 


Figure  5.71:  Comparison  of  experimental  and  theoretical  crossflow  mode  shape  at  x/c  = 0.40.  Rec  = 
2.4  x 106,  [6|  12]  roughness.  The  dots  are  simply  identifiers  and  do  not  represent  measurement  points 


Figure  5.72:  Comparison  of  experimental  and  theoretical  crossflow  mode  shape  at  x/c.  — 0.45.  Rec  = 
2.4  x 106,  [6|  12]  roughness.  The  dots  are  simply  identifiers  and  do  not  represent  measurement  points 
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Figure  5.85:  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.20,  Y = 0.85  mm.  Re 
2.4  x 106,  [6|  12]  roughness. 


Figure  5.87:  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  — 0.30,  Y = 0.9  mm.  Re, 
2.4  x 106,  [6 j 1 2]  roughness. 


Figure  5.88:  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  — 0.35,  Y = 0.9  mm.  Rec 
2.4  x 106,  [6 j 1 2]  roughness. 


Figure  5.89:  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.40,  Y = 0.9  mm.  Rec 
2.4  x 10°,  [6|  12]  roughness. 


Figure  5.90:  Power  spectral  density  of  spanwise  hot-wire  scan  at  x/c  = 0.45,  Y = 0.9  mm.  Rec 
2.4  x 106,  [6j  12]  roughness. 
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Figure  6.5:  Streamwise  velocity  contours  at  x/c  = 0.25.  Rec  = 1.6  x 106,  [6|  12]  roughness. 
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Figure  6.8:  Streamwise  velocity  contours  at  x/c  = 0.40.  Rec  — 1.6  x 106,  [6|  12]  roughness. 
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Figure  6.9:  Streamwise  velocity  contours  at  x/c  = 0.45.  Rec  — 1.6  x 106,  [6 1 12]  roughness. 
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Figure  6.10:  Streamwise  velocity  contours  at  x/c  — 0.50.  Rec  — 1.6  x 106,  [6|  12]  roughness. 
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Figure  6.11:  Streaimvise  velocity  contours  at  x/c  — 0.55.  Rec  = 1.6  x 106,  [6|12]  roughness. 
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Figure  6.12:  Streamwise  velocity  contours  at  x/c.  = 0.60.  Rec  ~ 1.6  x 106,  [6|  12]  roughness. 
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Figure  6.37:  Disturbance  velocity  contours  at  x/c  — 0.15.  Rec  = 1.6  x 106,  [6 1 1 2]  roughness. 
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Figure  6.42:  Streainwise  velocity  contours  at  x/c  = 0.10.  Rec  = 3.2  x 106,  [6|12]  roughness 


6.44:  Streamwise  velocity  contours  at  x/c  — 0.20.  Re, 
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Figure  6.46:  Streamwise  velocity  contours  at  x/c.  = 0.29.  Rec  - 3.2  x 106,  [6[  12]  rough 
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6.76:  Streamwise  velocity  contours  at  x/c  = 0.05.  Re 
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6.84:  Streamwise  velocity  contours  at  x/c  = 0.45.  Re , 
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Spanwise  hot-wire  scan  at  x/c  = 0.20.  Y — 0.9  mm.  Rec  — 2.4  x 106,  [6|36]  roughu 


Spanwise  hot-wire  scan  at  x/c  = 0.30.  Y = 1.0  mm.  Rec  — 2.4  x 106,  [6|36]  rough 
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.126:  Streamwise  velocity  co 


S.128:  Streamwise  velocity  c 


.130:  Streamwise  velocity  contours  at  x/c  = 0.30.  Re 


.132:  Streamwise  velocity  contours  at  x/c  — 0.40.  Re, 
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Figure  6.138:  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a span  of  99  mm  at 
x/c  = 0.25.  Rec  = 2.4  x 106,  [18|  12]  roughness.  The  dots  indicate  the  mean  of  the  profiles. 
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Spamvise  hot-wire  scan  at  x/c  = 0.35,  V = 0.85  mm.  Rec  = 2.4  x 106,  [18(12]  roughn 


panwise  hot-wire  scan  at  x/c  = 0.45,  Y = 0.85  mm.  Rec  = 2.4  x 


5.174:  Streamwise  velocity  contours  at  xlc  = 0.05.  Re, 


.176:  Streamwise  velocity  contours  at  x/c  = 0.15.  Rec 
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velocity  contours  at  x/c  — 0.35.  Re 
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Figure  6.182:  Streamwise  velocity  contours  at  x/c  = 0.45.  Rec  — 2.4  x 106,  [48|12]  roughness. 
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Figure  6.188:  Spanwise  array  of  100  mean-flow  boundary-layer  profiles  covering  a span  of  99  mm  at 
x/c  — 0.30.  Rec  — 2.4  x 106,  [48 1 1 2]  roughness.  The  dots  indicate  the  mean  of  the  profiles. 
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Figure  6.201:  Stationary  crossflow  mode  shapes  for  Rec  = 2.4  x 106  and  [48(  12]  roughness.  The 
symbols  are  simply  identifiers  and  do  not  represent  measurement  points. 
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Figure  6.222:  Total  and  individual-mode  disturbance  amplitude  and  amplification  factor  N for  Rec  — 
2.4  x 106  and  [48|  12]  roughness.  iV-factors  are  relative  to  the  point  at  which  the  disturbance  is  first 
detected . 
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Figure  6.223:  Stationary  crossflow  mode  shapes  at  x/c  = 0.10  and  x/c  — 0.45  for  Rec  = 2.4  x 106  and 
[A;  1 12]  roughness.  The  symbols  are  simply  identifiers  and  do  not  represent  measurement  points. 
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